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Theory of Oscillation Type Viscometers IT: 
A Thin Disk’) 


By Atronso Gir AzpritTia?) and GorDON FRANK NEWELL’), Providence, R. I. 


Abstract 


The small amplitude oscillation of a thin disk of finite radius suspended in 
an ‘infinite’ fluid is considered with a view toward its use as a viscometer. The 
fluid motion is analyzed under the assumption that the boundary layer thick- 
ness of the fluid is small compared with the radius of the disk but large com- 
pared with the thickness of the disk. Formulas are obtained relating the 
frequency and decrement of oscillation to the density and viscosity of the fluid. 


1. Introduction 


The slow oscillation of a disk suspended in an ‘infinite’ fluid by a torsion 
wire has been used by several investigators as a means of determining the 
viscosity of the fluid [1-5]*). The simplest approximate relation between the 
decrement of the free oscillation of the disk and the viscosity of the fluid has 
‘been obtained by assuming that the effects of the edge of the disk are negli- 
gible, i.e. by assuming that the drag per unit area exerted on the disk by the 
viscous fluid is the same as that which would be exerted on a disk of infinite 
radius oscillating at the same rate. For most practical applications, this ap- 
“proximation is, however, too crude for accurate viscosity measurements. 
Measurements of the decrement can be made with an accuracy ranging from 
about 0-1 to 1% whereas the error in the approximate theory may be any- 
where from a few percent-to 100% in typical applications. 
Several attempts have been made to improve the theory [1-5] by making 


corrections for the edge but in all cases these corrections have been of a semi- 


ed States Air Force through the Air Force Office 
of Scientific Research of the Air Research and Development Command, under contract No. AF 
18 (600) 1548. Reproduction in whole or in part is permitted for any purpose of the United States 
Government. Part of the work was done while one of the authors (G. F. N.) was being supported by 
a grant from the Alfred P. Sloan Foundation. 
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4) Numbers in brackets refer to References, Pp 
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empirical type and, though qualitatively correct, do not give predictions of _ 
well defined accuracy. The difficulty is primarily a mathematical one. Though 
the equations of motion for the fluid and for the disk can be accurately formu- _ 
lated they can not be accurately solved (except in certain limiting cases). 

It is the purpose of the following paper to analyze the equations describing 
the fluid motion and derive a relation between the decrement and viscosity 
assuming that the thickness of the disk is small compared with the boundary 
layer thickness of the fluid but the radius of the disk is large compared with 
the boundary layer thickness. This will extend the theory of the infinite disk 
mentioned earlier and include an approximate influence of the edge. It will 
however be obtained in such a way that the accuracy of the formulae and their 
range of applicability are both well defined. 

A general formulation of the theory of oscillation type viscometers along 
with a list of notation is given in Part I [6]. We shall only very briefly outline 
here the formulae of Part I which are relevant to the present problem. 

The mathematical difficulties arise mainly from trying to find solutions, 
w(é, 4, s) of the equation, equations (13), (14) of Part I, 


0?w 3 Ow 0?w 
Oe cE OR ay @) 
subject to the boundary conditions 
w(é, 7, s) = 1 on the surface of the disk. (2) 


We need not necessarily find the complete solution of equation (1) but we must 

obtain enough information regarding it to evaluate the integral, equation (16) 

of Part I, Bs (Pan Ow(E 
ee Ne peck a eS 

D(s) = 25 /| ga SO TS! do. (3) 


ve 


A 


After any possible transients have decayed, the angular displacement «(z) 
of the disk is a decaying oscillation of the form, equation (18) of Part I, 


a(t) = AeA?" cos(w vt + yp): (4) 


The motion of the disk and that of the fluid are related in that the values of A 


and w are obtained by solving the transcendental equation, equation (40) of 
Parr: 
(S + Ao)? + 1+ D(S) =0 (5) 


for a pair of complex conjugate roots somewhere near S ~ + 7 and evaluating 
« and A (real numbers) from the complex equation 


S=(ti-A)o. (6) 
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The fact that D(s) involves the viscosity » means that equations (5) and (6) 
define a relation between w, A and v which would be used to find » from 
measured values of w or A or both. 
Once D(s) has been found, the determination of » from known values of 
A and/or w can be done by a method of successive approximations. Although 
this step may be somewhat tedious it is straight-forward and presents no 
serious difficulties. The most serious problem is that of finding an accurate but 
manageable expression for D(s) in terms of the various physical parameters 
describing the fluid and the disk. It is this problem which will be of primary 
concern here. 


Figure 1 
Coordinates of disk. 


The geometry of the disk is shown in Figure 1. We let 2 % and &, be the 
thickness and the radius of the disk respectively, measured in units of the 
boundary layer thickness 6. This differs from the notation for the cup in Part I 
only in that these same quantities now refer to outside dimensions of a disk 
rather than the inside dimensions of a cup. The condition that the thickness 
of the disk be small compared with the boundary layer thickness which in 
turn is small compared with the radius of the disk means that 


NM K<1<&.- (7) 


Under these conditions, the equations for w must describe a situation 
consistent with the known fact that as one leaves the disk, the amplitude of 
the fluid motion decreases and |w| is appreciably less than 1 at a distance 
from the disk of order 6 or a distance 1 measured in units of 6. The surfaces 
of constant w are nearly parallel with the top and bottom surfaces of the disk 
along most of these surfaces and deviate from planes only within a distance 
of the order of one boundary layer thickness from the edge. . . 

To obtain a first approximation to D(s), one is justified in neglecting the 
influence of the edge because the integral in equation (3) extends over the 
entire surface, of which the vicinity of the edge is only a small part. Large 
errors in 0w/0n only near the edge will give only relatively small errors in D(s). 
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We, therefore, neglect the dependence of w(&, 7, s) upon & for € < & and 
choose®) 


exp[— s¥2 (7 — )] for 7 > 1% §&< Eo , | 
Wy = | exp[+ s¥4 (4 + mo)] for < —m, € < £9; (8) 
(0) : for é => & ’ | 


which is a solution of equation (1) except for a discontinuity at § = &). Substi- 
tution of w, into equation (3) gives 


9 d8x Sf 58? 


I 


D(s) ~ Dols) = =m si, (9) 
in which m = 0 6° 4 2/I is the similarity parameter of [1]. 

The nature of this approximation is such that equation (8) is an accurate 
representation of w for &; — > 1, 1. e. sufficiently far from the edge, in fact 
the error, w — Wy, decreases exponentially as one moves inward from the edge. 
The error in Dy, arises almost entirely from making a very bad approximation 
to w in the immediate vicinity of the edge. 

The fact that the error is confined to the vicinity of the edge can be used 
to considerable advantage. To an observer sitting near the edge and viewing 
the influence of the edge over distances of the order of one boundary layer 
thickness, the center of the disk is very far away and he sees the edge moving 
with nearly a linear oscillation in a direction perpendicular to the radius vector 
from the origin rather than with a torsional oscillation. 

The mathematical consequences of this can best be seen by choosing a 
coordinate system with origin at the edge. We let 


= EF, = E, (10) 
and rewrite equation (1) in the form 
0?w 3 Ow 0?w 
ime Cine Or ea oe 1p) 


We anticipate that 0?w/dx?, 0?w/d7? and dw/0x are all of order 1 over most 
of the region near the edge. The term (& — x)-! dw/0x is, however, of order 
é' <1. Far from the edge (x > 0) both 02w/0x? and dw/dx are small in 
accordance with the p-evious approximation. For a next approximation we, 
therefore, can neglect the first derivative term in equation (11) and define 


5) Although we are interested in solutions of equation (1) for complex values of s and par- 
ticularly for s ~ i, it is convenient in deriving formulas to consider s as a real positive number and 
s!/2 as its positive square root. Formulae for complex values of s can always be obtained by appro- 
priate analytic continuation into the complex plane. : 
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W(x, yn, Ss) to be a solution of the equation 


0?w 0?2w 
ea TiN = SW, . (12) 


The solution of equation (12) with w,=1 on A is still a very difficult 
problem, but it is possible to find an exact solution for the special case 7) = 0. 
This is done in section 2. In sections 3 and 4 we will obtain corrections for the 
curvature of the edge, the error due to neglecting the first derivative term in 
equation (13) and in section 5 we will obtain a correction for a non-zero but 
small value of 7p. 


2. Edge Correction for Infinite Radius and Zero Thickness 


We consider here the solution of equation (12) near the edge of a disk of 
very large (infinite) radius and zero thickness, i. e. the solution for which w,=1 
on a semi-infinite line (Figure 2). The solution is obtained by the Wiener Hopf 
technique’). 


ft 
47) 
| 
en = 
on : W,=1 
ee eee aa PS 5 I eee ee he 
x —», 
| 
| 
| 
Figure 2 


Geometry of thin disk. 


We first define the Fourier transform of w, with respect to x as 
lee) 
w*(p, , 8) = | Ax e428? w(x, 7, S) . (13) 
It is clear that this integral will not exist for = 0 or for Imagp > 0 because 


we expect that w, > Wy = exp(—s¥/?| |) for x > + oo. We assume, however, 
that the integral exists in some region of the lower half plane, say for 


Imagp = — ¢« with e > 0. The inversion formula for equation (13) in this case 
is given by ee 
il C , 
W,(%, 4, 8) = or / dp e+tP@ w*(p, n, S) . (14) 


6) Most of the work described in this section was taken from unpublished notes written by 
G. F. Carrier. A brief description of this is contained in [7]. The results were also recalculated by 
L. N. PERSEN. 
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The Fourier transform of equation (14) gives the differential equation 


0? wi 

“Fr ee ee 
and the most general solution of this equation that remains finite for | 7 | > co © 
is 


wt = A(b, s) exp[— (s+ 7" |y|] (15) 


in which A(f, s) is some arbitrary function to be determined so that w, satis- 
fies the boundary conditions at 7 = 0. 
In addition to the boundary condition 


2,(%,0,s)=1 forx>0, (16) 
_ the condition 
aoe I gt fare (16a) 
On n=0 


is implied by the symmetry of w, with respect to reflections across the line 
7 = 0. The boundary conditions tell us nothing about w,(x, 0, s) for x < 0 or 
0w,/07|,—o for x > 0, in fact the latter quantity is what we particularly wish 
to find in order to evaluate D(s). 


Suppose we define ; 


eee [a e-*P2 w,(x, 0, s) (17) 


—oo 


to be the Fourier transform of the (unknown) function w, along the semi- 
infinite line x < 0 and 


+00 ie.) 
ine OW %,-7, 5) i spe OWil%, 7, S) 
Wp, y= [avensoe 20%) Faggian ale 
(p, s) : ae oF 0. dx e au vals. (18) 


to be the Fourier transform of the (also unknown) function dw,/dy along the 
top side of the infinite line » = 0, (7 = 0,), or equivalently along the semi- 
infinite line x > 0. 

These Fourier transforms can also be found in terms of A(p, s) through 
equations (15) and (16). We thus obtain two relations among the three un- 
knowns /,h and A, namely 


Alps) = 16.8) — 5, = (6+ PAG) =AG,3). (19) 


Although these are the only equations we can obtain from the boundary 
conditions, they do not describe all the information we know about these 
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functions. If we insist that w,(x, 0, s) > 0 faster than e** as x > — oo, for some 
é> 0, then f(p, s) defined by equation (17), is analytic in # for all # with 
Imag? = — «. On the other hand, h(, s), is defined and analytic for Imag p<0. 
We also insist that /(p, s) > 0 at least as fast as order |f|~1 for Imagp > + oo 
and h(p, s) > 0 for Imagp > — oo. 

_ These requirements on h and / along with the one equation relating / and 
f obtained by eliminating A from equation (19), 


hp, s) = —(s + 9°)" [#6, 8) — |], (20) 


uniquely determine both / and 7. The fact that /(p,s) is analytic for 
Imagf = — « means that the analytic continuation of /(p, s) into the upper 
half plane must have a simple pole at the origin and a branch point at p =7/s. 
The fact that h(p, s) is analytic in the lower half plane means that the analytic 
continuation of f(p, s) has a branch point at = —iV/s. These are the only 
singularities that equation (20) will allow. Equivalently, we could rearrange 
the terms of equation (20) into the form 


St ethan = F Wa ip) — 5H] — VE — ip) I, 5) - (200) 


We see that the right-hand side of equation (20a) is analytic and bounded 
at least for Imagp = —e. The left-hand side is analytic and bounded for 
Imag < 0 and vanishes for Imag~ > — oo. Although equation (20a), as 
derived above, holds only in the strip 0 > Imagp = —«, the analytic con- 
tinuation is defined everywhere. We conclude that both sides of equation (20a) _ 
must be analytic and bounded everywhere in the p-plane and vanish for 
Imag ~p > — co. However, the only function which is analytic and bounded 
everywhere and vanishes at oo is the null function and therefore we conclude 
‘that both sides of equation (20a) must be identically zero, Lae: 


Ae ahem We ape (21a) 
and 
i((vVs— 4 p)t? — s*4] 
t(, s) ae (Vs—i py? p : (21b) 


Knowing h(p, s) explicitly, we can now work our way backwards and 
obtain A(f, s) from equation (19), wf from equation (15) and finally w,(x, 7, s) 
from equation (14), 

+00—1E 
1 s/t exp[i p » — (s + P*)"? In] 
w4(%, 4, 8) = Bet / d b (Vs — 1 py? : (22) 
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The integrand of equation (22) has a pole at # = 0 and branch points at 
p= +iys. If x > 0, we can close the contour with a semi-circle as shown — 
by the solid curve in Figure 3a, and fer x<0 by a semi-circle shown by 


p=/V¥s 
branch 
point 


\ p--i¥s” / 
e branch Zp 
WK point 7 
~ zt x a 

a 
X) p=iVs 
0 
Pa 
X) p=-iys 
. y p 
b 41 
lt 
|| 
Figure 3 


Deformed contour integrals in p-plane. 


the broken curve. The two contours can then be deformed into those shown 
in Figure 3b. For x > 0, there are two contributions to w,, one from the pole 
at =O and one from the branch cut in the upper half-plane. For x < 0, 
there is only the branch cut in the lower half-plane. The contribution to w, 
from the pole at p = 0 is quite interesting for it gives 


| exp(— ys ||) for = 0: 
Ds f i Meas (23) 
or x ; 
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which we recognize as the same function wy, given by equation (8) of the 
previous section for 7) = 0. 

The contribution to w, from the two branch cuts describes the distortion 
of w, caused by the edge. By changing the variables of integration so that 
p=iysy for x > 0 and p = —i/sy for x < 0, we can obtain integrals along 
the real line and so produce the relation 


vi yy Syne sin [Vs (y2 — 1)!? ||] exp(— y's y x) 


= : ( ) 


co d - ‘ 
i y (y = cos [js (y? — 1)¥# | 7 I] exp(+/s y x) 
; forx< 0. 


We are not able to evaluate this integral explicitly but we can evaluate 
0w,/0n at 7 = 0 which is a measure of the drag force per unit area exerted by 
the fluid on the disk and which appears in the integrand of equation (3). We 
find that 


eed = sti + (xs) exp(— v/s) —eric(xys)"?}, (25) 
n=0, 
in which erfc is the complementary error function. The first term, — Vs, is 


the contribution from w, and the remaining terms describe the influence of the 
edge. For small x, equation (25) gives a square root singularity, 


Ow, Pay earth! (x 4) 3 for | x /s | < ip (25a) 
On \n= 0, 
whereas for large x, the influence of the edge decreases at least exponentially, 


Ow, 
On 


mw ya{i 5 28 ys exp(— x /3)} for |xy5|S>1. — (25b) 


n=0, 


Our next approximation to D(s), D,(s) is obtained by substituting w, into 


equation (3), 


So 
: 


Ow, OW 
D,(s) = Dols) —4._ms eal. dx (& — x) | es Es eke (26) 


0 


We have just seen that the integrand in equation (26) decreases very rapidly 
with increasing x and by the time x reaches E,, it is dominated by a factor 
exp (— &, /s). We have assumed that Leys (= however, and we are in on 
position to account for errors of this order of magnitude. We cannot Hope to 
evaluate D(s) more accurately than to within some finite power of ny , hever 
to within an exponentially small error. The upper limit of integration &, in 


| 
i 
| 


. 
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equation (26) can therefore be replaced by co with no significant loss of 
accuracy. 

D,(s) can be found by substituting equation (25) into (26) but a simpler and 
more direct scheme is to express D,(s) directly in terms of h(p, s) and its 
derivatives. From equations (18), (21a) and (23), we see that 


r § fOay Ome _. a Of is 
[ane On On fee apn ile. sha p oe (27) 
: air in go-to ea 
(n+ 1) 2n+1 sni2 


This gives the values of the integrals in equation (26) if we let & > oo in the 
integration limit and expand (£) — x)? in powers of x. We thus obtain 


Dy(s) ~ ms [1 +2 (Ey3)4 — > EoV8* + > EoV9)~* — ze EoV9)-4]- (28) 


The error in equation (28) is exponentially small [proportional to 
exp(—é, )/s)] but the error resulting from approximating D(s) by D,(s) is of 
order (&|/s)-?. Only the first two terms of equation (28) are significant, 
therefore. 


3. Curvature Correction, Zero Thickness 


Assuming still that the thickness of the disk is zero, the chief source of 
error in the analysis of the previous section arises from the use of equation (12) 
instead of equation (11); we have neglected the curvature of the edge by using 
effectively a rectangular coordinate system rather than a cylindrical one. 

Since the first derivative term of equation (11) is important only near the 
edge at most, the next approximation is obtained by replacing & — x by & 
in this term and choosing w,(x, s) as the solution of 


O*w, 3 OW, | Owe 
Ox? en) Oa On? 


SW, (29) 


with the same boundary conditions as w,, equation (16). 

The solution of equation (29) can be obtained by almost exactly the same 
procedure as used in the previous section. The only difference up. to and 
including equation (20) is that (p? + s) is replaced by (p2+ 376 &'+ s). The 
latter expression has roots at 7% and — i 8 with 


= 9 \1/2 3 Mea; 9 \1/2 3 
a=(stae) —ge, B= (sta_)" + aR (30) 


instead of at the points + 7 /s. The two roots are on opposite sides of the real 


: 
: 


\ 


» § 


i 
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axis, however, and a separation similar to equation (20a) is possible, namely 


t h(p, s a ‘ 
ate Me = Fi — P-G-i AIG). 1) 


Since the right-hand side is analytic in the upper half plane and the left-hand 
side analytic in the lower half plane, we conclude that 


Fan ae BEM i p)ile (2) 
und 
ores SC —1 1/2 
(x, 7, s) = ae J dp expl? p 7 e — et + s)!/? | n]] es 


This is the same type of expression as equation (22) and from it we obtain 


OW, 
On n=0 


pete Vs {1 + (x «)-1? exp (— x x) — erfc (x x) 1} . (34) 


J,(s) is finally evaluated in the same way as D,(s) and differs from D,(s), 
quation (28), only in that (& «) replaces (& |/s). If we expand &, « in powers 
f (€) /s)1, we obtain for D,(s) the expansion 


Da(s) ~ ms {1 +2 EyV3)-1+ > Eo)? — | EvV8) 
(35) 


2 eye}. | 


onsidered as an approximation to D(s), the first three terms of D,(s) are 
ignificant instead of just two terms as in equation (28). 


4. Variational Method 


The procedure used in the previous sections can undoubtedly be extended 
urther for 7) = 0 but with increasing complexity at each stage and no satis- 
actory extension to 7) + 0 appears obvious. We now shift the direction of 
ttack to find a variational method. The chief advantage of this is that D(s) 
ill be expressed as an integral which has a stationary value with respect to 
mall changes in the function w(&, 7, s). D(s) will then be relatively insensitive 
) small errors in w and so we will be able to make better use of the approxima- 
ons to w derived in the previous sections. 

We begin with equations (1) and (2). Multiply equation (1) by €* w and 
itegrate over the volume V outside the disk. 


feceo[Be + E+ Be e]-0 
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with dr = 22 dé dy a cylindrically symmetric volume element of V. If we 
integrate this by parts (or use the divergence theorem) and use the fact that 
w = 1 on A, equation (2), and w = 0 at oo, we obtain the equation | 


face [(iey + (SE) eee] =[oome 6g 
V A 
with do = 27 dé. 


This equation is valid subject to the condition that w is a solution of 
equations (1) and (2). Comparing the right-hand side with equation (3) we 


see that - 
Diarra tet esl oe ties 6 
V 


If we should ask under what conditions this integral has a stationary value 
with respect to changes in w, we will find that the Euler equations give back 
equation (1). In other words the value of D(s) given by equation (37) is insen- 
sitive to small variations in w from values of w which satisfy equations (1) 
and (2). 

We can integrate equation (37) by parts, more or less reversing the proce- 
dure leading to this equation, assuming w = 1 on A, but without necessarily 
assuming that w is a solution of equation (1) and arrive at the form 


5 ie 2 2 
D(s) = ae s bares oe soy dt 2 w Ee " : : se = nt —s vl » foe 
A i 

If w is a solution of equation (1), the integral over V vanishes, but if w is 
only an approximation to the solution of equation (1), the integral over V 
describes a correction to the first term that considerably reduces the errot 
resulting from taking only the first term, as was done in the previous sections. 

As an illustration, we can apply this directly to the approximations w ~ w, 
and w ~ w, for 7) = 0. Using the fact that w, and wy, satisfy equations (12) 
and (29) and that the first term of equation (38) is D, and D, respectively 
equation (38) can be written in the form 


D(s) ~ Dy(s) + 6 ms &5* | dx (€ —»)* | dy w, 9, (39a 
D(s) ~ Dals) +6 ms 5% | dx x (& — x)° [ dn w, ie (39b 
: 


The upper limit on the x integral can be made infinite with no significant los 
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f accuracy and the functions w, and w,, from equations (22) und (33) used. 
he resulting four-fold integral can be evaluated by reversing the order of 
tegration. The procedure is straight-forward though somewhat tedious and 
gives the results 


D(s) ~ m s3/? [1 + 2 (€ st4\-1 + = (Eq su2)—2 + . (E, s¥2)-3 


i | (40a) 
sensed, | 
D(s) ~ m s3!2 [1 + 2 (E su2)41 + > (E, sl2)-2 — = (Ey su2)-3 | 
40b 
+ > Eo sty + ---] abe 
2 0 


The interesting feature of equation (40a), which is based upon w ~ w,, is 
he fact that it gives the first three terms of D(s) correctly, as seen by comparing 
t with equation (35) whereas the original evaluation of D,(s), equation (28), 
fave only two terms correctly. Although D,(s) was correct to only three terms 
ve expect that equation (40b) which is also based upon w ~ wy is correct to 
ive terms. 


5. Thickness Correction 


By far the most difficult correction is that for a small but non-zero thick- 
ess of the disk, 7) + 0. This perturbation in geometry causes a change in the 
nalytic properties of w and causes singularities to appear at the two right angle 
urns of the edge instead of at the one 180° turn for 7) = 0. 

There are several ways in which we might proceed but most of the obvious 
‘ays lead to very complicated mathematical expressions which are almost 
npossible to handle. The scheme which would seem to give the greatest 
ecuracy with the least effort (but which itself is by no means trivial) is the 
lowing. 


\) We consider the problem for 7 + 0 together with that for 7) = 0. We denote 
y x’, n’ the spatial coordinates when 7 = 0 and by x, 77 the spatial coordinates 
hen 7) + 0. By means of a conformal mapping, we find a coordinate trans- 
rmation x'(x, 7), 7’(x, 7) which maps the boundary of the disk for 1) = 0 
to the boundary of a disk for 7) + 0. 


3) We take the function w,(x’, ’) for 7 =0 obtained in section 2 and 
ap it into the (x, 7)-space to obtain a function 


w(x, 4) = w,(x'(x,), n'(x, 0), (41) 


e 


: 
7 
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which will satisfy the boundary conditions, equation (2) exactly and which for 
Ny <1 will be an approximate solution of equation (1). t 


(C) We substitute w{(x, 4) into the variational form for D(s) and calculate 
the lowest order terms in %p. 


Figure 4 


Conformal mapping. 


(A) to find the transformation x'(x, 7), 7'(x,) is relatively simple. We 
think of « and 7 as real and imaginary parts of a complex variable 


Z=x+19 . (42a) 
and similarly 


z= x +47). (42b) 


The transformation 2’ = z'(z) can, if necessary, be found in two steps. One can 
first find a transformation W = W(z) that maps the geometry of Figure 4a 
into Figure 4b with the boundary of the disk mapping into the line V = 0. 


This is a standard Schwarz-Christoffel transformation determined by the dif- 
ferential equation 
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One can also map Figure 4c into 4b by a Schwarz-Christoffel transformation 
the same transformation with Ho = 0 


’ 


dz’ a. 
aw 


By combining these transformation we see that 


BE ee 2, 1/2 
tae eg) oe) 
und 
p= 0 2''2 (2 nya) — 2’)? + 2 7, ati sin-a( 27 alo : 
: (43b) 
, , . 1/2 1/2 
22 aah ing 28 all 2)" (2+ £2) 
The form of equation (43) for | z’| > 2 mp/7 is of particular interest. It gives 
er Sie Ves Lo Bays gp Oe 1 
eae? ee dau) = log > Ose (43c) 
‘rom this we see that for x’ > 7’ = O(1) 
, n Y] 240 Rien oe i eh Eotiel O 
x wx + S + = log one y ~ (4 — %) + if Valais ta 


‘he radii of the disk for 7) = 0 and for 7) + 0 as defined by this mapping are, 
herefore, not equal but are related by the equation 


Ae 
Sep aaa Ss “le 4 ae. hp ek ag (45) 


No 


B) An explicit expression for w}(x, 7) is difficult to obtain although it is 
nplicitly defined through equations (22), (41), (42) and (43). We can show 
owever that w satisfies the differential equation 


KA 4 9/2 


dz’ \2 , i 1/2 Rat 
Aged (4 — 2 n/a)? + 7] St bso) 


dz 


02wt 0? wy 
Oxn* On? 


S 


hereas we would prefer to have a solution, w(x, ), of equation (12). t. 

If we compare equations (12) and (46) we see that the only difference is im the 
efficients of s w, or s w;. For x’ > 2 yo/z, the difference in the coefficients is of 
der 4 1/a x <1, at most, as compared with other terms in the equation which 
e of order 1. At x’ = 2 m/z, 7’ = 0, the coefficient in equation (46) is infinite 
it one can show that for x’ ~ 27,/z and 7’ = 0, the terms on the left side of 
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equation (46) are of order (x’ — 2/7)? whereas the term on the right-hand 
side is only of order (x — 2n9/7)~1. Throughout most of the region where Me ang 
n’ are of order 7, the terms on the left side of equation (46) are of order 75 
whereas the terms on the right-hand side are of order 1. 

We conclude from this that the difference between equations (12) and (46) 
is small in the sense that, for 7) <1, the fractional difference in the second 
derivatives of w/ and w, is of order 7. This in turn guarantees that the first 
derivatives of wi and w, and the functions themselves will be nearly equal 
everywhere. 

A simpler approximation to w,(x, 1) than w;(x,7) would be w(x, 7) itself 
for 7) = 0. One can show that the difference between w,(x,7) for 7) = 0 and 
w(x, 7) for 0 < y) <1 is also small (of order 7) at most places). The advan- 
tage of w}(x,7), however, is that it has the correct type of singularities near the 
corners of the disk, gives a good approximation not only to w,(x,7) but also to 
its first and second derivatives (even near the corner) and, most important, it 
satisfies the boundary conditions exactly. 

Any improvement of w; by perturbation methods for example is certain to 
be very difficult since w{ is already quite awkward to handle. Our best chance ~ 
of evaluating the effects of a small 7 is to use w} in a variational calculation 
of D(s). 


(C) The evaluation of D(s) from equation (37) is the most difficult step. We 
first transform the volume integral back to the z’-space so that 


=e0 
D(s) =2ms éo* | dr! | dx’ [fg — x(x’, 7’) ]? | (47) 


(Ge) + Ge) +9 ae 


| 


If we substitute w,(x’, 7’) for w(x’, 7’) we can expect to determine correctly 
only the lowest order term in 7% oe in addition only the first two terms in the 
expansion of this in powers of 1, i.e. we can determine the first correction for 
a small thickness but infinite radius and the first correction to this due to a 
large but finite radius. We will concentrate on finding just these two things by 
discarding any parts of equation (47) that are of lower order and subtracting 
away the contribution for 7) = 0 which we have already evaluated. 

We saw from equations (23) and (24) that w, can be decomposed into two 
parts w, and w, — wp, the first of which has a long range in x > 0 and the other 
is confined to the vicinity 4 the edge. It is also convenient to separate D(s) 
into two parts D( and D( the first arising only from w, and the latter con- 
ieee terms arising from cross-products of w, — w, and w, and from (wy — Wo) 
alone 
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The easiest way to evaluate D( is to use a mixed coordinate system (x, 7’) 
so that 


m= 2 ms &,* |/ ax dy! | — x]§ s exp[— 2ys | 7’ |] [1 - ee ‘ Ss 
de | 


alsa 


(48) 


in which the limits of integration are defined by the conditions x < € and 
n(x, 7) > 0. 
Using equation (43a) we deduce that 


I 


The region of the (x, 7’)-plane, | z| = O(7), in which this relation is not valid, 
however, has an area of order 75 and, despite some singularities in the integrand, 
the contribution to D(®) from this area will certainly be of order less than 1p. 
We will, therefore, replace this factor in equation (48) by 2 everywhere. 

Equation (48) now depends upon 7 only through the limit of integration 
x'(x, 4’) > 0 which for 7) = 0 would read x > 0. From equation (43c), we see 
that x’(x, 7’) > 0 corresponds to 


dz 
dz’ 


‘ ee =24+ 0( 2) for [2 S>7,: 


: No 2|9'| x ge 
a> — 2 Pe IPE aA for | 4’ | > no - (49) 
Again the contribution to equation (48) from the region | z’ | = O(mp) is negli- 
vible and we can use equation (49) for all | 7’|. 

If we subtract from equation (48) the contribution to D(s) for 7) = 0 and 
let D and D\) denote the first two terms in an expansion in powers of £5 + of 
the part of D® that is proportional to 7, we find that 


PyOh ee SNe [1 y + log(=)| » DP= 05 (50) 


1 Tica a 


n which y ~ 0:5772 is the Euler-Mascheroni constant. An interesting feature 
of equation (50) is the fact that for 7) <1, the main contribution is of order 
7 logy. As a practical matter, however, 7%) is never so small that we can con- 
ider log 7, as a large quantity and it is usually better to think of log 7) as being 
of order 1. 

That part of the integrand in equation (47) that gives D\ decays exponenti- 
ly going away from the edge. Since we neglect errors of the type exp (—/'s &5), 
here is no necessity of restricting the range of integration and so we let 
9 > co in equation (47). The only way that 1 influences the value of D@® 
ow is through the fact that x(x’, 7’) and |dz/dz’ |? depend upon mp. If we denote 
yy D”) and D!) the first two terms in an expansion in powers of &5 ' of the part 


AMP IXa/8 


| 
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of D® that is proportional to %») we obtain from equation (47) the expressions: 


+00 +00 ! : ean : | 
DM. 2 ms2ty eee (wh — wi) (a — 1) | 
and ve 
+00 Gee 
D® 26 ms &* | dx’ | dr 
«fe (a + Cay eibree: cay 


Except in a small region near the origin, 


| dz |2 =—2 9 ¥" 
| da’| a (x72 + 7’2) 03} 
and equation (51) can be written in the form 
pase SG Jo je i! : oS? a (54) 
Tt ca x + 7’? 


The integrand of equation (54) is a function only of the parameter s but not 
No Or &). Even s can be eliminated from the integrand by rescaling the coordi- 
nates by //s in which case one can see immediately that D! is proportional to 
ms! y/E, with a rather complicated real integral as a constant of proportio- 
nality. 

In view of the fact that only an integral representation of w, is known, 
equation (54) involves a four-fold integration. Some integrations were done 
analytically but other parts were done numerically to obtain 


De® =m 532 10 0-73 . (55) 
SO 


Equation (52) can be greatly simplified by some rather elaborate partial 
integrations. The object in this is to eliminate all first derivatives of w, from the 
double integral at the expense of introducing second derivatives. At the same 
time, one does this in such a way that first derivatives of x’ — x are also 
eliminated in favor of second derivatives. The second derivatives of w, are then 
removed by using equation (12) while the second derivatives of x’ — x are 
removed by using the fact that x’ — x is the real part of an analytic function 
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and therefore a solution of LapLace’s equation. This leads us to a form | 


mo) = a pte fe ( x’ — x) jek = wo) 
a 1’ =0, 


On’ 
— ila [ae eee |e (56) 
0 
— fax fay x' s (wi — we) (fe i). 
100) 0 


By using some of the same types of approximations used above in evaluating 
D), the first and second integrals in equation (56) can be evaluated analytically 
with the first integral giving the term containing 7, logy) that appears below. 
The third integral is similar to the one which appears in equation (51). It also 
can be reduced to a constant multiple integral, part of which was evaluated 
numerically. We finally obtain 


Do) Sees Te. [log (= =) +16]. (57) 


6. Summary 


By combining the terms from equations (40b), (50), (55) and (57); we arrive 
at our final expression for D(s), 


: 2 3 3 | 3 
D(s) = m 8°? ‘ + E, sue 2 (E, st/2)2 4 (E, si/2)8 2 (& st/2)4 
- see [— log (779 st?) + 2-14] (58) 


+ ere say (— log (ro s*2) + 2-0) + - a 

As a consequence of equation (7), the power series in (& s"’*)-1 of the first line 
contains successively smaller terms and the error resulting from terminating 
this series after five terms is of order (£,s"/*)—*. If we consider the slowly vary- 
ing function log 7 to be of order 1, the second line contains the first two terms 
in a similar power series expansion of the first thickness correction. The error 
resulting from terminating this series is of order (79/&) (Fo sV?)-2. Finally there 
are errors due to neglecting higher power terms in 7. This error is expected to 
be of order (7/&>) (779 s"/”). 

If we denote the half thickness (in any desirable physical units of length) by 
h and the radius of the disk by R then h = 7) 6 and Rk = &, 6 with h and R meas- 
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ured in the same units as the boundary layer thickness 6. Equation (7) then 


implies 


h<Zs<R. 3 (59) 


This can be true only if #/R is very small, a severe restriction on the geometry 
of the disk. 

Once a sufficiently thin disk has been chosen, the decision as to whether or 
not the value of 6 is suitable depends, among other things, on.the viscosity of 
the fluid. Since the same disk would normally be used to measure viscosities of 
a number of gases, one question that will arise is: given the values of A and Rk 
(with h/R <1), for what values or range of values of 6 is equation (58) most 
accurate? If 6 is too large, then the error term of order | & s/?—>|x 6° will 
dominate and if 6 is too small, the term of order (7%o/&o) | 9 s¥/? | « 6+ will 
dominate. The probable error will be a minimum when these two errors plus the 
error of order (%/&) |&  s"/”|? are all comparable. This situation obtains for 


w= [(g)"]> m= 965" (60) 


in which case all errors are of order 5° = O[(A/R)°*). 

Since only the order of magnitude of the errors are known, it is impossible to 
put reliable bounds on the error although one can make a reasonable guess. 
There are three error terms in question and they may very well all have the 
same sign so as to add rather than cancel. Furthermore most coefficients in 
equation (58) are larger than one. A reasonable guess of the probable error un- 
der these most favorable conditions, equation (60), would be 5 (h/R)*? and a 
conservative estimate perhaps 10 (h/R)°*. Thus for h/R = 1/40, we can expect 
errors of about 1% and for h/R = 1/160, about 0-1% in the range of 6 given 
by equation (60). ' 

Regarding the possibility of improving the accuracy of equation (58) by 
computing more terms, it seems likely that one could, by some way or another, 
evaluate more terms of the first series in powers of (&)s'/?)-1 and one could 
certainly improve the accuracy of the second series (first order terms in 7) by 
using w, in the variational scheme of section 5 instead of w,. Since the evalua- 
tion of the last few terms of equation (58) has already become quite tedious, 
any extension of this type would not appear to be very rewarding. More serious 
than this, however, is the fact that we see no practical way at all of evaluating 
terms proportional to 75. We will discuss at the end of Part IV some improve- 
ments of a different sort, however. 

In practical applications, equation (58) is to be used in conjunction with 
equations (5) and (6) to determine the viscosity of the fluid from measured 
values of m and A. Equation (5) is complex and so gives two real equations 
which could potentially be used to evaluate two unknowns but in practice, one 
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usually considers only one equation, the imaginary part of equation (5), 
+ 2 (4 — Ay) =ImD {(+iw—dAo)}. (61) 


The imaginary part of D is easily obtained by writing s = (i — A) in the 
form |s| exp(+ 7) with 0 < »< zs0 that 


si =|s|cos(jy) +2 |s| sin (jg) 
and 


1 I, oe 
logs? = 5 log|s | See 


The viscosity appears in D through m, &, and 7 all of which are proportional 
to v'/? and it can be determined from equation (61) by a method of successive 
approximations. In most cases the fastest procedure would be to use an ap- 
proximate value of y to estimate & and 7p. One uses this to estimate the values 
of the third, fourth and fifth terms of equation (58) and log, leaving v un- 
known everywhere else in equation (58). Equation (61) then becomes a qua- 
dratic equation for the unknown value of v/*. A solution of this quadratic 
equation can then be used to re-estimate £ and 7). Successive approximations 
will converge very fast because the estimated y is used only in small or slowly 
varying terms. 

Equation (58) has been compared with the experimental data described by 
KESTIN and WANG [4] and by KEsTIN and MoszynskI [3] for gases of accurately 
known viscosity. In the range where equation (58) is supposed to apply, the 
theory and experiment agree to within the estimated errors both of which are 
about 0-5% in this case. More detailed comparisons will, however, be described 
elsewhere. 
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Zusammenfassung 


Es wird untersucht, inwieweit eine diinne aufgehangte Scheibe endlicher Ab- 
messung, wenn sie mit kleiner Amplitude schwingt, sich zur Messung der Zahig- 
keit verwenden lasst. Die Bewegung der Fliissigkeit wird unter der Voraussetzung 
analysiert, dass die Grenzschichtdicke klein ist im Vergleich zum Radius der 
Scheibe, aber gross verglichen mit der Scheibendicke. Die Formeln zeigen die 
Abhangigkeit der Frequenz und des logarithmischen Dekrements der Schwingung 
von der Dichte und der Zahigkeit der Fliissigkeit. 
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Fehlerabschatzungen bei Randwertaufgaben partieller 
Differentialgleichungen mit unendlichem Grundgebiet 


Von LotHar CoLiatz, Hamburg?) 


Randwertaufgaben bei partiellen Differentialgleichungen mit unendlichem 
Grundgebiet lassen sich oft durch einfache Transformationen in solche mit 
endlichem Grundgebiet iiberfithren. Da Aufgaben mit unendlichem Grundgebiet 
haufig vorkommen, ist es vielleicht berechtigt, Typen solcher Aufgaben zu- 
sammenzustellen, bei welchen sich Fehlerabschatzungen fiir Naherungs- 
l6sungen nach einfachen Randmaximum-Prinzipien durchfiihren lassen. Unter 
anderem zeigt sich dabei zum Beispiel, dass dies unter gewissen Voraussetzun- 
gen bei der (bei Umstrémungsaufgaben auftretenden) 2. Randwertaufgabe in 
mehr als 2 Dimensionen méglich ist. Die folgenden Abschatzungsmethoden 
funktionieren bei den Umstrémungsaufgaben der Potentialtheorie gerade in 
3 Dimensionen, versagen in 2 Dimensionen und bilden somit in gewissem 
Sinne ein Gegenstiick zu den Methoden der konformen Abbildung, die bei den 
genannten Aufgaben in 2, aber nicht in 3 Raumdimensionen anwendbar sind. 
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I. Elliptische Differentialgleichungen 
1. Die Differentialgleichung 
Vorgelegt sei die homogene Differentialgleichung 


i . 0U 
LU =— P(X) A,U + 27 On(%}) ax, + S(X,) U=0 (1) 
fiir eine Funktion U(X,,..., X,) der  reellen Veradnderlichen X; in einem 


offenen Bereich B des x, -Punktraumes, der fiir ein passendes R, alle Punkte 
mit R > R, enthalt; dabei sei 


|x| = be x] ” (2) 


7=1 
gesetzt. 
Ax bedeute den Laplaceschen Operator 


P,Q,,S seien gegebene Ortsfunktionen, die in B stetig sind. P(X;) sei in B 
positiv und besitze fiir R > co einen Grenzwert P > 0. Die Menge J’ der Rand- 
punkte von B bilde einen stiickweise glatten Rand. Durch reziproke Radien 
werden neue Koordinaten eingefiihrt: 


m= st, X=4, = D/H), rR=1, (3) 
,= 
und als neue abhangige Veranderliche wird 


‘ u(x,) = 72-” U(x, 7—*) (4) 
benutzt. Mit 


Pe) = P(%,7-7),  Oxlx,)=7> SOT), s(x,) = 7-4 S(x,7-*) (5) 


lautet dann die neue Differentialgleichung (mit 4, als Laplaceschem Operator 
mn den neuen x,-Koordinaten) 


Ag om 7 (X;) 
ee 6 
c=s—(n—2)y, y= 2, mh | 


oe tole ;) 4 =0 mit ie. 2m09,| 


3ei der Transformation (3) gehe B bzw. der Rand J’ in einen Bereich B, bzw. 
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ei ° . . . . ens 
I, iiber, und es sei B; bei Hinzunahme des Nullpunktes ein einfachzusamm 


hangender, beschrankter Bereich. 
Es wird nun vorausgesetzt, dass die Funktionen #, 9;, 5 und damit auch 
g, und c im Nullpunkt 7 = 0 stetig sind (bzw. fir r = 0 so definiert werden 


: 


1 


| 


kénnen, dass Stetigkeit bei y = 0 besteht) und dass iiberdies c >0 in B; und : 


p >0in B, +I; gilt. 


Betrachtet man nun nur Lésungen U der Differentialgleichung (1), die fiir 


R-> oo von der Gréssenordnung U = 0(R?-") sind, so hat w fiir 7 >0O nach 
(4) die Gréssenordnung wu = o(1/r). Dann hat wu bei y =0 eine hebbare Singu- 
laritat, das heisst, « kann fiir 7 = 0 so definiert werden, dass wv in ganz B; + I 
eine mit stetigen partiellen Ableitungen bis zur 2. Ordnung einschliesslich ver- 
sehene Lésung von (6) ist?). 


2. Die Randbedingung 


Zu der Differentialgleichung (1) trete eine Randbedingung auf J” hinzu. 


MU =A,(X,) U — A,(X,;) => =As3(X;) (auf I”). (7) 
Die A,(X,) sind gegebene Funktionen auf J”. 

Es werden folgende 3 Randbedingungen betrachtet: 
1. Randbedingung: A, =1, A,=0. 
2. handbedingung: 3A, —O07A,—1. 
3. Randbedingung: A, =1, A,+0: in diesem Falle werde A,=0 vorausge- 
setzt. 

In jedem Punkte von I" sei entweder die 2. Randbedingung oder eine der 
beiden Randbedingungen 1,3. vorgegeben. 

Es werde bei der 2. und 3. Randbedingung iiberdies angenommen, dass in je- 
dem Punkte von J’ entweder von vornherein (Fall 1) eine eindeutig bestimmte, 
ins Innere von B weisende Normalenrichtung definiert ist oder (Fall 2) eine 
Normalenrichtung so definiert werden kann, dass ein endliches Stiick die- 
ser Normalen N zu B +I’ gehért. Damit ist Bedingung (7) erklart. Bei der 
winkeltreuen Transformation (3) geht die Richtung von N in die Richtung der 
inneren Normale v auf J’, iiber, und es gilt 


0U Ou Or 
om et (my 2) yn 
ON a Ov G hae Overe (8) 


(Auch im Falle 2 wird durch (3) eine Richtung v festgelegt, und es gilt (8).] Die 


*) Dies gilt sogar unter etwas schwacheren Voraussetzungen. Vergleiche Pu. HARTMAN, On 


the Local Behavior of Solutions of Au = g(x, u, Vu), Comm. pure appl. Math. 9, 435-445, speziell 
437 (1956). 
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neue Randbedingung fiir w lautet dann 


g(%) My = g(x) w — a,(x,) 2 2% = 72 a,(x,) (9) 
mit 
a,(x;) =A,(x,7-?) (fir k =1, 2, 3) | 
und ay (10) 
| glx) = (4) = (m= 2) a(n) r 2. | 


Der Randoperateur M; ist durch (9) erklart, falls g(x,) + 0 ist. Es sei g(x,) > 0. 
Es gilt nun der erweiterte Randmaximumsatz: 
Eine Lésungsfunktion (x;) von (6) lasst sich in B, im Falle der 1. und 
3. Randbedingung unter den getroffenen (ence soutien und fiir a, 20 
‘durch den Maximalbetrag von M, ,y auf I’, abschatzen$) : 


|p(x;)| S Max|M; g| (giiltig in B;+J°)). (11) 


3. Diskussion der verschiedenen Randbedingungen 


Bei der Abbildung (3) geht die 1. Randbedingung wieder in die 1. Rand- 
bedineung uber, A,=4,=g=—1, A,=a,=0. 

Fiir die 2. Randbedingung werde 0R/ON =O (Bereich B als «Stern- 
bereich») vorausgesetzt, dann ist 07/0y <0 in den nach (3) entsprechenden 
Randpunkten von I’;. Es ist jetzt A, = a,=1, und (10) ergibt 


g(x) =—(n—2) 7 5. 


Im ebenen Falle (n = 2) wird g =0, und man erhalt wieder eine 2. Rand- 
bedingung auf /’;. Aber in mehr als 2 Dimensionen (n > 2) entsteht die 3. Rand- 
bedingung, und zwar mit dem <richtigen» Vorzeichen bei der Normalableitung, 
so dass die Abschatzungsformel anwendbar wird. . 

Die 3. Randbedingung geht bei der Transformation (3) wieder in die 
3. Randbedingung iiber, und zwar wegen 07/0v < 0 wieder mit dem richtigen 
Vorzeichen bei der Normalableitung. Rechnet man M;q@ nach den Trans- 
formationsformeln in Gréssen des Aussengebietes um mit 


so ergibt sich 


M,o= ae aay (4,0—4, a) (12) 


3) Vergleiche zum Beispiel L. Corrarz: Numerische Behandlung von Differentialgleichungen, 
2. Aufl. (Berlin, Géttingen und Heidelberg 1955), S.372; oder Schiffstechnik 4, 67—68 (1957). 
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Die Abschatzung (11) lautet daher 


OA 0®@ aes bs 
|O(X,)| S jane Max (7 |4,0- 415771) (giiltig in B+), — (13) 


und es gilt der 


A AS CATE 


Satz: Es sei B(X;) eine Lisung der homogenen Differentialgleichung (1) mit 
stetigen Koeffizienten P, Q;,, S von der Gréssenordnung ® = 0(R*—"), und es sev — 


P>0Oin B+I. Die durch die Transformation (3) mit reziproken Radien nach 
(5) hervorgehenden Funktionen p, q;,, s seien 1m Nullpunkt r = 0 stetig, und es 
gelte R8S => (nm — 2) D' XQ. 

k=1 


Ausser den obigen allgemeinen Voraussetzungen tiber den Bereich B und 
seinen Rand I" gelte bei der Randbedingung (7) entweder Ay=1, A,z=O oder 
A,=0, A,g=1. Im Falle der 2. und 3. Randbedingung set B ein Sternbereich 
mit OR|ON = 0. Dann fallt die Grosse 


A, -0OR 
R an Ze 


g=A + (= 2) 


aus, und es gilt die Abschatzung (13) bet der 1. und 3. Randbedingung und bet der 
2. Randbedingung, falls die Dimensionszahl n > 2 ist. 

Fehlerabschatzung: Ist V eine Naherungsfunktion fiir die (als existuerend vor- 
ausgesetate) Losung U der Randwertaufgabe mit der inhomogenen Differential- 
gleichung LU = H(X,) und der Randbedingung (7), wobet V die Differential- 
gleichung exakt erfiillt und die Gréssenordnung o(R°-") hat, so gilt fiir den 
Fehler ® =V —U die homogene Differentialgleichung LO =0. Ist D der Defekt 
in der Erfiillung der Randbedingung 


Ds AeA, ay, (15) 


so geht (13) tiber in 
= 1 es coined Paley 
\V-U| S par Max (" 5 ||) (giiltig in B+). ~ (16) 


Bet der 1, Randwertaufgabe ist g =1. 


4. Beispiele 


1. Erste Randwertaufgabe der Potentialtheorie fiir den Aussenraum eines 
Ellipsoids. 

Es sei U=1 auf dem Rande J” des abgeplatteten Rotationsellipsoids 
X?+ Y?4 2 Z%=1 und AU =0 ausserhalb des Ellipsoids. Im Unendlichen 
gehe U gegen Null. 


(4) 


se ee aN 
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Fiir eine Naherung V, die AV = 0 erfiillen muss, werde der Ansatz 


Vz a, Ve (17) 
k=1 


/ 


gemacht mit den Potentialfunktionen (0? = X2+ Y?2 und R? = 0? + Z?) 


1 e?— 22? , 3et— 24.92 224824 
: Kulp ee ge : = + 72 
(18) 
V.= a (5 0° — 90 0! Z2 + 120 02 Z4 — 16 Z9). 


Die Fehlerabschatzung (16) besagt hier 
|V —U| <5 Max|RV—R|  (giiltig auf und ausserhalb I’). (19) 


Die Konstanten a;, sind so zu bestimmen, dass auf /’das Betragsmaximum von 
RV —R moglichst klein ausfallt (man zeichnet sich dazu zweckmassig den 
Verlauf von RV,, fiir k =1,..., m langs eines Ellipsoidmeridians auf). Man 
erhalt : 


Benutzung von Naherungsausdruck V 


2y2 
0,901 V, +0,05858 V, 
0,9005 V, + 0,075 V, + 0,0038 Vy 
0,90013 V, + 0,075 V,+ 0,005 26 V; + 0,000053 V, 


V, und V, 
V,, Vz und V, 
V,, Vo; Vz und V, 


1 1 
V; 0,8535 V, (a= Sr ‘) 


Benutzung von Fehlerabschatzung A Eeepie die 
ee — 10 eal 
V, V—U|s 2 0,147 U — 0,854| S< 0,147 
V, und V, V—-U|s a 0,041 U — 0,784| < 0,041 
V,, Vz und Vz, V—U/s e 0,0158 U — 0,7809| < 0,0158 
V,, Vo, V; und V, V—-U\s = 0,006 63 U — 0,7837| < 0,00663 


2. Zweite Randwertaufgabe eines Umstrémungsproblems. 

Es sei U das Geschwindigkeitspotential fiir die Potentialstromung um das 
Ellipsoid 7’: X2+ Y2+Z?/2 =1, das in Richtung der Z-Achse mit der Ge- 
schwindigkeit 1 angestrémt werde. Auf J’ gilt dann OU/ON = 0. Zieht man das 
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Potential u — Z fiir die nicht durch den K6rper gestérte Potentialstro6mung aby 
so liegt fir U = Uz U genau eine 2. Randwertaufgabe der oben betrachteten — 
Art vor; es ist AU =0 ausserhalb J’, im Unendlichen verschwindet U, und auf 
That man 0U/ON = —OZ/ON. Fiir eine Naherung V werde der Ansatz (17) 


mit den Potentialfunktionen 


| 
| 
: 
| 
: 


Z Yd, 9 5 
W=%: y= 242-60); Ve= fr @Z4— 40240 +1504; 
ie = (16 Z$ — 168 Z4 9? + 210 Z? p* — 35 0) 


[Abkiirzungen o, R wie bei (18)| 


durchgefiihrt. Die Fehlerabschatzung (16) lautet jetzt‘): 


1 
|V —U| S$ Max! 


I 


R OV OZ \ | 
eee a A 


Man erhalt mit diesen wenigen Ansatzfunktionen noch etwas grobe Resultate; 
es soll hier auch nur die Methode illustriert werden. 


Benutzung von Naherungsausdruck Fehlerabschatzung 
V, 0,9V, 
V, und V, 0,7V,+0,1V, 
V,, V. und V, 0,65 V,+ 0,1 V, + 0,032 V, 
V,, V2, V, und V, 0,63 V, + 0,1 V,+ 0,037 V,+-0,011 V, 


Es sei bemerkt, dass nur der Einfachheit halber ein Rotationskérper mit axia- 
ler Anstrémung gewahlt wurde, die Methode ist natiirlich fiir beliebige (nicht 
notwendig rotationssymmetrische) Kérper, sofern sie Sternbereiche mit stiick- 
weise glatter Randflache sind, und fiir beliebige «Anstellwinkel» verwendbar. 


4) Meinem Assistenten, Herrn Dipl. Math. Alfred Meyer, danke ich herzlich fiir die Durchfiih- 
rung numerischer Rechnungen. Er hat ausserdem eine entsprechende Fehlerabschatzung fiir die 
kompressible Unterschallstrémung um einen rotationssymmetrischen K6rper bei axialer Anstré- 
mung aufgestellt, was demnachst veréffentlicht werden soll. 
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II. Parabolische Differentialgleichungen 


5. Monotoner Charakter der Aufgaben 
| Die Differentialgleichung laute 


TU =U,— F(X,, ag om U,, U;.) = 6 (21) 
mit den Abkiirzungen 
U Ow. mA — 00 
er OK? tee CO Xeea a 


Es tritt zu den Veranderlichen Xj, ... , X,, gegeniiber Nr. 1 noch die Variable ¢ 
(etwa die Zeit) hinzu. 

Im (X,, ..., X,,)-Raum sei B ein offener, zusammenhangender Bereich mit 
stiickweise glattem Rand J”. Auf B seien die Anfangswerte von U vorgegeben. Es 
sei B=B+T, und Berl ys 3 seien die Punktmengen des (X,,..., X,,, t)-Rau- 
mes, fiir die X,,..., X,, zu B bzw. I’ bzw. B gehoren und 0 <t <¢#, mit 
festem ¢, gilt. 

Die Funktion F sei monoton nicht fallend in den U;; bei fester Wahl der 
Argumente X;, ¢, U, U;, das heisst, aus a; = 6; (fir 7 =1,..., 1) soll 


F(X,,t, U, U;, a) 2 F(X;,, t, U, U;, ) (22) 
folgen. 
Auf I, sei fiir 0 << ¢< t, wieder die Randbedingung (7) gegeben, und es 
werden wieder die 3 Randbedingungen wie in Nr. 2 betrachtet mit denselben 
Voraussetzungen iiber A,, A,und Ag. (Es lassen sich durch die folgenden Betrach- 
tungen ohne weiteres auch Randbedingungen erfassen, bei denen der Rand I, 
mit ¢ variiert ; jedoch sei der einfacheren Schreibweise wegen hiervon abgesehen. 
Auf B seien die Werte von U gegeben; das ist in (7) fir 4; =1, Ag=0 enthalten.) 

Wenn B beschrankt ist, hat der Operator monotonen Charakter®): Fiir 
2 Funktionen U, V (es wird hier stets vorausgesetzt, dass die auftretenden 
Funktionen U, V stetige partielle Ableitungen so hoher Ordnung besitzen, wie 


sie in den Uberlegungen vorkommen) mit 
TV<TUinB, und MV <MU auf Bund J, (23) 
silt V <U in ganz B,. Hat die Differentialgleichung (21) die speziellere Form 
TU =U, — g(U) —H(X,, t, U;, U;;) =9 (24) 
g(U) stetig differenzierbar, 4 in den U;; monoton nicht fallend], so folgt bei 


5) L.Coxiarz, Fehlerabschatzungen fiir Ndaherungslosungen parabolischer Differentialgleichungen, 
Anais Acad. Brasileira de Ciencias 28, 1-9 (1956). 
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beschranktem B aus 
|\TU| <e,, |TV] Se,inB, und |MU-MV|SéaufB+T: | 


Fir g’/(U) <0 gilt 


|U —V| $64 (e+e) ¢ in B,, (25) 
und fiir 
g(D) Sy “anit ys 0 
gilt 
~ =< prt yt &2\ le 26) 
[Os (One eee (26) 


Ist tiberdies V = 0 eine Lésung von (24), so gilt im Falle g’/(U) <0 sogar der 
gewohnliche Randmaximumsatz: Eine Funktion U mit TU <0 (bzw. 20) 
kommt auf dem Rande B +I’, der oberen (bzw. unteren) Grenze aller Werte, 
die U in B, annimmt, beliebig nahe. 

Erstreckt sich der Bereich B ins Unendliche, so hat man ihn durch passende 
Transformation so in einen endlichen Bereich iiberzufiihren, dass die eben 
genannten Sadtze anwendbar werden. Hierzu ist nicht immer die Inversion (3) 
brauchbar. Es liegt nahe, allgemeiner anzusetzen 


u(x;, t) = a(x;) U(d, (x), t) mit Noe Op ee (27) 


Es sei a(x;) > 0. Bei 6; werden Ableitungen nach dem Argument durch Striche, 
bei a und uw Ableitungen nach x; bzw. ¢ durch Indizes 7 und ¢ bezeichnet. 
Man erhalt 
pre pets) nyOP gis 
Cty wale Mann? eels GE 
fs Ad Li 
jig bie 


(28) 


(a2, — (WY) + 2a.4,8,) uy + (aa, Bf +2420; — aay, dul, 


Dabei geht (21) in eine neue Gleichung iiber: 


U, 


Tu=u,— I(x; t, U, U; YU; 


laa 


wobei f dieselbe Monotonieeigenschaft beziiglich der u,; besitzt wie F beziiglich 
der U;;. 

Man hat also nur nachzupriifen, ob die Randbedingungen nach der Trans- 
formation die geforderten Voraussetzungen erfiillen. Dabei geht das Verhalten 
von U im Unendlichen ein. Ist zum Beispiel gefordert, dass U fiir R ->+ co die 
Gréssenordnung U = 0(R*~") hat, so hat bei der Inversion (3), (4) w die Gréssen- 
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ordnung o(1), also den Grenzwert 0 fiir y +0; man hat dann die 1. Rand- 
bedingung w = 0 fiir 7 =0 und 0 <t<4,. Die weitere Diskussion verlauft ent- 
sprechend Abschnitt I. 


f 


6. Speztalfall von nur einer Raumkoordinate 


Es werde der Spezialfall n = 1 herausgegriffen. 
1. Es liege ein «halbunendlicher» Bereich vor: B sei das Intervall 


Ap F100, 


wobei X,>0 angenommen werden darf (X steht fiir die Veranderliche X,). 
la) Es werde die folgende Transformation verwendet: 


= oo u(x, t) =xU(=, i), eet iM Uren U,.. (29) 
Verlangt man von U Beschranktheit fiir X >cound 0 << #,, so hat man 
fir uw die 1. Randbedingung u(0, 7) =0 auf dem Randteil x =0,0<t<4%,. 
Die Randbedingung 
oU 
A,(t) U — Axl) (53°) = Aad) 
geht tiber in 


[4.4 — 4a] «+ 40) 2 5% = 45). 


Die 1. Randbedingung (A4,=1, A,=0) wird also wieder eine 1. Randbedin- 
gung, und die 3. Randbedingung (A4,=1, A,=0) eine 3. Randbedingung, 
wobei auch das Vorzeichen bei 0u/0x = —0u/dv die geforderten Voraussetzun- 
gen erfiillt, sofern man nur 1/x — A, > 0 sichert, was durch Wahl von X, bei 
beschranktem A, stets méglich ist [es muss X» > A,(X) sein]. Dagegen hat 
0u/0x bei der 2. Randbedingung nicht das gewiinschte Vorzeichen. 

1b) Man benutzt dann die folgende Transformation. Bei 


1 1 
de an u(x,t) = (a — x) u(>. i 
und der 2. Randbedingung (A,=0, A,=1) fiir X =X, mit X,>0 erhalt 


man als transformierte Randbedingung 


also die 3. Randbedingung, wobei das Vorzeichen bei 0u/0x die obigen Voraus- 
setzungen erfiillt, wenn nur a — x, 20 gilt, also@21 /X» gewahlt ist. 


: 


128 : LotHaR COLLATZ ZAM E 


2. Es seien lings der ganzen X-Achse (—oo < X < oo) die Werte von U | 
vorgegeben, und fiir X¥ > +00 sei U(X, #) fiir jedes feste ¢ beschrankt. Bei 


\ 


der Transformation 
X =tgx, u(x,t) =cosx U(tgx, 2), 


Uy =u sinx + 4, COsx, 
= 3 
Uyy = (u + 4, ,) cos*x 


hat man u(+<2/2,¢) =O und fir w(x, t) eine 1. Randwertaufgabe. Allerdings 
tritt hier in der transformierten Gleichung neben w,, das Glied uw auf, so dass © 
zum Beispiel bei der Warmeleitungsgleichung Ux, = kU, nach (26) die Fehler- 
schranke in ¢ exponentiell anwachst. (Man kann erreichen, dass kein solches 
Glied mit w auftritt; der Faktor von w in (28) verschwindet im Falle 1 =1 
fiir a(x) =1/[c, + cy. b,(x)] mit beliebigen Konstanten (c,, c,); die Transforma- 
tion (29) ordnet sich hier ein.) Doch ist dem Verfasser keine brauchbare Trans- — 
formation bekannt, die in den Fallen 1b und 2 das exponentielle Anwachsen der 
Fehlerschranke vermeidet. 


Summary 


Types of boundary value problems of partial differential equations for infinite 
domains are discussed which can easily be transformed in such a manner as to 
allow estimations of error (for approximate solutions) similar to the boundary 
maximum principle. First, second und third boundary value problems for the 
outer domain of linear elliptic and certain linear and nonlinear parabolic differ- 
ential equations are examined. For elliptic differential equations one of the 
results is that the second boundary value problem for more than two dimensions 
can be included. The estimates of the paper can thus be applied to problems of 
flow around some object, not in the case of two but of three dimensions. This is 
in a certain sense a counterpart to the conformal mappings method which is 
successful for two but not for three dimensions. Numerical examples show that 
estimations of error can easily be carried out. 


(Eingegangen: 8. Oktober 1957.) 
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Introduction 


Une équation différentielle étant donnée, quelle équation aux différences _ 


faut-il lui faire correspondre pour obtenir des valeurs approchées de ses solu- 
tions? La recette classique pour le remplacement des dérivées successives d’une 


fonction g(x) par des différences 3a, y(x;), est simple et facile a appliquer ; mais 
1 


elle repose entiérement sur des propriétés des polynomes; cette recette ne tient 
aucun compte du fait que la fonction p satisfait précisément a une equation drffé- 
rentielle connue. C’est pourquoi |’équation aux différences classique ne fournit 
en général pas les valeurs exactes de w, méme lorsque l’équation différentielle 
est simple et se laisse résoudre élémentairement. 

Dans la premiére partiz de cet article, je me contente d’étudier quelques cas 
simples, ot la solution exacte est connue. Cette étude préliminaire conduit a 
plusieurs méthodes générales. — Celles-ci sont ensuite appliquées, dans la 
deuxtéme partie du travail, 4 la construction de bonnes équations aux différen- 
ces dans des cas non triviaux. — Dans ce travail, j’insiste davantage sur les 
méthodes générales que sur telle ou telle formule approchée obtenue grace a 
elles: d’ot la place considérable accordée a la premiére partie, qui ne traite que 
de problémes «triviaux». — Je ne m’occupe ici que d’équations linéaires. 

Les principaux résultats ont été annoncés dans deux notes [6, 7]?) aux 
Comptes rendus de l’Académie des Sciences de Paris. — J’exprime ma vive 
gratitude au Centre National de la Recherche Scientifique a Paris et au Fonds 
National Suisse pour la recherche scientifique (par l’Ecole Polytechnique Fé- 
dérale): c’est grace a leurs bourses (1954-55, resp. 1955-57) que le présent 
travail a pu étre réalisé. 


PREMIERE PARTIE: 
Equations aux différences exactes 
I. Equations différentielles ordinaires 4 coefficients constants 
§1. Equations aux différences exactes. — Fonctions de neeuds et fonctions de cellules 


1.1 La solution générale de l’équation différentielle est connue: il est dés lors 
clair que les méthodes qu’on va développer ne sont intéressantes que dans la 
mesure ou elles se laissent transposer a des cas non triviaux. 

Soit L{p]=0 1 équation différentielle donnée dansun intervalle 1:0<x< nh: 
on impose en outre des conditions aux extrémités: nous y reviendrons. Inté- 
ressons-nous aux valeurs de w aux points 0, h, 2h, ..., nh (A est la «maille») et 


2) Les chiffres entre crochets renvoient a la Bibliographie, page 179. 


ies 
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cherchons 4 quelle équation aux différences elles satisfont exactement. Pour ; 
cela, nous considérons d’abord la méthode la plus immédiate: : 
L’énumération des solutions fondamentales de L{[y] = 0: on construit une 
équation aux différences satisfaite par chacune d’elles. 
Cette facon de procéder est tautologique: quand on connait la solution géné- 
rale de l’équation différentielle, nul n’est besoin d’équations aux différences! 


fa 


Nous écrirons l’équation aux différences 


Oy P(X) + OP (Xo + A) + nD (%o — A) + Hon Y (Xo + 2 A) | (1) 
+ ony (%—2h)+---=0, | 
ou un nombre fini de coefficients « sont différents de zéro, sous la forme 
Lig] = <Li, y> = 0°) 
avec 
Li = a Ox) + %, Ok + hy 9 x4—h) T Xp Cae on + Gop O(z,—2h) ee 


6(x,) est la mesure de Drrac au point Xp. 

On voit que les équations aux différences classiques sont exactes pour les 
solutions des équations différentielles (2), (4) et (6), qui sont des polynomes. 
Dans les autres cas, les équations classiques sont données par les quelques 
premiers termes du développement de l’équation exacte pour h +> 0. 

Remarque. Ainsi s’explique en partie la réputation des équations aux diffé- 
rences d’étre inutilisables pour le calcul des valeurs propres supérieures: pour 
un h donné, les coefficients des équations classiques sont d’autant plus mauvais 
que 7 est plus grand. 


1.3 Condition d’encastrement 


1.31 Comment faut-il traduire dans le langage des équations aux différences 
la condition d’encastrement 


p(0) = 9'(0) = 0 (8) 


pour une fonction g(x) solution, dans un intervalle 0 < x < a, d’une équation 
différentielle donnée L[p] = 0? — La recette classique est simple: on prend un 
point virtuel « = —h et l’on y pose 


p(—h) = H(A) ; (9) 


le raisonnement heuristique suivant y conduit: la premiére dérivée de y qui 


3) <f, 2> désigne le produit scalaire file) a(x) dx. 
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puisse ne pas s’annuler en 0 est la seconde: au voisinage de 0, on approxime 
p(x) par c x, qui est une fonction paire. — Critique: la relation (9) n’est exacte que 
si p est elle-méme paire; c’est le cas si L[g] =p”. Mais en général ce n’est pas 
le cas: l’expression finie de la condition aux limites doit dépendre de L’équation 
différentielle (c’est du reste-évident). 


1.32 Cas ou L{g] = d*y|dx* = 0. Solutions fondamentales: 1, ah, (x]h)®, (x[h)3; 
seules les deux derniéres satisfont (8) ; la condition (9) est inexacte, car la fonction 
(x/h)? ne la satisfait pas! On cherche une relation m(—h) = « y(h) + B @(2 hy), 
satisfaite par (x/h)? et (x/h)*; donc le vecteur (—1, «, 8) est paralléle au produit 
mctoriel (1, 1,4) x(—1, 1, 8) = (4, —12,2);a=3, B= —1/2. 


oh) =3 oh) — > 92h) (10) 


est la condition exacte d’encastrement correspondant a l’équation d'p|dx4 = 0. 
Cette relation équivaut, vu l’équation finie g(—h) — 4 (0) + 6 g(h) — 4 @(2h) 
+ o(3 h) = 0, a l’équation exacte 


9 (h) — 5 g(2h) + G3 h) = 0; (10' 


tandis que (9) équivaut a l’équation inexacte 7 m(h) — 4.9(2h) + 9(3 A) = 0, 
qui n’est vraie qu’asymptotiquement pour / > 0: sacrifice inutile. 

Remarque. Cet exemple simple montre aussi qu'il est artificiel d’exiger que 
léquation aux différences donne un systéme symétrique en méme temps que 
Véquation différentielle: car ce n’est pas le cas avec la condition exacte (10’). 


1.33 Cas on Lig] = d*p/dx* — 44 = 0. Solutions fondamentales satisfaisant (8) : 
Ch(7y x) — cos(y x) et Sh(n x) —sin(y x); écrivons de nouveau C = Ch(yh), 
etc., nous obtenons 


2S C=—2 se S—s 
g(—h) = 1 + er a ae p(2h) . (11) 


La relation (10) est le début du développement de (11) pour 7 > 0 (tandis 
que (9) est le début pour / > 0, et n’est exacte pour aucun 7). 


1.4 Fonctions de cellules (ou d’intervalles). La valeur p(%») de la fonction y en 
un point x, n’est qu’un cas particulier d’une fonctionnelle locale Wil = <P, 
produit scalaire de m avec une distribution Y a support compact: si ¥ = 6(,,), 
Yip] = p(x). Nous considérerons le cas oti le support de W est dans un inter- 
valle T: OX x%<h; a chaque intervalle [,: x, 5% < %p4, = %, +h de lon- 
sueur i correspond une distribution Y;,: la translatée de W ayant son support 
dans I,,; nous définissons une fonction de cellules (ou d’intervalles) « = u, par 


Pégalité 
: eth eis, PD. 
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Cas particuliers: (a) Si Y% = di) est la mesure de Dirac, u,(I,) est égale a 
la fonction y Al’origine du segment J: c’est la fonction de neeuds des paragraphes © 
Tet ee | 

(b) Si ¥ = —0’, les valeurs de w sont celles des dérivées de 9p. 4 

(c) Si Y est la fonction caractéristique de I (Y= 1 sur J, Oailleurs), alors” 


Uy Ty) =/¢ dx . 
Ik . 


Considérons maintenant une équation différentielle L{g] = 0 a coefficients | 
constants; cherchons une équation aux différences LU[w] = 0 satisfaite par la 
fonction d’intervalles u(I,) = <Y%,, p> pour toute solution m de L{y] = 0! On 
peut naturellement énumérer les solutions fondamentales de L{y] = 0 et cal- 
culer les fonctions d’intervalles u correspondantes. Il suffit cependant de re- 
marquer que <Y;,, p> = (W(x), p(x + x,)) pour déduire de (1) la relation 


Lu] = 9 U(Mp) + 4, U(L) + wy, U(L_4) + oer u(Ip) +--+ =0, (12) 


ot la somme est en réalité finie. Nous voyons que, quelle que soit la distribution 
WL a les mémes coefficients que L*. Bien entendu, cela n’est vrai que pour la 
«traduction fine» de l’équation différentielle et non pas pour les conditions aux 
limites! 


§ 2. Méthode: construction d’une distribution 


2.1 Les considérations suivantes sont indépendantes de l’hypothése que les 
coefficients de L{q] = 0 soient constants. 

Nous revenons au probleme initial (§ 1.1) de déterminer une équation aux 
différences L*|p] = 0 satisfaite par toute solution qm de Lig] = 0. L’équation 
aux différences cherchée n’est autre qu’une combinaison linéaire des équations 
différentielles relatives a tous les points d’un segment: nous multiplions done 
scalairement l’équation L{p] = 0 par une distribution T a support compact: 

= <T, L{y]> = <L[T], p>, ou L est l’opérateur différentiel adjoint 4 L. Nous 
voulons donc avoir, conformément a ieiat 

L{T] ra Ly = Xo 9 (x) + Oy Ox +h) + Oy Oem + My Ox 42h) a Ses (1”) 

les coefficients a, a,,... (en nombre fini) sont a déterminer en sorte qu'il 


existe une telle distribution JT. Mais, comme nous allons le voir, la distribution 
T elle-méme n’est pas sans intérét. 


2.2 Si l’équation différentielle donnée est L[g] = f(x), nous obtenons immeé- 
dvatement sa «traduction discréte»: c’est l’équation aux différences 


Lily] = <T, p>. (13) 
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Pour autant qu’on ne s’intéresse qu’aux valeurs de m en un nombre fini de 
points, nul n’est donc besoin de connaitre exactement la fonction f(x): seuls 
jouent un rdle les produits scalaires <7, />, qui sont en nombre fini. 


2.3 Cette méthode permet également de traduire par des équations finies 
les conditions linéaires aux limites imposées a la fonction g(x). Supposons que 
Ja fonction @ soit cherchée dans un intervalle 0 < «<b =~mnh, et que l’on im- 
pose au point 0 les conditions aux limites suivantes, en nombre fini: 


Lite] = Bio (0) + Bi. (0) + --- =8;, 
Lig] = Bi PO) + Bi, 9 (0) +--+ = B;, 


PSs ellie) @ se else 6 =v « 6 + 2 ‘piv 6 wis © 6 ‘8 jeje 9 0. (0% 0 


ou toutes les sommes ont un nombre fini de termes. Nous construirons alors des 
distributions T telles que ces expressions données appararssent: 


L*(q) = <L[T],9> =) 90) +, P(A) +772, 9(24) +--+ GQ Lily] t+ o Lolglt+-:- 
= Yo PO) + ¥, P(A) + °° pasa 


c’est-a-dire, si nous écrivons L? = B;)(—1)° 6 + B;,;(—1)! 6’ + --- , etc. (6 = dq), 
i= LT) = 0 J.) eee s din) T Yar On) Sa aaa ae Ts 1B oe Ri, 


les coefficients y et ¢ doivent étre déterminés en sorte qu'il existe une telle 
distribution T. Si l’on a L{p] = f(x), l’équation (13) reste valable pour notre L 
défini selon (1”). 

Les conditions aux limites se traduisent par autant d’équations aux diffé- 
rences complémentaires. Le cas le plus banal est la condition aux limites 
y(0) = a (a donné). 

Tout cela est valable, que les coefficients de L[p] soient constants ou 


variables. 


B4 E xemples 


2.41 L’exemple le plus trivial est fourni par l’«équation» g(x) = f(x): Popé- 
rateur «différentiely L = 1 est de degré zéro; L = 1; on veut avoir 1: T= 6;,,, 
donc T= 6,4); LE =L[T) = T= 6y,, Vou naturellement 9(%) = f(%0)! 
(L’équation aux différences ne concerne que la valeur en un point.) 


2.42 Exemple un peu moins trivial: Lip] = oa) = fe): LIT)=—T" 
= Oy Orn) + %—a Oem F devant étre A support compact, JT = 0 en dehors du 
segment [x — #, x]; T est constante a l’intérieur de ce segment: nous y posons 
med? dota =1, «,=—1; Li = —(x,—n + Ox); Péquation aux diffé- 
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rences (13) devient 
= 9 (9 — h) + ole) = | fx) de 
Xy—h 


2.43 Lip] = y'(x) + 9 v(x) = f(x); alors 


LT) ==—T' +9 T = 0904) + ty Oy -9s 
nous prenons 
0 en dehors de [%» — h, Xo]; 


Tig) = 
ert" Y dans (45 ty Xols 
dou a, = el et a-, =—1; donc 
~ 9 (=f) + €™* lag) = | f(x) en dx 


ty —h 
(en accord avec la méthode de variation des constantes). 
2.44 Pour l’opérateur L = d?/dx*, on construit la fonction 
0 pour *#%—h; 
T(x) = K—KX+h pou %74-~ha% S&S 4%; 
=i Xp h pour % Se Sx kh; 
O pour  *% 2% +h; 
dot L¥ = L[T] = bu, — 2 62) + S(x,+% - 
2.45 Pour l’opérateur L = d?/dx* + 7, on construit la fonction 


0 pour *=*,—h; 


- sin[4(% — x» + h)] pour %,—hS%x 5%; 

TX) = 
Les 

; sin[y(—*% + %)+h)) pour x S*Ss%+h; 


0 pour * 2x +h; 
dott 18g = nas a => Ce h — 2.cos (y h) O(x4) + O (45 +h) . 
2.46 Pour l’opérateur L = d*/dx4, on construit la fonction 


4 pour |2| 22h; 


(2h — |€|)* 


T(x) = pour A= |é| =2h; 


6 
2h = AES — |&|)9 ; 
Fepeceetien) pour |é| Sh, ob =x —%; 


dot Li = L[T)= O(z,—2n) — 4 O,—» + 6 Oy, — 4 Oe, 4m + Ol, 42) * 


soon ye ee ta son “peony ecianecas emmmanmnanammaiiaeas 
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Montrons ici, a titre d’exemple, la construction détaillée de T(x): T étant 
a support compact, T(x) =0 pour |x —x,| > 2h; comme L[T] = a*T |dx4 
doit étre une combinaison linéaire de mesures de Dirac 6 et non pas de 6’, 6” ou 
6”, T(x) doit étre continue ainsi que ses dérivées premiére et seconde; T" sera 
discontinue aux seuls points %)— 2h, %)—h, %, X%)» +h, %)+ 2h. Désignons 
les intervalles 


Li Kye =. yeah — hs Ie ay — We oS Ky: 

le Hee KG Ee Li t+th Sx Sr++ 2h; 
et posons %, = x — (%) — 2h); x, = %— (%)— A); Xz =H —Xq; Hy = X— (Xp +h). 
Mans chacun des intervalles J, 1,, 1,, 4, L[T] = dT /dx¢ = 0, T y est donc 
un polynome du troisiéme degré. — Pour x, << 0, nous avons JT = 0; donc en 
m0: 7 = 7 =—1" =; d’ow dans J,: T = const - x3; nous normons cette 
constante en sorte d’avoir 77” (x, = +0) =1: 


1 


sponte: 


aes 
En x, = h—O: 


1 cas Meares: 
Fe Wie TD IA Dheh BE" 1; 


nous devons donc avoir dans J: 
1 1 gucate, 1 1 
Tah + > Wag t Shag t k= = (xy + H+ (z - a) x3, 
ou la constante k reste a déterminer. En x, on a %, = hi: 
Ge Pe (ee. 2 n Win 
Pes (G+ A) 2, he ([+3A) A Perit 6 hh ee Tan be 


Dans J, et I,, nous posons naturellement T(x) = T (2 %) — x); la continuité de 
T’ en x, impose donc que T’(x)) = 0, c’est-a-dire k = —1/2. Nous avons alors 
Bons)t: 2” =1; dans 7,: T” =6k=—3; dans J,: T” = +3; dans I: 
T” = —1; nous obtenons bien 


Ly = Tage) ie 6 (4,—2h) = O (x,—b) + 6 D(x.) ert O(5+h) r O (x5 42h) c 
2.47 Pour l’opérateur L = d*/dx* — yj, on construit de méme la fonction 


0 pour |é| = 2h; 


VAS 
iS) 
> 


Pos Sh [24 — [2|)] — sinly (24— |B) 1 aang 
T(x) = Zoe {8+ Ch fy Ge [E/N] — (C +29) Shy JE) 


—s-cos[y (h— |é|)] + (2C + ¢) sin[y (h — |E|)3 
pour |é| <A, 
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ob é=x—2x,, S=Sh(nh), C=Ch(nh), s=sin(nh), c=cos(nh); 


Vow LF = LiT)= On on 2 (CPO 0e Se (2574 Ge) 0.,\ 
—2(C+¢) 0.45 + Ou,zaee 


2.48 Condition d’encastrement p(0) = y’(0) = 0 avec lopérateur L = d*/dx'. 


| 
| 


Nous étudions la fonction w dans un intervalle 0 < x < mh. Les distri- | 


butions T et les équations aux différences relatives a ]’intérieur de l’intervalle 
ont été construites au § 2.46 ci-dessus. y(0) = 0 est maintenant connue; mais 
il faut une seconde équation finie pour traduire la condition d’encastrement. 


Conformément au § 2.3, nous construisons une distribution T telle que L[{T] soit — 


une combinaison linéaire de 6), O(0), On), O20)» On): 

Désignons les intervalles 1:0 <% S h ie bh SS ZN; 1 2k eee 
Nous posons T = 0 en dehors de ces trois intervalles. 

Dans J,: de facon analogue au § 2.46, on a 


i 
T= = (3 h = Pa 
d’ou dans J: 
1 a 
T =~ (3h—x)>+ (k—Z) 2h— x), 


la constante k, étant a déterminer. 
Dans J,: T et T’ doivent étre continues en 0, donc nulles; 


L = heh HR 
I] nous reste a déterminer k,, ky, kg par les trois conditions que T, T’ et T” 


soient continues en «=h: on obtient kj = —7/12, ky = 3/2, kg = —11/12; 
donc dans J,: 


3 5 1 35 mM 11 " 
Le a x 12 ee Le ana ea "(Ge 0) S32? 
dans J. 
T==(3h-x8—+ 2h—08, TY=1; 
dans [ 
i 3 Ja ehae al eee 
dot 
ger oh es ad*T 11 9 
Lp = LT | = por = pr Oa ey Oi Gap eee 
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‘Si d*p|dx* = 0, léquation aux différences est donc 


0 = Lylgl => 90) — 3h y'0) + 99H) — 2 g2m) +9 BH) 


= 9 o(h) — 5 g(2h) + g3h) | 


en accord avec (10’). 

Si d*py/dx4 = f(x), il suffit de remplacer a gauche 0 par <T, f. 

Toutes nos formules sont naturellement en accord avec celles que nous avons 
obtenues par énumération aux §§ 1.2 et 1.3; la connaissance des distributions 
T permet en outre de transposer immédiatement nos équations aux différences 
pour l’équation différentielle L[g] = f(x). — D’autre part, ces distributions 
peuvent étre utilisées avec succés pour le calcul approché de problémes a 
coefficients variables (calcul local des perturbations, cf. chapitres III et IV). 


2.5 Comparaison avec d’autres méthodes aux différences 


2.51 Notre méthode par construction d’une distribution a plusieurs points 
communs avec celle de H. SASSENFELD, que R. ZURMUHL ({12], p. 404 ss.) 
appelle Quadraturverfahren: (a) les deux méthodes font usage de l’équation 
différentielle donnée en tous les points d’un segment; (b) l’intégration par 
parties intervenant dans la méthode de SASSENFELD correspond a notre 
passage de l’opérateur L a l’opérateur adjoint L. (c) Cependant, contrairement 
a SASSENFELD, nous construisons la distribution 7 qui seule peut fournir 
léquation aux différences exacte; c’est pourquoi nous ne sommes pas amenés 
a une équation intégrale. 


2.52 Contrairement aux deux méthodes comparées ci-dessus, celles que 
L. Cottatz et R. ZuRMUHL appellent Mehrstellenverfahren et Differenzenver- 
fahren hoherer Annaherung n’utilisent l’équation différentielle qu’en wn nombre 
fini de points. Le dernier procédé a le défaut important de faire intervenir dans 
chaque équation finie plus d’inconnues qu'il n’est nécessaire. — Nous verrons 
cependant (§§ 7 et 9) que le «Mehrstellenverfahren » donne des résultats souvent 


proches de ceux de notre méthode. 

2.53 Il est essentiel, pour la deuxiéme partie de ce travail, que notre 
méthode fournisse l’équation aux différences exacte lorsque les coefficients sont 
constants. Cet avantage n’est pas réalisé par les autres méthodes: c'est pourquoi 
je peux dire qu’elles consentent des sacrifices inutiles. 


2.6 Une loi de composition pour les opérateurs aux différences 
Soient 
LO*[g] = So), p (vo + mh) et L*[p] = D/ mh P (Xo + mh) 
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deux opérateurs aux différences a coefficients constants ; je définis leur produit 
de composition comme celui de deux distributions: ; 
co 
2 x 12 1 Zz A 
Tee LY*)[g] = Ds Ob p(x tkh), ou Oy) = Pais Ot ond 
m=—co 
les sommes n’ayant en réalité qu’un nombre fini de termes. 
Loi de composition: Soient L" et L®) deux opérateurs différentiels linéatres a 
coefficients constants; alors 
(L! Lee 25 Los * Les; 


de plus, si T;™ assure le passage de L® a L®* et T,2) de L°) a L)* (cf. § 2.1), 
alors la distribution 
Tae) = Tw * Tre 


obtenue par produit de composition4), assure le passage de LY L°) a (L® L®)¥. 
Démonstration: Par hypothése, 


OT) =LO* et L° [Te] = Le*; 
donc 


L® 


(2) [dpa * 2,e)= LT a] * L?T,@] = Lo* * LP i+ : 


Exemples: Reprenons les cas traités au §1.2 (et au § 2.4): il suffit de 
décomposer 


et 


pour voir que les opérateurs aux différences exacts (4), (5), (6), (7) peuvent 
étre obtenus par composition a partir de (2) et (3). Par exemple: 


= inh »& 
(4. + 1 ii)” = Oo ah e! yh oh 
et 
ia Sead \ te é 
sei fs, nee —inhk S 
(as in) = — Gq) +e" Ow ; 
donc 


a ; 
(aor + 1)" = 8:4 — 2.008 (% h) by + Spy - 


On an 24. également la loi de composition pour les distributions T construites 
au 


4) Si T,(x) et T,(x) sont des fonctions, (a Ca == Hh T,(€) Ty (x — &) d&. 


= Sico 
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_ 2.7 Généralisation: distributions To conduisant a des équations aux différences 
_ pour des fonctions d’intervalles 


: 2.71 Comme au §1.4, nous considérons une distribution Ya support dans l’inter- 
valle [0, 2), et nous nous intéressons aux valeurs de la fonction d’intervalles h 
définie par “,(I,) = <Wy,y>; ou J, est un intervalle x, <x«< nay = Xt h 
du réseau et Y;, est la translatée de Y, A support dans J,. 
De méme qu’au § 2.1, nous multiplions scalairement l’équation différentielle 
L{g] = 0 par une distribution T™ a support compact: 


0 Sr, Lig) =< Lil), os 


mais a présent nous voulons avoir 


EB eet aa op oy oe oo Beye 


et les coefficients «4 sont 4 déterminer en sorte qu’une telle distribution T? 
existe. Nous aurons alors 


: 0= LE[u,) Sol u(h) + a7 u(h)+ a9, u(L,)+- 


Si L[q] = f(x), il suffira de remplacer 4 gauche 0 par <T", f>. 

On remarque (comme au § 2.2) que les valeurs u(J;) sont déterminées dés 
qu’on connait ces produits scalaires <7", f>, qui sont en nombre fini. Il n’est 
pas nécessaire de connaitre exactement /(x). 


2.72 A lintérieur (traduction de l’équation différentielle) : 
Soit T la distribution assurant le passage de L a L;* (pour les fonctions de 
neeuds, cf. 2.1); L étant a coefficients constants, la distribution 


TUO= TV, 
obtenue par produit de composition, assure le passage de L a L;,; et les coefficients 
Z 
oo ol, ob, ... de LP sont tdentiques aux aX, %p, py ahs ++ fal Be 
En effet: 


ee LT) 4 Y= fo Oe oy 05,45) + os) ea aM ta, +--. 


L’équation aux différences est donc, a l’intérieur de Vintervalle considéré, 
indépendante de la distribution W- cela, en accord avec le § 1.4. 

2.73 Ilen va tout autrement pour les extrémités: on doit maintenant construire 
des distributions TS (indépendamment des 7), telles que (je garde les notations 
du § 2.3) 


LP =L(T T= ye ye ty et tt oe lets 


on doit prendre les coefficients y~ et £7 en sorte qu’une telle TU existe. 
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2.74 Exemples de conditions aux limites. Je me borne au cas ou VY est la fonction 
cavactéristique de Vintervalle I, =[0,h): P(x) =1 pour0 S% <h, 0 ailleurs. 
(a) Equation donnée: dy/dx = f(x) pour x > 0, condition en O: (0) = 4. 

d 


ae dnt sO 
ee je cherche T“ telle que a YE ee Oe 


c’est-a-dire: TS doit étre A support compact, continue sauf en 0, et de dérivée 
premiére constante dans (0, A), nulle ailleurs. Je prends donc 


Te =h—x«x. por O< xh, Oaillews; y=1, C=—h; 


d’ou 


(E5, fy =<To, pS = CT, 2 = 4 prXr-h Oy, >= 4h) — BF; 


(c’est d’ailleurs trivial). 

(b) Equation donnée: d*p/dx? = f(x) pour 0 <x < uh; condition en 0: 
y(0) = a. L = L = d?/dx?; je cherche TU(x) telle que TO” = y, ¥o+ y, ¥4+ 6 b@) 
(I, est Vintervalle h < x << 2h). TO doit étre continue, T™’ continue sauf en 0. 
Je prends: 


se i he. 
Leto = = 2 Wise pl ees ee 
Alors 
Pe! == 3 + ek (0) » 
dott 


pe KL Eg ye Oe op 5, ye ee 


Si f =0, on peut done écrire u_, — 2 up + uy, = 0 avec (cellule virtuelle) 
U_y~=2ha— uy (réflexion al extrémite). 

(c) Equation d?p/dx? = f(x) pour 0 < x < mh; condition en 0: ’(0) = 0. 
Je cherche T(x) telle que TO” = y, % + y, ¥, + ¢ d). TU doit étre continue 
sauf en 0, 7’ — € 6 continue partout. Je prends: 


r 1 - 
Ri TI=5 (e- 2h; Re Tl s- lea. 


4 


Ona 


Le ey ae 693 <LOR= To, pO SI = 4, Ee 


Si f= 0, on peut écrire w_,— 2% + u,=0 avec u_, =u) — h2b (réflexion). 
(d)_ Si Péquation donnée est d4p/dx4 = f(x) pour 0 < x < nh, avec les 
conditions en 0: p(0) = a et ’(0) = d, on traduit ces deux conditions par deux 
equations aux différences: l’une («centrée en Jj») fait intervenir Uy, Uy et Us; 


——_ 
2 


i 
7 


el \ re i 


“"— 


_ Vol. IX a, 1958 Contribution a la méthode des équations aux différences 143 


autre («centrée en J»), a9, w,, uy et us. En particulier, sia =) = 0 et f =0, 
on obtient 
85 23 5 23 


eee i hy ==) et our =z 


UW —4u,+u,=0. 


Il reste bien entendu que les coefficients ainsi déterminés ne sont valables 
que lorsque Y, est la fonction caractéristique de l’intervalle J,. 


2.75 En vertu du § 2.72, on a de nouveau ici (LM L®)9 = 1100 4 12/0 en 
définissant ce produit de composition comme au § 2.6; mais il n’est pas vrar 
que Tr) 2) = Pow * Tra): le§ 2.72 montre que To = T « ¥, donc 


O cs <= = Leh Fee ea 
Tr 1) = Tape) ie oa Ta) : T ) PS Ta) * Te) ae Pray # To): 


§ 3. Méthode récurrente 


Contrairement aux méthodes exposées ci-dessus, celle dont nous allons parler 
est indépendante d’une connaissance préalable des solutions de 1’équation 
différentielle. De plus, elle n’est nullement restreinte, par son principe, a 
l’étude d’équations linéaires (bien qu’ici nous nous bornions a ce cas). — Cette 
méthode jouit d’une grande autonomie vis-a-vis du calcul différentiel. Nous 
’appliquerons, par la suite, 4 des problémes aux dérivées partielles. 


3.1 L’idée de base est la comparaison des opérateurs aux différences Lj et Life 
(ou bien Lj? et Lf). —Essayons de déterminer une suite d’opérateurs discrets 


Li, Life, Lif, ..., Lifon, ... tels que les solutions Uj, Uyjo, ... des equations 


L¥[u,]} = 0, Lifo[Upjo] = 0, .-. coincident en tous les neeuds relatifs a la maille h 
(opérateurs cohérents)! Si, lorsque N > oo, léquation aux différences 


— Libon[ Myon] = 0 devient équivalente a une équation différentielle L{~| = 0, nous 


saurons que la solution y du probléme continu prend exactement les mémes valeurs 


que U;, en tous les neuds ou celle-ci est défime. 
En d’autres termes: L’éqguivalence de l'opérateur discret Li et de l’opérateur 


continu L s établit en trois éapes: 
(a) Passage récurrent de Ly a Liisa; 
(b) Résolution de la récurrence: trouver une expression des coefficients de L#§ 


telle que la récurrence soit satisfaite; 
(c) Passage ala limite pour h > 0: Lit > L. 


Remarque générale: La difficulté pratique de la méthode rvéside dans la deu- 
xieme éape (b). On peut considérer cette éape (résolution de la récurrence), 


pour laquelle une recette générale fait défaut, comme correspondant au probleme de 


LV intégration d'une fonction donnée. 
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3.2 Il est souvent plus facile de calculer la récurrence Lj > Li,,. Nous parta-_ 
geons l’intervalle J considéré en 2 n intervalles J,..., 2, de longueur 4; nous ~ 
obtenons 27+ 1 nceuds, dont 2 —1 sont intérieurs 4 J. Appelons ™ la 
fonction de nceuds solution de L¥#[m,] = 0; 7’appelle uz, la restriction de uj, aux 
neuds «pairs» (dont la distance aux extrémités est un multiple entier de 2 h). 
Il faut déterminer des équations aux différences L#,[u;,] = 0 satisfaites par u2;! 

Comme nous ne considérons ici que des équations linéaires, chaque opérateur | 
L#, sera construit comme combinaison linéaire de plusieurs opérateurs Lj. Nous 
avons ici l’ analogue discret du § 2.1. Au lieu de passer du continu au discret, nous 
passons maintenant de la maille / a la maille 2 h. — Les coefficients de la com- 
binaison linéaire cherchée peuvent étre eux-mémes considérés comme les valeurs 
d’une fonction de neuds T,, (correspondant a la distribution T du § 2), a support 
compact. Nous complétons notre formalisme discret par analogie: 

(Up, Uy, > = D3 v,(x;) u(%;); <L¥o,], > = Op L#lmy) - 
neuds x7 

Nous avons 


0= (7, Liu) = ATi u,> = LF, [M4] , 


a condition que L[7;,] s’annule sur tous les nceuds «impairs». On construit 7; 
en sorte qu’il en soit ainsi, et l’on obtient L*, en fonction de L*. — Si l’on avait 
Ly [u,] = v,, UL faudratt remplacer a gauche 0 par <T;,, v;,>. 


3.3 Exemples 
3.31 L’équation 


dp 
art A= 0: Je pose L¥[u,] = —u, (x —h) + A, u,(x) =0, 
ott A, est a déterminer. (a) La récurrence est immédiate: 

L¥[M.,] = —™, (« — 2h) + A,, u(x) =0 


avec A,, = Aj. (b) On satisfait la récurrence en posant A;, = e*”, la constante 
a étant a déterminer. (c) Passage a la limite h > 0: 


U(x) =p, u%(x—h)=p—ho'+ OW), A,=e™=1+ah+ O(R), 
donc 
O= Liu) =—p+hg'+(1+ah) p+ O(n), d’ot wy’ +ap=0; 


pour obtenir l’équation voulue, il faut donc poser a = n. L’opérateur L¥ est 
déterminé par A, = e", 

| A titre de comparatsons : (x) Equation aux différences classique: A, =1 + nh; 
c’est la «fonction linéaire tangente» a la fonction exacte A, = e”", 
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(B) «Mehrstellenverfahreny: 
h? 


(x — B) = glx) — ho'(s) + > p"(a) + O04) | +2 


y (x —h) = Sats 
d’ou par combinaison linéaire: 
2 (x —h) + ha! (x —h) = 2 q(x) — hg'(x) + O(F); 


on utilise l’équation différentielle w’ + 4 y = 0 aux points x et x —h: 


(2 — 4 h) p (x — h) = (2+ uh) yx) 


et l’on pose donc 


n? h? ie Ve ; 
hes Oa ; 
c’est la «fonction bilinéaire osculatrice» de 


2 778 fps 
mho1+yh+ 4 ee m gebee 


3.32 Corde vibrante: Cet exemple est particuliérement important, car il me 


servira de modeéle. 
Equation différentielle : 


we tp =0 (0) = p(t) = 0 
ee P= , avec 90) =gp(il)=0. 


Je pose 
L¥*[u,] = , (x — h) — A, u,(*) + 4, (4 + h) = 0 


(A, est a déterminer). 
(a) 1°) Méthode: combinaison linéaire directe: [v, = ™,(x), etc.] 


I 
. 
— 


U, 94 — A, Uy 1 dts 


u a Ay U, + Uy 4 =) as 


x—h 


U, — Ay Uy in + Ue s2, =O] -1- 


: 
ine age teen = QO avec Aj, =A, -~ 2. 


ZAMP IXa/10 
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2°) Construction d’une fonction de noeuds 7;, telle que Es (Th) @& LEG) 
soit différente de zéro aux seuls neeuds x — 2h, x, x + 2h: 


x—2h x x+2h 
| onarte 
Eth) = 1 As, I 
T, = Lf aed 
dou 
A,,= A, —2. 


(b) Résolution de la récurrence: A»,/2 = 2 (A,/2)? — 1; comme cos(2 «) 
= 2cos?« — 1, la récurrence est satisfaite si l’on pose A, = 2 cos(k h). 

Remarque: Si le probleme physique n’était pas continu mais discret, la 
grandeur physique étant non pas 7 mais 4,, on aurait pour / = ne: 


A; 
Aj. = 20s (n arceos“;" ; 


soit A,,, = 2 T,(A,/2), T,,(x) étant le n-iéme polynome de TCHEBYCHEFF. 
(c) Le passage a la limite # > 0 donne ¢” + hk? mg = 0: on rejoint |’équation 
différentielle donnée en posant k = 7; Ly est déterminé par 


A, = 2 cos(7 h) . 


Calcul des valeurs propres, en choisissant par exemple 4 = 1/4: Compte tenu 
des conditions aux limites wu) = u, = 0, nous avons trois équations linéaires 
homogénes pour 14/4, Uy/2, Ug/q: 


—A Uy)4 + Uy)9 = 0 | —A 1 0) | 
Uy —A My t Ug, =0, dct O=| 1 —-A 1|=A (2— A?); 
Uy, — A tg, = 0 od ae 


donc cos(7/4)=0 ou + 1/V2, Q= Mm avec m= 172,35, 6, 1 Ue eee 
c’est-a-dire m = 0 mod 4. 

L’absence des valeurs propres n= 4N x ne doit pas nous étonner: car les 
fonctions propres Wyn correspondantes s’annulent aux cing neuds x = 0, 1/4, 1/2, 
3/4, 1 et induisent donc la fonction de neuds triviale U,-1)4 = 0. — Ce phénoméne 
(«éclipse» de certaines valeurs propres) est fréquent. — Comme la fonction 


propre fondamentale n’a pas de noeuds a l’intérieur de l’intervalle, la fréquence 
fondamentale ne s’ éclipse jamais. 


A titre de comparaison: 


n4 h* n® hé igre 
2, 360 “ 


2 cos(y h) = 2 — 72h? + 
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a) Equation aux différences classique: 


A, =2— 77 hi; 
6) «Mehrstellenverfahren»: 
tee ey he eee ae es ie ie? 
Ay 12+ 72h? ata hinge a tiger 


3.33 Barre vibrante: 04p/dx* — n* y = 0. Pour l’intérieur: Je pose 
Li [ul] = 4, ,+ B,u,_,+A,u, + Bu, .,+ 4 


K+2h° 


(a) Construction de 7, telle que 


x—4h x—2h M3 K+2h x+4h 

rane ] ] ] 
LT) = Le(T,) =1 By As, Boy 1 
T, = core hf pees aaa 


d’ou les formules de récurrence: 
Ay, =A 2B+2 tt Byy'2 A, = Be. 

(b) Résolution de la récurrence: la réside toute la difficulté. La récurrence 

est satisfaite si l’on pose 
Ai =4Ch(kh) cos(kh) +2 et- B,=—2Ch(kh) — 2cos(kh). 

(c) Le passage a la limite (4 >0) conduit a l’équation différentielle 
d4p/dx4 — k* py = 0; on rejoint donc |’équation donnée en posant k= 4: Lj est 
déterminé par 


A,=4Ch(yh) cos(nh)+2 et B,=—2Ch(yh) — 2cos(nh). 


3.34 Les conditions aux limites peuvent étre traitées de facon analogue 
(mais les calculs sont en général longs et la résolution de la récurrence (b) trés 


difficile). 
Exemple: Corde vibrante (cf. § 3.32 ci-dessus) a extrémités tenues élasti- 


quement: 
gy (0) — x p(0) =0. 


Je pose 
Xi, Uy + %, = 0 (X, a déterminer) . 
(a) Eliminons uw, avec tg — Ay Mj, + Yon = 0: 


(1+ A, X,) m+ %,=0, dor X,,=1+4+4,X,=1+4+ 2X, cos(yh). 
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(b) On peut résoudre cette récurrence en posant : 
X, = —cos(nh) + asin(n h) . 


(c) Le passage a la limite h > 0 donne y'(0) + 4 7 (0) = 0: il faut donc poser 
a= —x/n et nous avons 


X, = —cos(y h) — 7 sin (7h) . 


(Selon la recette classique, on aurait posé X, = —1 — h x.) 
Interprétation par la construction d’un nceud virtuel en x = —h: 


uU_,—-A,U+%y=90 avec 4,= [cos (x h) — a sin (7 n)| u, 
Extrémité fixe: 4 = oo, 4) = 0; extrémité libre: x= 0, #1, = %. 


3.4 Récurrence avec des fonctions de cellules 


3.41 Pour lintérieur de Vintervalle considéré, les relations de récurrence sont 
identiques a celles qu’on obtient pour des fonctions de neuds. — En effet: Appelons 
u,, une fonction de cellules pour la maille . Soient I et I) , deux h-intervalles 
adjacents, et 1?” = J UI |; nous définissons u,, par 


Mads hee =a u,(L Wy + B u,( (Ee oe 


Nous supposons que ™, satisfasse une équation aux différences LO[w,| = 0. 
Combinons linéairernent ces équations exactement comme nous l’avons fait 
(pour Lj‘) au § 3.2: nous obtenons des relations du type: 


j 
Oty My (LE) + hay, MALE op) + Eon Mp og) + Oy, (IM 4,) + = 0, 
~ (h ~ h ~ I > h 
Oly My (Tet a) + Say My(L0) 5) + Eo, M(LO,) + Oy, M(Te 54) +--+ = 0, 
OU %, X24, % 2h, X4y, +. Aépendent des coefficients de LO exactement comme 


les coefficients de L*, dépendent de ceux de L*. 

Multiplions la premiére de ces équations par «, la seconde par f et addition- 
nons: nous obtenons une équation L5}[w.;,] = 0 avec les mémes coefficients. 

N. B. Nous savions déja (§ 1.4, g 2. 72) que les coefficients de LO sont les 
mémes que ceux de Lt. — I] n’en est plus de méme aux extrémités! 


3.42 Exemple: corde vibrante (cf. § 3.32). — Je prends de nouveau (cf. §2.74) 
pour ¥, la fonction car actéristique de l’intervalle d,,5 done 
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Pour Vintérieur: Je pose 


Ly [mu] = u,(1@_,) - A, uy, (1) + (EO), \iee 0 


x +h 


Je dois combiner linéairement plusieurs équations de ce type; interprétons les 
coefficients de cette combinaison linéaire comme les valeurs d’une fonction de 
cellules T,: 


i a LF (m)> = LATA U,> = Lihue, ] 


a condition que l’on ait: 


x—2h Os a ie x+2h ee x+4h 
| | 
pie (ee Wale ee 
eee i A SAG 1; 
dot 
i= ul A,+1 A,4+1 i 


et par conséquent A,, = Aj — 2, en accord avec le § 3.32. 
Extrémité 0 tenue élastiquement: g'(0) — x w(0) =0. 
Résultats du calcul: Construisons une cellule virtuelle 7_,: —h <x < Oet 
posons l’équation u(I_,) — A, u(Jo) + u(Z,) = 0; il faut poser 


_ 4 sin(y h) — x[1 — cos(y h)] 
ae ‘1 sin(y h) + x [1 — cos(y h)] ud) - 
Extrémité fixe: ~(0) = 0, % = co, u(I_,) = —u(d). 
Extrémité libre: y’(0) = 0, x = 0, u(I_y) = + u(h). oy ak 
Dans ces deux cas particuliers, on a un simple principe de réflexion. 


3.43 Remarque générale: Le phénoméne de l’éclipse (§ 3.32) se produit égale- 
ment pour les fonctions de cellules: lorsqu’une fonction propre g(x) induit une 
fonction de cellules identiquement nulle u,: < Y;, p > = 0. 


II. Premiers essais de transposition 4 des équations aux 
dérivées partielles 


§ 4. Le probléme de Cauchy pour un demi-plan 
Dans le demi-plan supérieur y > 0, construisons un réseau quadratique dont 
les droites soient inclinées a 45°, les nceuds inférieurs reposant sur l’axe des x. 
Appelons h V2 le cété des carrés. QM, sera le carré de centre (x, y), c’est-a-dire 
de sommets (x —h, y), (« +h, y) (x, ¥—A), (%Y +4); T™) sera le triangle 
rectangle isocéle de sommets (x — h, 0), (« + A, 9), (x, h). 
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4.1 Pour lVintérveur: 


Pour les fonctions w, définies aux nceuds du réseau, soit Lj lopérateur aux 
différences défini, pour un carré OW pat 


* ee bx. = 
LO yl%] ale he, y + U, +h,y Ux y—h Us y+h 4 


ot j’écris u,, pour u(x, y). - En d’autres termes: l’opérateur L, fait corres- 
pondre a chaque fonction de nceuds une fonction de carré. 

(a) Récurrence: Partons d’un nceud arbitraire O sur l’axe des x, et définis- 
sons, dans le réseau précédent, les noeuds «pairs» comme ceux qu’on atteint en 


parcourant, a partir de O, des arétes doubles du réseau (de longueur 2 V2 h). = 


Nous considérons maintenant le réseau quadratique de cétés = 2 V2 h, dont les 
neeuds sont «pairs». Appelons u,, la restriction de uw, a ces noeuds; nous re- 
marquons que 


* ae ; me 
L¥ yx, y Men] a Uy 2h, y as Woh y ae U,, y—2h Uy yh 


est la somme des expressions L#[w,| relatives aux quatre carrés dont ge”) est 
la réunion: 
* ze pee fe * * 
Le ee lo ronety v eg Ey es + Le v-h x Che yth-: 


Par conséquent, si w, satisfait l’équation aux différences Li[a,] = 0,(Q), ot 
On est une fonction de cellules, il s’ensuit que ws, satisfait 1’équation 


LF Mon] = Oo,» ot 02,(Q°”) = s) 0,(Q) : 
ot) “E Ql2h) 


(b) La récurrence est satisfaite si l’on pose 
en(Q™) = k |] ox, y) dx dy, 
Qt’) 


la fonction k o(x, y) étant arbitraire. 

(c) Passage a la limite h > 0: Admettons qu’il existe une fonction p(x, Vv) 
deux fois continiment dérivable et telle que uw, , = = p(x, y) aux noeuds du ré- 
seau; développons: 


Lil] = Wi, — Pyy) + OM); 0,(Q™) = 2 oh? + O(A’) , 
c’est-a-dire 


0 = Lil) — 0,(Q) =? @,, — yy — 2% @) + O(}3); 
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d’ot, si nous posons k = 1/2, 


Pxx— Pyy = O(% ¥) et 0,(Q == | el, y) dx dy. 
g(t) 
4.2 Au votsinage de l’axe des x: Définissons encore Vopérateur L* pour les 
triangles T™ par 
if 1 : 
Ly. [uy] 2 U,—h-0 ~ ca U, 11,0 - Un 1 : 

(a) Récurrence: Soit (x, 0) un nceud pair; on vérifie immédiatement que 
| Be aa UT eon oF dE eeae ai LS cela signifie que, si Li; 21M] fa on( TO), 
il s’ensuit 

Qan(Te”) = On(T ey) + On(T Ps) + 04(QM,) - 


(b) La récurrence est satisfaite si l’on pose 
x+h 


o,(7) = a 5 if om) y) poy 
it) 
(c) Passage a la limite h > 0: 
Li {m,] = —hg, + OF); 9,(T%) =2oh+ OF), 


c’est-a-dire 
0 = L¥[u,] — 0,(7) = —h (gy, + 20) + O(F), 
d’ou 
1 
G(X) = an Py (x, 0) « 
En posant / = y, nous retrouvons donc la formule bien connue: 


p(%, ¥) = = [p(x — y, 0) + y (x + y, 0)] 


x+y Vv 


se XG 0) dé — 5 / » (Pex — Pyy) U aN . 


n=0 |gé—x| < y=—n 


§5. Traduction finite de I’ équation (23 = ») (ss + 9 a) p(x, ¥) = Qo 
Bek Reraderene un opérateur aux différences du type 


L¥{u] = A, Uy y +, Bi (Uz—n,» a Wy shy a Ux, y—h “s na y+n) 


ar Uy jh y—h oa Uy py +h at Un sn y—h Ag Uy, +h, v+h? 


ow les coefficients A, et B, sont constants. On peut caractériser l’opérateur 
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i par le schéma: 


ea Ba eee 
Ey tie ee 
i Gamal ee el 


(a) Cherchons a établir une récurrence h > 2h: nous voulons combiner 
linéairement quelques équations du type ci-dessus en sorte d’obtenir une équa- 
tion du méme type pour la maille 2h. 

Résultat: la récurrence n’est possible (exactement) que si A, = Bj; elle 
est alors donnée par By, = 2 — Bj et Ag, = B3,. On voit — et c’est un point 
essentiel — que la condition de récurrence A = B? se laisse elle-méme itérer. 

(b) J’écris la formule de récurrence pour B sous la forme 


elle est satisfaite si l’on pose 
B,=—2cos(nh), A, =4-cos* (7h) . 
(c) Le passage a la limite / > 0 conduit a l’opérateur 
O2 o? 
ht (So - 1?) lees ai ?) + O(hS) . 


Remarque: Si au contraire, partant de l’opérateur différentiel, on veut 
construire un opérateur aux différences correspondant, on retrouve immédiate- 
ment celui dont nous sommes partis, avec les valeurs ci-dessus de A, et B,: 
il suffit d’utiliser le § 3.32 et l’analogue pour deux dimensions de la loi de 
composition du § 2.6. 


5.2 Si l’équation donnée est 


0? 02 
(ier + i?) le oy r?) PY = 0p = const, 


nous posons l’équation aux différences sous la forme L{[w] = h4 0, (0, est une 
fonction de nceuds constante): 

(a) La récurrence est donnée par ,;, = @, cost (7 h/2); 

(b) elle est satisfaite si l’on pose 0, = (k/h*) sin*(y h/2); 

(c) le passage a la limite 4 > 0 monte qu’on rejoint l’équation donnée a 
condition de poser k = 24 99/74, c’est-a-dire 


= ( Ege TL ) 
Cn = Co h Din . 


En particulier, si 7 = 0, on a @, = 0) pour toute maille h. 
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DEUXIEME PARTIE: 
Calcul approché 
III. Equations différentielles ordinaires: Problémes aux limites 


§ 6. Méthode: calcul local des perturbations 


6.1 La méthode classique des perturbations ne s’applique avec succés que 
lorsque l’équation différentielle donnée est, dans tout l’intervalle considéré, «voi- 
sine» d’une équation dont on connait la solution exacte. — Au contraire, nous 
ne faisons pas cette hypothése globale; si les coefficients de l’équation différen- 
tielle sont suffisamment réguliers, ils varieront peu dans chaque petit intervalle 
(formé d’une ou plusieurs mailles): dans un tel intervalle, notre équation dif- 
férentielle «non perturbée» sera a coefficients constants. Nous construisons 
ainsi d’abord un systéme d’équations aux différences «non perturbées»: sa 
solution sera déja plus précise que celle de l’équation différentielle non per- 
turbée du calcul global. 

Calcul plus précis: Le terme perturbateur contient la fonction inconnue; il 
sera pris dans le membre de droite de l’équation différentielle; la perturbation 
de lV’ équation aux différences sera donnée par son produit scalaire avec une distri- 
bution T (cf. § 2.2), et calculée avec une précision a choisir. 


6.2 Le procédé que je propose ici se situe donc entre le calcul global des per- 
turbations et la méthode aux différences classique; cette derniére conduit a un 
partage fin de l’intervalle malgré une équation aux différences médiocre. — 
Si lon veut un résultat précis, il y a lieu de faire porter son effort non seulement 
sur le nombre des équations (diminution de la maille h), mais aussi sur leur 
qualité. 

Nous avons en méme temps |’avantage de disposer d’équations aux diffé- 
rences exactes pour le cas de coefficients constants (grace aux méthodes du 
chapitre I), c’est-a-dire en l’absence de «perturbation». 

Remarque: Les idées appliquées ici aux équations aux différences sont ana- 
logues a celles de Cu. BLANc ([2], p. 187 notamment), qui, pour le cas de coetf- 
ficients lentement variables, construit une fonction de GREEN solution (sauf 
en un point) d’une équation proche de |’équation adjointe. 


§ 7. Exemple: y" + (1+ x)y=—1 avec y(—1)=y(1) =90 


On cherche la valeur yo) de y en x = 0. 
On comparera avec CoLLatz [5], pp. 136 (équations aux différences clas- 
siques), 154 (équations aux différences portant sur 5 points), 158 («Mehrstellen- 


Perce! 
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verfahren»), 176 (calcul global des perturbations), 195 (méthode de RITZ), : 


211: au moyen d’une série de puissances, COLLATZ montre que 


v9 — 0,9320536| < 4,2 - 10-7; 


ainsi que ZURMUHL [12], pp. 398 et 403. 


7.1 Préparatifs: Nous partageons l’intervalle (—1, 1) en un nombre pair 2/h de 


mailles. Soit x) un nceud intérieur; pour toute distribution T a support com- 


pact, nous avons 


CT" + (ba 8) Toy) = <T,y" + 1 eer = = 7, = ev Phe 


Je pose 
1x =o) (n> 0) ,et. b= xX 


et je choisis naturellement pour T(x) la fonction indiquée au § 2.45: 


| 0 pour |é| 2h; 
T€é) Bean Acie ee 
| = sin{n (a — |é|)]} pour |é| SA; 
alors 
Pp Die Oy ay 2icos (4 h)iOgy dey s 
d’ot 


Vao—h — 2 cos (7) h) Vx, + Vey th eas ie (2 Xo § air e*) uy + 13 : 


I] s’agit de calculer approximativement le produit scalaire a droite: 


3 
' 1 
Y x9 zi “Gy Vx Wr Og"); V x5 = Qh (Yash = Vx.—n) —- O(h?) ; 


Va = = my, eke. 
Z 2 
(L,1>= na [1 cos(y h)]; <1 €> = <F,.€*> = 0 (symetrie): 


; ye n> ns : 
CTS?) = eae Say ae) len ee ue + Off) . 
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Si, dans le produit scalaire, on néglige les termes en h!, on obtient Véquation 
aux différences non perturbée 


Vek COSA) Y. + ¥, Se = <1, 1 = — [1 — cos( h)] . 


7.2 Résultats numériques (comparés a ceux des autres méthodes) : 


p Bot ione Equations | 
Méthode Perturbations locales, q ayes aon 
en négligeant les termes en aux Gitterenecclmatallon: 
différences 
ee pour verfahren» 
a | OFdinalres : 
Maille A ns 18 5 points 
B= 1 0,851 1 0,93003 | it - 0,8571 
1 
yo = a 0,922 0,9358 |0,932105| 0,967 0,9288 | 0,9276 
i 
= ar 0,927 0,9328 |0,932059| 0,949 0,9314 | 0,9312 
ut 
h= = 0,938 0,931 945 
| 
h : 0,9336 0,93205 
oe = , , 
10 


Remarque: Cet exemple a été choisi a cause des ouvrages cités; il est clair 
que, si (1 + x?) était remplacée par une fonction a comportement moins ré- 
gulier, les trois méthodes «classiques» (a droite du tableau) seraient d’autant 
moins adéquates a traduire le probléme fidélement dans la langue des équa- 
tions aux différences: car elles ne font intervenir que les valeurs de la fonction 
aux noeuds du réseau. — Ici (comportement régulier), le «Mehrstellenverfahren» 
s’avere excellent et fournit sensiblement la méme précision que la méthode 
des perturbations locales ou l’on néglige les termes en /°. 

A titre de comparaison: CoLLatz ([5], p. 196) a obtenu par la méthode de 
Ritz: avec deux termes y, ~ 0,9334; avec 3 termes yo ~ 0,932070, ce qui est 
excellent. 

Conclusion: Si l’on veut augmenter la précision, il est bon de perfectionner 
l’équation aux différences avant de réduire la maille. - Dans le tableau ci- 
dessus: les méthodes classiques aux différences progressent verticalement; au 
contraire, le calcul global des perturbations progresse horizontalement (partie 
de gauche de la 1 ligne); il est préférable de choisir une progression de type 


diagonal. 
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IV. Equations différentielles ordinaires: Problémes aux 
valeurs propres 


§ 8. Méthode 


Nous allons construire des équations aux différences presque exactement 
comme au § 6; la difficulté supplémentaire réside dans le fait que nous ignorons 
les valeurs propres exactes: il s’ensuit que Jes fonctions auxtliaires T(x) dé- 
pendent elles-mémes du parameétre inconnu A. — On est ainsi conduit a des égua- 
tions transcendantes. 

Cette circonstance est évidemment désagréable; mais il faut remarquer 
qu'elle est dans la nature méme du probléme. D’autre part, elle présente un 
léger avantage: on obtient du méme coup, comme racines d’une méme équa- 
tion transcendante, des valeurs approchées pour plusieurs valeurs propres. 


§ 9. Exemple: y"+A(1+ x2)y=0 avec y(—1) = y(1) =0 


On comparera avec CoLLatz [4], p. 297 (équations aux différences ordi- 
naires), et [5], pp. 235 et 240: 


2.71034 = A; SS 2A: 


ainsi que ZURMUHL [12], p. 403-405 («Mehrstellenverfahren ») et p. 461. 


9.1 Préparatifs: Je pose 7? = 4 (1 + x) (ot A est inconnue) ; I’équation donnée 
peut étre écrite sous la forme 


y+ wy = —A (x2 — x2) y. 


Ecrivons x — x) = & et multiplions scalairement cette équation par la distri- 
bution T déja utilisée au § 7.1: nous obtenons I’6quation aux différences 


Veoh — 2 COS (Hh) Vet Vee ART, (2%) E42 Se 
on a 1cl Vi =—y7? y,,; pour le reste (OLS, CD ESS ste) ere § 7.1. —Il est 
clair qu'il faut prendre d’autant plus de termes du développement que le 
produit 7 / est plus grand. 
Si l’on néglige les termes en h!, on obtient l’équation «non perturbée» 


Vig—h — 2 COS(H A) y, + Ven = Oh); 


dou naturellement une équation transcendante pour A. 
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9.2 Résultats numériques (comparés a ceux d’autres méthodes): Je me borne 
icl aux premiéres valeurs propres impaires (qui correspondent a des fonctions 
propres paires). 


Perturbations locales, Equations raced 
Valeurs Méthode en négligeant sae « Mehr- Wee be 
Ces les termes en différences bcd ,| H- BUckNER 
aha! —————_}—___—__—____+J ordinaires | Y°"*""™ [3] 
Maille ht h 
Ve es DI 2,47 2 2,4 2,286 
We = 2,221 2,180 2,025 2,209 2,189 
Zz i os 2,19 2,173, 2,078 2,186 2,179 
1 72 
1 
k= 2,129 
ame 
1 
jigs = 2,1774 
5 


20,1 i253) 17,21 16,84 


[sss (Sn 


ae |, 61,7 - & 
fe 
Ve a 55,5 
As 1 wt 
» 
1 
x 


Les valeurs approchées calculées par la méthode des perturbations locales avec 
la maille / = 1/2 sont racines d’équations transcendantes. 
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On voit que le «Mehrstellenverfahren» et le remarquable procédé de 
H. BicKner [3], qui ne donnent que relativement peu de travail, sont trés 
satisfaisants pour le calcul approché de la premiére valeur propre Ax, et cer- 
tainement de A,; par contre, notre méthode a l’avantage de fournir une vue 
d’ensemble sur un nombre considérable de valeurs propres. 

A titre de comparaison: CoLLatz ([5], pp. 235 et 240) a obtenu par la mé- 
thode de Ritz (avec deux termes) 


ApS 2177 486et Ay = 18300 


évaluations extrémement satisfaisantes. 
Je remarquerai de nouveau ici que le «Mehrstellenverfahren», pas davantage 


que l’équation aux différences ordinaire, ne peut tenir compte d’un comporte- 
ment anormal des coefficients de l’équation entre les nceuds du réseau. 


V. Equations aux dérivées partielles: Problémes aux limites 
$10. Quelques remarques générales 


10.1 Construction de l’équation aux différences: Il n’est (hélas!) en général pas 
possible de construire une équation aux différences «exacte»; on cherchera 
alors une équation aux différences fournissant une approximation la meilleure 
possible. - La méthode récurrente nous rendra des services appréciables au 
prochain chapitre, et nous comparerons ses résultats 4 ceux de la méthode 
classique et du «Mehrstellenverfahren». 


10.2 Si le domaine considéré n’est pas convexe, on utilisera avec avantage les 
fonctions de cellules (voir des exemples au chapitre suivant). Notamment, si 
l'on travaille avec des fonctions de nceuds, on introduit une grave erreur en 
prenant un coin rentrant du contour comme nceud du réseau. 


10.3 La traduction discréte des conditions aux limites pose un probléme trés 
difficile. - On peut transposer, par analogie, les régles (exactes pour une di- 
mension) obtenues au chapitre I. Nous le ferons au prochain paragraphe pour 
la condition d’encastrement: nous allons montrer que notre régle est meilleure 


que la recette classique (comme on l’a vu, cette derniére n’est méme pas exacte 
pour une dimension). 


S11. Exemple: Equilibre d'une plaque hexagonale encastrée, 
uniformément chargée 


Probleme aux limites: AAy(x, y) = 1 avec » = 0g/0n = 0 sur le contour 
d'un hexagone régulier de coté 1; on cherche la valeur de gp au centre. 


7 


“ 
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_ 11.1 Méthode classique: L’opérateur aux différences caractérisé par le schéma 


ci-dessous (CoLLaTz [5], p. 366) remplace Vopérateur différentiel 9 h4 44/16; 


et l’encastrement est traduit par la recette (9) du § 1.31. 


Extrapolation: on admet que (0) — um» est approximativement proportion- 


nel a h?; d’ou, si l’on a obtenu u™) et u\) resp. par les mailles h, et hg, la for- 


mule approchée 
ie un) — ht ye 


ee eggs eg 


11.2 Notre expression de l’encastrement: Nous nous inspirons du cas unidimen- 
sionnel, analogue au § 1.3: toutes les solutions F(x) de 


{ee = a = const be TO} LF (Or 0 | 


satisfont a 


Pour simplifier, nous négligeons, dans notre probléme a deux dimensions, la 
variation de 0*@(x, y)/On4* le long du contour de l’hexagone: nous appelons a 
cette quatriéme dérivée normale. 

Nous utilisons de nouveau l’opérateur aux différences classique (voir le 
schéma ci-dessus). . 

Extrapolation: nous ne connaissons pas la valeur de a; il est facile de l’ex- 
trapoler a partir de deux mailles successives /, et h,, en sorte d’obtenir une 
méme valeur approchée wu) = u\") (on a deux équations linéaires a deux in- 
connues @ et Up). 

Voici quelques exemples illustrant la construction de nceuds virtuels en 


s’inspirant de (10"): 
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Méthode classique: Notre expression de l’encastrement: 
h=1 h=1 h=12 
Susi b b 
3u+2 me Ne 
pan W ATIVAN 3 
DADA 
pe AALS 
DAA” 
0 uty JL b= 3v-! 
bb _ twa 
4 4” 64 


11.3 Résultats numériques: 


Notre expression 


Maill ee ag at ae Méthode Extrapolation 
aille 


en posant enectras classique classique 
a=0 polant a 
ee 0,01875 | 0,031 25 
1 0,01022 0,00946 
h= — 0,01119 0,01491 
A 0,01010 0,01021 
h= a 0,01040 0,012 30 | 


VI. Equations aux dérivées partielles: Problémes aux valeurs propres 
§12. Réseau triangulaire (neuds) ou hexagonal (cellules): membrane vibrante 
Nous voulons traduire l’équation 
Ag(x, y) + 2 ox, y) = 0 
par une équation aux différences 
L¥(u,] = 0 


pour les valeurs u, de m aux neuds d’un réseau de triangles équilatéraux de 
cété A («maille»). 


12.1 Nous allons procéder par récurrence (de 4 4 2h), en nous inspirant du es 
DOLENtIaN bye bane Cy, Cg, ... les valeurs de wu, (aux nceuds ainsi désignés du 
réseau) ; nous caractérisons l’opérateur Ly par le schéma suivant (a dec ite): 
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c’est-a-dire que nous écrivons 
6 
* — 
L* (a) = A, a = ) aes 
ae 


Nous avons donc: 
O= L¥[a]) =A, a => b (14) 


O= Sis = —6a+ (A, — 2) Dadi ca2 Na. (14') 


(a) Il faudrait combiner linéairement ces deux équations en sorte d’obtenir 
une relation du type 


et 


O= Li(a) = A,,a—)" ¢; (15) 


c’est malheureusement impossible. — I] faut donc sacrifier quelque chose lors 
du passage de (14) et (14’) a (15), et on n’obtiendra ainsi que des approxima- 
tions pour les valeurs et fonctions propres. 

On évalue facilement l’expression 


HW Ag|, +13 *) 


ew 


Zi [a] = (0, + by + Bs) — (ba + b4 + bq) =(3 4 + 


WAg|, +h +) = h Z[g), , 


Bl) Ww 


_ (3 a+ - 
ou Z[], dépend uniquement des dérivées de ¢~ en a. 
Calculons maintenant 
Z¥([Z¥(m)]. = Ze [br + 03 + 45] — Zi [ba + ba + Oc] 
=) c+ 2)7b-(6a+2)/ 4); 
en itérant le résultat précédent, on obtient donc 


64425543 0-2 a=W Z Zip], = Ol). (16) 


ZAMP IXa/11 
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Nous pouvons maintenant éliminer > b et Dd entre (14), (14’) et (16), 
ce qui donne 


A? , 
(A -24,) 4-7 = 00); (15’) 
cette évaluation remplace (15). Nous devons donc choisir 
2 
A. AS (17) 


N.B. On peut aussi déterminer A,, et les deux coefficients d’une combinai- 
son linéaire de (14) et (14’) en sorte que celle-ci différe de (15) par une quantité 
petite d’ordre le plus élevé possible en 4: par un développement en série (en 
utilisant par exemple le tableau IX de CoLLatz [5], p. 516), on obtient égale- 
ment (17). 

(b) Il s’agit 4 présent de résoudre la récurrence (17): on peut l’écrire sous 
la forme 

2A,,=A,—44,= (Apre2)* = 45 
soit 


Aun? 22 (ta?) ‘a 


La relation cos(2 «) = 2cos?« —1 permet de poser (A; — 2)/4 = cos(k h), 
soit 
A, = 2 +4 cos(k h) ; 


la constante k (indépendante de h/) est a déterminer. 
(c) Passage a la limite h > 0: 


on rejoint donc |’équation 


+ O(A4), 


| 
| 
a 


Ay +Agp=0 
a condition de poser 
V3a 
k — Wit 


Comme la récurrence n’est pas exacte, on n’obtient (pour une maille finie) 
U L > ’ ra a : 
qu'une valeur approchée 4 pour A. L’opérateur L¥ est finalement caractérisé par 


A, = 2+ 4cos(5-34). (18) 


Is : O L . soe 
L’équation aux différences fournit un nombre fini de valeurs possibles pour 


A n> pour chacune d’elles, on obtient en général une infinité de A, dont les pre- 
miéres sont effectivement utilisables (voir l’exemple, § 12.7).5) 


5 : P fp ae dee 
~ °) Si, pour un A, déterminé, les J étaient complexes, on pourrait prendre leurs valeurs absolues 
|4 | comme valeurs approchées de quelques A. 
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La récurrence est exacte si y est telle que Z# Zi(u,] = 0 partout et pour 


tout /; c'est le cas notamment si Z*[u,] = 0. Exemples: 


12.11 Fréquence fondamentale d’une bande de largeur a: 


ante Au—2u=0, A=2, Ay = exact . 
0 0 0 
u 
0 0 0 


B 12.12 Fréquence fondamentale d’un triangle équilatéral de cété a: 


h= =; Awa O57 “w= 0 donc 0=4=244008($]/4), 


= 16 xz 
(==. = CLA 
Ay 3 a2 


(Voir aussi la sixiéme partie du tableau numérique au § 12.7.) 


12.2 Autre méthode: on connait, pour toute fonction @ satisfaisant a 


Agy+Ag= 0, la «formule de la moyenne» 
Fe 1 ra A 5 : 
HbYA) ple) = ae f vot he dd =x tiv); 
0 


le terme de droite est égal a 


5 ; I . 
ED) tet he Py 401k"), 


m=0 


(19) 
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par conséquent, sil’on caractérise Lf par A, = 6 Jo(h Va), on a L#[u] = O(1°); 
au lieu de (18), on peut donc poser 


Ve uAClinik 3 (20) 


et I’équation aux différences L¥[w,] = 0 déterminera des valeurs approchées A; 
pour les A;. — Pour A, + 0, on ne peut obtenir qu’un nombre fini de valeurs 
approchées réelles [contrairement a (18)]. 


12.3 Méthode classique (COLLATZ [5], pp. 366 et 368): On pose 
A,=6-5 1A, (21) 


et l’on prend les valeurs propres discrétes A; comme valeurs approchées des 4;. 
En effet, si l’on pose A, = 6 — 3 h?A/2, on a L#[u,] = O(n’). 


12.4 «Mehrstellenver fahren» (COLLATZ [4], p. 297): On pose 
06 —18h2 A mee "Gate 


Ay = , Sot Ass . : 
16+ 2A NN gt 


m 
d’ot des valeurs propres approchées A;. En effet, avec 


96 —18-h2A 
162 


AN 
on a L#[u,] = O(h'). 


12.5 Comparaison des développements en série. — On voit facilement que, des 
quatre méthodes indiquées, la méthode classique est la moins précise: comme 


indiqué, Lj[w,) = O(h*) pour celle-ci, = O(h*) pour les trois autres. — Posons 
AVA =x; 
2+4 cos( 3 x) Mi = ate aa 5 ya > xe + O(x8) 
Z i See 1280 ; 
6 Iylx) 6S Sip sate 2) Ga Os 
0 LEV: 384 
3 3 
fee ae ay ee 
6 a * = 6) x x, 
96 — 18 w? 3 3 3 
16+ #2 bd a Ser mei yan 


(On constate que les deux premiéres expressions sont trés voisines.) 
Remarque: Si le bord de la membrane est élastiquement lié, on pourra 
l’exprimer en s’inspirant du § 3.34. 


a ll 
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12.6 Réseau de cellules hexagonales: On peut reprendre purement et simple- 
ment les expressions (18), (20), (21) ou (22) de A,; il faut seulement remarquer 
que, si h© est le cété des hexagones, h = V3 h©. Pour exprimer les conditions 
au contour, on construit des cellules virtuelles extérieures au domaine, en 
sinspirant du § 3.42 (principes de réflexion pour la membrane a bord fixe et 
pour la membrane libre). 


12.7 Exemple: Membrane hexagonale de cété 1, a bord fixe (p = 0) 

En méme temps que quelques valeurs propres approchées calculées par les 
procédés aux différences, j’indique des évaluations obtenues par des méthodes 
de G, P6LYA: 

(a) G. Pértya [8]: des régles générales énoncées dans cet article, il s’ensuit 
que les six premiéres valeurs propres A,, ..., Ag de notre membrane hexagonale 


‘satisfont les inégalités 


1 cercle dé rayon 
ds < sa AY le de vay sal 
a=1,...,6; y= rayon conforme (intérieur) maximum = 0,89850 pour notre 
hexagone. Comme on va le voir, ces évaluations par excés sont excellentes,; pour 
~ > 6, le membre de droite donne une valeur approchée de 4;. 

(’) Considération heuristique: avec 7 = 2,4048, on a (cf. POLYA et SZEG6 
(iL \p.. 268) 3-3 Vay y =1 pour le cercle, = 0,9964 pour le carré, = 0,9860 pour 
le triangle équilatéral; on a lieu de penser que 7+ Vay y > 0,9964 pour l’hexa- 
gone régulier. 

(8) G. Porya [9]: L’hexagone régulier étant un «pavé», 


k = surface 3 y3 


(y) G. Porya [10]: Pour tout «quasi-trapéze» (par exemple le trapéze qui 
forme la moitié de l’hexagone) de hauteur h et de «longueur moyenne» m, on a 


gee = L(="), 


ou ¢= L(g) est la fonction inverse de q = Vie 124. Viz 2 ee | eee 

avec r entier, 72 <¢ < (r + 1), les racines étant prises positivement. — Dans 

les cas qui nous intéressent ici, 0< kh/m < V3, donc 1<¢<4 et r=1, 
c’est-a-dire ; ~ 
It 1 

(6) Le produit (A,- surface) est minimal pour un cetcle et vaut alors 

2,40482 =; par conséquent, si les lignes nodales de la k-iéme fonction propre 
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LO‘TE > ¥ > OL‘OE 


CLS = Vi SS 11 


te Crk 
sonbrystiney 


COL > Y > £669 


Up ae 
(VATOg ap Sapoy}9ut) “tp Gg : 


sogyooidde simefea 
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partagent l’hexagone en m, sous-domaines, on a 


2 
= 2,40482 2 


= an 


ny, = 6,9929 n, « 


(ec) A Vaide du quotient de RayLeIcH, PéLrya démontre facilement que 


Les nombres A et 4 indiqués par une astérisque * ont été calculés a partir 
du méme A, qu'une valeur propre approchée précédente; cela tient au fait 
que les fonctions cos et J) ne sont pas biunivoques, et ne se produit jamais 

m 
dans les colonnes de A et A. 

Les valeurs approchées A ont été pour la plupart calculées, pour cet exemple 
précisément, par L. CoLiatz ([5], p. 369); cette méthode classique est, comme 
on voit, bien moins satisfaisante que les trois autres; la formule (22) est du 
méme auteur, et vaut mieux. 


On remarquera les bonnes approximations fournies par J et A notamment 
pour les valeurs propres supérieures, en accord avec la remarque faite au § 1.2. 


12.8 Remarque sur la méthode des distributions et la «formule de la moyenne» 


12.81 On peut retrouver la formule (19) en construisant une distribution (ici: 
fonction) adéquate 4 support compact: on s’inspire du § 2. Prenons %) comme 
origine et introduisons les coordonnées polaires 7, 0; nous définissons 


0 pour" =F: 
"12 GYD OVD —MOYD AVA) poor rs 


ou J, et N, sont les fonctions de BEssEL et de NEUMANN. 
On vérifie que 


AT PAT Sih yA) oer f Du.» 20, 
6=0 


OU do et du, sont des mesures de Dirac: l’une 4 lorigine, l'autre au point de 
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coordonnées polaires h, 8. Donc 
‘Ap +29, T> = <p, AT +1 T> = JhY/2) 90) — 5 f ph, 0) dO, 
WiC 
d’ot. la «formule de la moyenne» (19) pour toute fonction @ satisfaisant 
Ag +A = 0. 


12.82 Si au contraire l’équation donnée est Ag + A = f(x, y), la formule (19) 
est donc remplacée par 


Th V4) 910) — a— f oh, 8) d0 = <T, fy; 
6=0 


de cette formule exacte, on déduit immédiatement une formule approchée 
analogue a (20). 


12.83 Pour l’équation de Poisson Ap = f(x, y) (probleme aux limites), on 
construit la distribution a support compact 


| 0 pour rZh; 
T(r, 6) 1 , 

wee In— .pour nh; 

alors 
1 a5 
4 ERE eer’ Ou, 9) 29 , 
6=0 

d’ou 


1 2% 
(T, fy = ‘Ag, T> = «g, AT> = 90) — az f lh, 8) 48; 
6=0 


on peut en déduire immédiatement la formule approchée 


i " 
ay, 8 HT, fF) + OR"). 


§ 13. Mémes réseaux: plaque vibrante 


L’équation est AAg(x, y) —A p(x, y) = 0; il sera seulement question ici 
de plaques encastrées: yp = 0y/0n = 0 sur le contour. 


13.1 Je décompose l’opérateur différentiel AA-—A=(4+ V2) (A — Va) ; puis 
m’inspirant de la «loi de composition» du § 2.6, je construls le produit de 
composition des opérateurs aux differences correspondant a —3 hh? (A + Va)/2 
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et A —3/? (A — VA) /2. Comme on l’a vu au § 12.1, ceux-ci sont de la forme 


avec 
a ae 
Ay == 2 4 cos( 2 nV) 

pour le premier, et donc 


Ca 
4,=2+4ch(13 2 Vi) 
pour le second, A étant une valeur approchée pour /; posons 
ssi ite 
Y= sh h Vi 


et composons les deux opérateurs discrets; nous obtenons l’opérateur L¥* 
suivant: 


avec maintenant 


A, = (2+ 4 Ch») (24+ 4cosr) +6, B,= —2—4 Chy—4cosp, | 
Sal cheese (23) 
LCRA | 


I est clair que l’équation aux différences fournit une équation transcendante 
pour y, avec les désagréments que cela implique; mais cette circonstance n’est 
pas artificielle, elle est (comme au § 8) dans la nature du probléme (déja pour 
la barre vibrante) 
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13.2 Conditions aux limites: Je me contenterai d’exprimer l’encastrement en 
m’inspirant de la régle (10) du § 1.32 (pour une dimension): soient a, b, c les 
valeurs de u, en trois points intérieurs au domaine (voir la figure ci-dessous) 


’ 


nous définissons #, = a’ en un neud virtuel extérieur par 


ee ee (24) 


I] vaudrait mieux, en principe, s’inspirer de (11) (§ 1.33) [dont (10) est la 
limite pour 7 — 0]; mais (11) entrainerait des calculs trés compliqués. 


13.3 L’opérateur classique aux différences, correspondant a 


as 


ag (44 A), 


est caractérisé par le schéma ci-dessous, avec 
mio = ps 
k 16 ; 


A étant une valeur approchée pour A. 


Méme si l’on utilise cet opérateur classique, on vemplacera avec avantage 
lV’ expression classique de l’encastrement (a’ = a) par la régle (24) ci-dessus: les 
calculs n’en sont pas sensiblement augmentés, mais les résultats sont bien 


meilleurs, comme le montrera l’exemple suivant. 
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13.4 Exemple: Plaque hexagonale encastrée, de coté 1 
(Comparer avec CoLLaTz [5], p. 371.) Calcul approché des deux premiéres 
valeurs propres correspondant a des fonctions propres symétriques. 


Méthode du § 13.1. bs A 
ee : Contour: selonle § 13.2 Moubode aioe Méthode classique 
propre Maille i ; améliorée au contour Cc 
? (équation mere (COLLATz) 
approches transcendante) i 
awe 106 53 32 
h= a 159 136 104 
Ay 2 
— = 158 149 129 
el 468* = zs 
A» h= = 2080 589 524 
sum ' 
WS 3 2237 1397 1150 


L’astérisque * a essentiellement la méme signification qu’au tableau du 
§ 12.7; ici: a partir d’un seul nceud intérieur, on se permet de calculer deux 
valeurs propres approchées; il est clair que c’est trés osé. 

Le calcul de la seconde colonne donne essentiellement le méme travail que 
celui de la troisiéme; l’amélioration est sensible, surtout pour les valeurs ap- 
prochées de 4,; par contre, pour dy, la seconde colonne est presque aussi 

sym 
inexacte que la troisiéme: l’essentiel est ici lopérateur aux différences utilisé, 
plutét que l’expression de l’encastrement. » étant plus grand, Chy et cos» 
sont tres différents des premiers termes de leur développement en série. — Le 
calcul de la premiére colonne est malaisé (équations transcendantes). 


§ 14. Réseaux quadratiques: membrane vibrante 
Nous traduisons cette fois l’équation 
Ag(x, y) + 4 p(x, y) = 0 


approximativement par une équation aux différences pour les valeurs m, de » 
aux noeuds d’un réseau de carrés de cété h. 


l4 dy Soients 4, 0,; bgy%: sie he Ch; aa womee lesevaleuss ide U,; TOUS 
caractérisons l’opérateur aux différences Lj; par le schéma 


HDA ete Fo oe 
7 
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‘c’est-a-dire 


4 
Lt{a] =A,a—35,. 
1 


Nous avons donc 
Caf} \ai— A,@— 57d (25) 


"AS See ai oy) Gan) fee Ad) ee (25’) 


(a) Relation a récurrence: L’idéal serait de combiner linéairement ces deux 
équations en sorte d’obtenir une relation du type 


Om Leal Ag, aac; (26) 


et 


c’est impossible. - Nous voulons sacrifier le moins possible dans le passage de 
(25) et (25’) a (26), c’est-a-dire nous cherchons a déterminer Ag, et des coef- 
ficients « et 6 en sorte que l’expression 


Q = L4,[a] — a L¥[a] — 6 SY L¥[b] soit = O(n") 


avec le plus grand u possible. Pour cela, j’utilise par exemple le développement 
en série donné par Corxatz [5] (tableau VIII, p. 515). 


Q= (A, 2 4, +48) 2+@—-A, fp) Dot B—-1) Yet 2p Ye 
= pl, (Any 4+ (4 A,) (+4 A) + HP AQ|, [-4 +4 + (8— Ay) B+ 0008); 


pour annuler les deux premiers termes du développement de Q, je choisis 


a = (A, — 8) B+4 
(ei 
Aj, 434 (4p 4) + (4-4) 8, (27) 


a 


le coefficient 6 restant /zbre. 
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Remarque: 11 faut le reconnaitre: il y a de l’arbitraire a vouloir annuler 


exactement les deux premiers termes du développement de Q: la récurrence — 


est ici moins satisfaisante que dans le cas du réseau triangulaire (§ 12.1). 
(b) Pour résoudre la récurrence (27), nous l’écrivons, si B + 0, sous la forme 


2 


FEE 2 ! 1=2(Fos 2641) —1; 


la relation cos(2 a) = 2 cos?« — 1 permet donc de poser 
FA, 28 +1 =cos(k h):; 


la constante k reste a déterminer. 
(c) Passage a la limite h > 0: 


Hehe ee — yl, — #2 Ag|, + O(n). 


On retrouve, a la limite, l’équation différentielle 4g + A g=0 a condition 


de poser k = VB aA. L’équation aux différences est donc finalement donnée (en 
fonction du choix de f) par 


& 4,—28 +1=cos(aVBa), (28) 


1 étant une valeur approchée pour A. 
Si au contraire 6 = 0: la récurrence (27) devient A, — 4 = 4 (A, — 4), on 
la résoud en posant A, — 4 = x h®; le passage a la limite 4 > 0 donne 


0= A,a—S"b = xh? g(a) — h® Ag(a) + O(n) ; 


on retrouve, a la limite, Ag + A gm =0 a condition de poser x = —A. On ob- 
tient donc, si 6 = 0, l’ expression classique 


A,=4-Arh, (28’) 
A étant une valeur approchée pour A. 


14.2 Extrapolation du coefficient 8: A cause des termes d’ordre supérieur qu'on 
a négligés dans le développement en série, aucun f ne donne, en général, A=A 
exactement pour toute maille /; on constate cependant qu'il est essentiel, pour 
l’évaluation, de calculer avec un bon coefficient 6: on peut l’extrapoler a partir 
des valeurs A,,, A,,de A pour deux mailles successives h, et hg: on détermine 
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6 en sorte que les deux mailles donnent le méme 1; en d'autres termes, B et-A 
sont déterminés par un systéme de deux équations transcendantes du type (28). 


Utilsation des mailles h et 2h: Il est aisé de calculer f et J a partir de A, 
et A,: la formule de récurrence (27) donne 


butane, 2 12 
pores (29) 


et A est donnée par (28). 


14.3 Fonctions de cellules sur un réseau quadratique 


Un raisonnement sur les fonctions de carrés donne naturellement, pour 
Vopérateur 


les mémes formules (27), (28). — Conditions aux limites: on construira, par ré- 
flexion, des cellules virtuelles extérieures au domaine, en s’inspirant du § 3.42. 

Ainsi que je l’ai remarqué, il y a lieu d’utiliser les fonctions de cellules quand 
la frontiére du domaine considéré présente des comms rentrants. 


14.4 Exemples 


14.41 Carré de coté 1. Fréquence fondamentale: A, = 2? ~ 19,739. 
A laide de fonctions de cellules: 


(28) donne A = 27" exact. 

A Vaide de fonctions de nceuds, on obtient A,j). = 0 et Aig = 2V2, dou 
de nouveau f = 1/2 et Ay = 2 x? exact. 

Comparaison: la formule classique (28’) donne A = 8 pour h ‘a let A = 16 
pour h = 1/2; l’extrapolation classique (admettant que A — A soit proportion- 
nel a h?) donne A ~ 18,67. 
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14.42 Quelques résultats numériques: 


| Méthode 
eta faa 
Valeur propre | “tions | Maille | 4 - (28’) | Extra- 
approximée Fe, (B = 0) | polation 
classique | classique 
t AD 8 
1 ie A, cellules 18,67 
1 0 16 
0 + 
4 | A, | cellules 9,33 
2 8 
—2 6 
1 ay 3 cellules me 11,79 
7 V2 10,34 
h=1] 1 3 
A, ie 1 5,447 
, eh fares 5 |2:7913] 5,100 4,835 
4 : 5,985 6,18 
h = = |3,380 | 5,73 5,58 : 
i) 
1 h= COR ees 2 |3,556702 24 
- : 4,961 77? 47,2 
' |noeuds | 4 = —: V2 |4,74320 41,37 
; 4,998 712 49,18 
Cf. Cottarz [4], h= 3,262 |4,94 7? 47,23 
p. 292 8 
1 * 
1 a i —2 {14,2272 ; = 
; < i 12,54 sr? = 
* | noeuds | 4 = = —//2 |12,10 2? 86,63 
1 12,97670? 126,46 
Cf. Cotxarz [4], h= 8 2,18 |12,82 70° 116,5 
p. 292 


L’astérisque * a le méme sens qu’au § 12.7. 


ZAMP 


5 707 
~ 49,348 
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14.5 Conjecture: Dans les exemples ci-dessus, toutes les valeurs de 6, calculées 
par extrapolation, étaient supérieures ou égales A 1/2; 6 = 1/2 fournissait la 
valeur exacte de la fréquence fondamentale de la membrane carrée, et des 
bornes inférieures Re < A, pour les autres exemples. 

Je conjecture que, lorsqu’on travaille avec des fonctions de cellules (carrés) 


et que l’on pose dans (28) 6 = 1/2, c’est-a-dire A, = 4 cos(h /i/2), ae soit une 
évaluation par défaut pour A,: A, = A). 

J'ai démontré qu’il en est ainsi si, au lieu de procéder par réflexion sur le 
contour: 


Uu 


mais cette recette est inadéquate pour un calcul numérique efficace. 
Il semble bien que, d’une facon générale, la formule 


Ay = 4 c05(h \/2) (¢= 3) 


soit, pour des fonctions de cellules, meilleure que la formule classique 
A, = 4 — A fh? (qui correspond a f = 0). 


§ 15. Réseaux quadratiques: plaque vibrante 
Nous considérons une plaque vibrante encasirée le long de son contour 1”: 
AAp AG s0, == =0, sur J: 


nous travaillons ici avec des fonctions de nceuds. 


15.1 Nous retenons, pour l’opérateur différentiel —h? (A + uw), lopérateur aux 
différences correspondant (28) avec 


ZAMP IXa/12 
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a —h? (A — p) correspond donc l’opérateur discret obtenu en remplacant cos 
par Ch. — Comme au § 13.1, je m’inspire de la «loi de composition» du § 2.6: 
je construis le produit de composition des deux opérateurs aux différences 


considérés, avec “= VA: j’obtiens l’opérateur discret suivant, qui correspond 
a h4 (AA — 2): 


A, = 4-16 Chécosé ; 
by 


= —4(Ché + cosé), (30) 
is E 4 
moos sb 


A étant une valeur approchée pour A. — I] est clair qu’on est amené ainsi a ré- 
soudre une équation transcendante en &. 


avec 


15.2 Conditions aux limites: Comme au 13.2, je me contenterai de transposer 
la régle (10) du § 1.32 pour exprimer l’encastrement: soient a, b, ... les valeurs 
de u, aux noeuds intérieurs du domaine; nous définissons 


Uy, = a = 3a — b/2 en un neud virtuel extérieur au domaine. 


15.3 L’opérateur classique aux différences correspondant a h* (4A — 4) est 
(A étant une valeur approchée pour A): 


De nouveau ici: Méme si l'on utilise cet opérateur classique, il y a avantage 
a remplacer l’expression classique de l’encastrement: 


a’ = a par la régle ci-dessus: a’ = 3. a — b/2. 


Se 
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15.4 Exemple: Plaque carrée encastrée, de cété 1 


ARONSZAJN et DoNoGHUE (([1], part I, p. 34) ont démontré: 


120318 SY) 13003 


nie ; : ; 
Voici les valeurs approchées fournies pour A, par les trois méthodes: 


Méthode du § 15. 1; 
BS A513 Méthode classique 


, contour: selon le § 15.2 : 
Maille Beatie cee améliorée au contour Méthode classique 
2 dante) selon 15.2 
Ee ae ee ee ee ee ee ee 
il 
A= oF 791 512 384 
1 
h= cw 1076 891 648 
1 
h= vs 1195 1068 829 
ai 
h= = 1245 1154 949 


Notre expression de l’encastrement joue un role essentiel, bien qu’elle ne com- 
plique pas les calculs. 


[1] 


[8] 
[9] 
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Summary 


Let a linear boundary value problem or an eigenvalue problem be given, in 
one or more dimensions. I look for such a construction of a corresponding dif- 
ference equation that, in the case of ordinary differential equations with constant 
coetficients, its solutions be exactly those of the continuous problem. — This 
‘trivial’ case is first investigated and several methods are pointed out; these are 
then applied to: (a) ordinary differential equations with variable coefficients 
(‘method of local perturbations’); (b) vibration problems for membranes and 
clamped plates. 


(Recu: le 10 septembre 1957.) 


Kennzeichen atmospharischer Eispartikeln 
i, Der 


Graupeln als Wachstumszentren von Hagelkérnern 


Von Ro.Lanp List, Weissfluhjoch-Davos!) 


1. Einleitung 


Beobachtungen und Messungen an iiber 60 Sammelproben von Graupel- 
und Hagelfallen der Jahre 1953 bis 1957 sollen Aufschluss geben iiber das 
Wachstum kleiner Eispartikeln zu grossen Hagelkérnern. Die verschiedenen 
Entwicklungsstufen, ausgehend vom Eisbildungskern, sind gemiss dem Schema 
in Figur 7 in folgende Haupt- und Nebenphasen zu unterteilen: 

A. Wachstum durch Sublimation: vom Eisbildungskern zum Schneekristall. 
B. Wachstum durch Anlagerung von Wassertropfen: vom Schneekristall zum 

Hagelkorn. 

a) vom Schneekristall zur Reifgraupel; 

b) von der Reifgraupel zur Frostgraupel oder Klareispartikel: 

c) von der Klareispartikel oder Frostgraupel zum Hagelkorn. 


1) Bidgendéssisches Institut fiir Schnec- und Lawinenforschung. 
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Phase A, Ba) und Bc) stellen dabei ein eigentliches Gréssenwachstum dar, 
Phase Bb) besteht in erster Linie in einer Verdichtung. 

Wahrenddem viele Arbeiten iiber Eiskerne ausgefiihrt, zum Beispiel [1-3]?), 
und Wachstum und Aufbau der Schneekristalle umfassend durch NAKAYA [4] 


Eisbildungskern 


Schneekristall 


Sublimation 


schnelles Gefrieren 


Gréssenwachstum 


Reifgraupel 


schmelzen teilweises Schmelzen langsames Gefrieren 


Wassertropfen 


Verdichtung 


q 
oO 
n 
o 
aq 
Qy 
wn 
S 
= 
p 
n 
Gq 
ie} 
o 
3 
~ 
OQ 
=) 
a 
q 


Wachstumsweise durch 


Gefrieren Gefrieren 


Klareispartikel TypII Frostgraupel Klareispartikel Typ I 


Anlagerung von Wassertropfen 


langsames oder schnelles Gefrieren 


ie | ae” 


Hagelkorn 


Gréssenwachstum 


Figur 1 
Wachstum vom Eiskern zum Hagelkorn. 


behandelt wurden, sind in den Kenntnissen tiber gréssere atmospharische Eis- 
kérner diverse Liicken festzustellen. Das Ziel der vorliegenden Arbeit ist es 
nun, diesem Mangel zumindest teilweise abzuhelfen. 

Bei der Charakterisierung von Eispartikeln traten in der Folge verschiedent- 
lich Differenzen auf mit der Nomenklatur, wie sie von der «World Meteoro- 
logical Organisation» vorgeschlagen wird. Der Grund liegt darin, weil neues 
Beobachtungsmaterial solche Abweichungen wiinschenswert erscheinen lasst. 


2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 192. 
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2. Reifgraupeln 
2.1 Charakteristerung 


Aussehen: Reifgraupeln («soft hail, graupel or snow pellets») sind weisse, 
opake Partikeln kugelsektorartiger (kurz konischer), rundlicher oder unregel- 
missiger Form mit einem Durchmesser, der 7 mm selten tiberschreitet. 

Entstehung: Reifgraupeln entstehen durch Anlagerung von Wolkentropfen 
an einen urspriinglichen, durch Sublimation gewachsenen Eiskristall. Das An- 
gefrieren des eingefangenen Wassers erfolgt relativ rasch am Ort des Auf- 
treffens (vergleiche auch mit Mason [3)). 

(Die deutschsprachige Bezeichnung «Reifgraupel» ist insofern irrefiihrend, 
als es sich nicht um durch Reifansatz, das heisst durch Sublimation gewach- 
sene Teilchen handelt — eine Erkenntnis, die erst neueren Datums ist. Der 
Name ist damit nur in historischem Sinne als Bezeichnung fiir spezielle Eis- 
partikeln aufzufassen.) 


2.2 Ausserer Aufbau der Reifgraupeln 


2.21 Entstehungsphasen und Formen. Bei der Auswertung der Graupelfalle 
vom 28.Juli und 2. September 1957 zeigten sich Teilchen verschiedener 
Wachsstumsstadien, so die in Figuren 2 und 3 dargestellten, mit gefrorenen 


se TA ; ; 
Figur 2 Figur 3 


Vergraupelter Schneestern Nr. 1. Vergraupelter Schneestern Nr. 2. 


Wassertropfchen bedeckten Schneesterne oder eigentliche konische Reif- 
graupeln, wie sie in Figuren 4 und 8 zum Ausdruck kommen. Dass die Form 
mit der Partikelgrésse wesentlich zusammenhingt, zeigt die Tatsache, dass 
praktisch alle Teilchen mit einem Durchmesser von mehr als 2-3 mm kugel- 
sektorartig befunden wurden. Der Grund hierzu diirfte darin liegen, dass bei 
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den kleineren Reifgraupeln teilweise noch die Struktur des ursprtinglichen ’ 
Schneekristalles zu erkennen ist, die durch die aerodynamisch bedingte Ak- 
kreszenz beim Fallen nach und nach verwischt wird. Die damit entstehenden 
grosseren Teilchen weisen daher meist symmetrische Formen auf, wie sie auch 
durch Winkelmessungen erhartet wurden. Kegel6ffnungswinkel an 80 koni- 
schen Graupeln vom 23.August 1954 ergaben Werte zwischen 80 und SOS. 
ebensogross waren die Ubergangswinkel zwischen Kegel- und Kugelflacheu. 


Figur 4 
Kugelsektorartige Reifgraupeln, Achse senkrecht zur Bildebene. 


Wahrend letztere bei weiteren Beobachtungen konstant blieben, betrug der 
mittlere Kegeléffnungswinkel bei Sammelproben vom 28. Juli 1957 rund 10. 
L. v. Bucu [5] stellte bereits 1813 fest, dass konische Reifgraupeln mit der 
Spitze nach hinten fallen, eine Ansicht, die von ARENBERG [6] nicht geteilt 
wurde. Spiter angefiihrte Photos von Hagelkérnern und Frostgraupeln zeigen 
jedoch die Berechtigung der alteren Behauptung. Da bei den vorkommenden 
' Luftgeschwindigkeiten eine Vereisung immer gegen den Strom stattfindet, er- 
gibt sich die Aussage, dass der primare Eiskristall in der Graupelspitze zu 
suchen ist. Die Figuren 5 und 6 bestatigen dies, wird doch dort bei Scharfein- 
stellung des Mikroskopes auf eine Ebene senkrecht zur Kegelachse vielfach 
eine Andeutung einer dendritischen Struktur gefunden. Diese speziellen Kenn- 
zeichen bleiben manchmal teilweise tiber die ganzen Graupelkegel erhalten, 
wie es Unregelmassigkeiten der dieselben abschliessenden «Kreise» zeigen (zum 
Beispiel Figur 4). : : 

Sind die Wachstumszentren keine eigentlichen Schneesterne, so wird es 
kaum mehr méglich, den Typ des Ausgangskristalls zu eruieren. Dies spielt 
speziell beim Versuch, auch auf Grund der Wachstumsbedingungen der ver- 
schiedenen Schneekristallarten auf die Ausgangsbedingungen des Reifgraupel- 
Wachstums zu schliessen, eine Rolle. Die hierfiir von NAkaya [4] angegebene 
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Figur 5 Figur 6 
Graupelspitze mit Andeutungeines Schneesterns. Graupelspitze mit Andeutung eines Schneesternes. 


Methode befriedigt jedoch nicht, da die angefiihrten Daten den Tatsachen 
nicht immer zu entsprechen scheinen. (NAKAYA fordert fiir das Auftreten von 
Dendriten Temperaturen zwischen — 13 und — 18°C; da aber wahrend der in 
Frage kommenden Graupelfalle bei einer Wolkenhéhe iiber der Station von 
maximal 200 m Bodentemperaturen zwischen + 1 und + 2°C gemessen wurden, 
scheint dies ganz unméglich.) 

2.22 Gréssenverteilungen. Als Beispiele von Gréssenverteilungen seien in 
Figur 7 drei Graupelfalle dargestellt. 
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Grossenverteilung dreier Reifgraupelfalle. 


Kurve a betrifft die Sammelprobe vom 19. August 1954, 16.25 Uhr, bei 
total 548 Partikeln. 


Kurve 6 betrifft die Sammelprobe vom 23. August 1954, 13.30 Uhr, bei 
total 4912 Partikeln. 


Kurve c betrifft die Sammelprobe vom 23.August 1954, 15.45 Uhr, bei 
total 1226 Partikeln. 
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Figur 8 Figur 9 


Diinnschnitt einer Reifgraupel. Diinnschnitt einer Reifgraupel. 


2.23 Grosse der Aufbauelemente. Die Grésse der Aufbauelemente, das heisst 
der angelagerten Wolkentropfen, ist aus den friiher erwahnten Figuren 2 und 3 
ersichtlich. Ihr Durchmesser liegt im Mittel um 20. Ein merkbarer Gréssen- 
unterschied zwischen direkt an den urspriinglichen Kristall gelagerten Teil- 
chen und den im Endstadium angewachsenen Partikeln besteht, wie eine Aus- 
messung verschiedener, rund 0,35 mm dicker Diinnschnitte (Figuren 8 und 9) 
ergab, nicht. Spezielle Angaben tiber Diinnschnittechnik siehe [7]. 


2.24 Dichten von Reifgraupeln. Die Diinnschnitte geben in durchscheinen- 
dem Licht Aufschluss iiber die in Eispartikeln eingebauten Luftblasen; sind 
diese gegenseitig so gekoppelt wie die bei Reifgraupeln auftretenden Kapillar- 
systeme, so kénnen die Eiskérner direkt als eigentliche Eisgeriiste angespro- 
chen werden. 

Eine direkte Dichtemessung durch Wagung und Bestimmung des Auf- 
triebes in einer Fliissigkeit wird damit praktisch illusorisch. So wurde zur Er- 
mittlung des Volumens geometrisch einfacher Graupeln eine Ausmessung vor- 
gezogen. Die berechneten Dichten betrugen fiir den Graupelfall vom 28. Juli 
1957 im Mittel 0,7 g/cm, fiir die Partikeln vom 3. September 1957 rund 
0,5 g/cm3. Diese Werte, die mit einem Fehler von + 20% behaftet sind, liegen 
im allgemeinen hoher als die von andern Autoren angefiihrten Werte. Na- 
KAYA [4] zum Beispiel gibt in seinem Standardwerk den allgemeinen und kon- 
stanten Wert von @ = 0,125 g/cm® an und lasst dabei die Abhangigkeit des 
spezifischen Gewichtes von der Packung der Wolkentrépfchen unberiicksich- 
tigt, die ihrerseits von der Gefriergeschwindigkeit des an ein Eispartikel an- 


gelagerten Wassers beeinflusst wird. 


2.25 Zonen gleicher Wachstumsbedingungen. Eine starkere Anderung der 
Wachstumsbedingungen von Reifgraupeln hat eine andere Dichte zur Folge, 
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so dass zusammenhangende Zonen gleicher Eiskonzentration als Wachstums- 
phasen angesprochen werden kénnen. Ansatze hierzu zeigt Figur 9. Die darin 
abgezeichneten Kugelschalen-Ausschnitte sind gestaffelt beziiglich der Grau- 
pelspitze, die auch damit als Wachstumsausgangszone oder Wachstumszentrum 
angesprochen werden darf. 

Diese Kennzeichen sind aber speziell bei Vergleichen mit Bedacht zu inter- 
pretieren, da bei gleichen Umweltsbedingungen das Wachstum von Eiszonen 
von der Fallgeschwindigkeit der betrachteten Kérner abhangt! 


2.3 Kristallographie von Reifgraupeln 


2.31 Kennzeichen, die auf die Richtung des értlichen Wachstums deuten. Das 
Interesse an der Kristallographie natiirlicher Eispartikeln basiert auf Hin- 
weisen, die in polarisiertem Licht beziiglich deren Wachstum sichtbar werden 
kénnten. Diinnschnitte zeigen denn auch die Anordnung der zu beobachtenden 
Eiseinkristalle — bei Reifgraupeln sind darunter eher Zonen gleicher Orientie- 
rung zu verstehen (Figuren 10 und 17). Irgendwelche Gesetzmassigkeiten tiber 


Figur 10 Figur 11 
Dunnschnitt durch die Achse einer Reifgraupel, Dinnschnitt senkrecht zur Achse einer 
mit polarisiertem Licht. Reifgraupel, mit polarisiertem Licht. 


Wachstumszonen und Wachstumsrichtungen sind aber infolge der geringen 
Anzahl von Kristalliten in Reifgraupeln nicht festzustellen, wohl aber, wie 
spater gezeigt wird, bei grésseren Teilchen. 


2.32 Anzahl der wirksamen Eisbildungskerne. Unter der Annahme, dass ein 
einzelner Einkristall nur bei Anwesenheit eines Eiskernes entstehen kann, ist 
der Schluss nahe, dass die Zahl der Einkristalle gleichgesetzt werden kann mit 
der Zahl der wirksamen Sublimations- oder Gefrierkerne. Es gilt damit, die 
Zahl der Kristallite in einem Reifgraupel aus verschiedenen Diinnschnitten in 
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der Kornachse (zum Beispiel Figur 70) oder senkrecht dazu (Figur 17) zu 
_ interpolieren. Das ausgewertete Material lasst beziiglich dem Graupelfall vom 
_ 28. Juli 1957 ein Einkristallvolumen von rund 1-5 mm? erwarten, mit andern 
_ Worten kommen bei einer gemessenen Eisteilchendichte von 0,7 g/cm® pro 
Gramm Wasser 500 Eisbildungskerne. Daraus kénnen wir auch die Eiskern- 
_ konzentration der Erzeugerwolke errechnen, die unter der Annahme eines 
_ Wassergehaltes von 3 g/m? rund 1,5 wirksame Teilchen pro Liter Luft ergibt. 


3. Frostgraupeln 
3.1 Charakterisierung 


Aussehen: Frostgraupeln («small hail») sind Partikel kugeliger, kegelartiger 
oder unregelmdssiger Form. Ihr Aussehen ist teilweise weiss und opak; die 
- Oberfliche ist immer zonenweise glasartig. Der Durchmesser iiberschreitet 
_ selten 7 mm. 


Entstehung: Frostgraupeln entstehen durch Anlagerung von Wassertropfen 
an Reifgraupeln; der Gefrierprozess findet langsam statt, so dass das ange- 
_ schwebte Wasser Zeit hat zu zerfliessen. 


3.2 Ausserer Aufbau von Reifgraupeln 


3.21 Entstehungsstadien. Ein Vergleich mit der Charakterisierung des 
Wachstums von Reif- und Frostgraupeln zeigt, dass der Unterschied in der 
_ Gefriergeschwindigkeit liegt: Bei grosser Unterkiihlung des angestrémten 
_ Wassers erfolgt ein schnelles Gefrieren direkt am Ort der Anlagerung, wobei 
die Tropfen ihre Form mehr oder weniger bewahren; langsames Gefrieren bei 
geringer Unterkiihlung lasst jedoch dem Wasser Zeit, sich zu verteilen. Eine 
_scharfe Trennung der Graupelarten kann trotz aller Zwischenphasen ausge- 
fiihrt werden, wenn die augenfallige zonenweise Verglasung der Oberflache bei 
teilweisem opakem Innern als Hauptcharakteristikum fiir Frostgraupeln an- 
genommen wird. Der Ubergang einer Reifengraupel in eine Frostgraupel sei 
anhand von Figur 12 erklart, wobei zu bemerken ist, dass das Korn streng 
kugelsektorartig war, der Schnitt jedoch nicht in der Kegelachse gelang. — Die 
Ubergangsgraupel stellt ein Wachstumsstadium dar, wo der Habitus noch auf 
eine Reifgraupel schliessen lasst; die gut sichtbare teilweise Verglasung sowie 
der Diinnschnitt zeigen jedoch, dass die urspriingliche Reifpartikel in warmere 
Zonen mit weniger stark unterkiihlten Wassertropfen gelangt sein muss, da die 
Kapillaren von der Gegenstromseite, vom kugeligen Flachenteil her zum Teil 
ausgefiillt worden sind. Dafiir spricht die deutlich zu erkennende Dichte- 
abnahme gegen die Spitze hin. - Nach Beginn der Verdichtung der Partikel 
bei konstantem Volumen erfolgte ein weiteres Gréssenwachstum an der Stirn- 
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Figur 12 Figur 13 
« Ubergangsgraupel». Frostgraupel. © 


seite, wo eine blasenarme Schale auf Frostvergraupelung hinweist. In einem 
weiteren Stadium wiirde ein Zerfliessen der angelagerten Wassertropfen in die 
Kegelflache erfolgen, gekoppelt mit einem Vollsaugen der Kapillaren auch von 
diesen Seiten her. Die so entstehende Form lasst sich durch einen einer Kugel 
aufgesetzten Kegel charakterisieren. 

Urspriinglich kugelige oder unregelmassige Reifgraupeln behalten ihre 
Form beim Ubergang zu Frostgraupeln bei. 


3.22 Grdssen- und Formvertetlungen. Das Untersuchen von Frostgraupeln 
wird erschwert durch ihre nassen Oberflachen, die beim Sammeln ein gegen- 
seitiges Verpappen der Korner und Diinnschnittbilder gemass Figur 73 zur 
Folge haben. Beim Graupelfall vom 18. Juli 1954 auf Weissfluhjoch waren die 
Partikeln jedoch relativ trocken, so dass es ohne Miihe gelang, Gréssen- und 
Formverteilungen abzuklaren. Diese sind in Figur 74 dargestellt und zeigen 
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“ zwei Maxima bei 4,0 und 6,5 mm grésstem Durchmesser, eine Eigenheit, die 
_ von zwei zeitlich etwas gestaffelten Niederschlagen herrithren kann. In der 
bi gleichen Darstellung ist auch der Anteil der kugeligen, kegeligen und unregel- 
_ massigen Formen festgehalten. 

a Eine Unterscheidung nach der Durchsichtigkeit erfolgt fiir kugelige und 
% kegelige Partikeln in den Figuren 15 und 16, wobei fiir konische Frostgraupeln 
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Figur 16 Figur 17 
Durchsichtigkeitsverteilung kegeliger Kegeloffnungswinkelverteilung kegeliger 
Frostgraupeln. Frostgraupeln. 


der Sitz der opaken Zonen jeweils in der Gegend der Spitze zu suchen ist. Die 
Maxima der undurchsichtigen Teilchen liegen jeweils bei den kleinern Gréssen, 
wahrenddem die halbdurchsichtigen eine Haufung bei grossen Durchmessern 
aufweisen. Bemerkenswert ist weiter, dass bei den kugeligen Frostgraupeln 
vereinzelt ganz durchsichtige, sogenannte Klareispartikeln festgestellt werden 
konnten. Figur 17 gibt uns weiter Aufschluss iiber die Offnungswinkelvertei- 
lung fiir die kegeligen Formen. Die Streuung in einem Bereich von 60° ist er- 
wartungsgemiss wesentlich grésser als bei den untersuchten Reifgraupelfallen 
mit 10°. Der Grund hierfiir diirfte darin legen, dass der Kegelwinkel nicht nur 
von den Wachstumsbedingungen, sondern speziell auch vom Wachstums- 
stadium abhangt. Die in der gleichen Figur vorgenommene Unterteilung nach 
Durchsichtigkeitsgrad gibt keine speziellen Hinweise. 


3.23 Dichten von Frostgraupeln. Die vorgenommenen Dichtemessungen er- 
gaben Werte zwischen rund 0,8 und 0,9 g/cm®; sie liegen damit sehr nahe bei 
der Eisdichte, die 0,915 g/cm? betragt [8]. 


3.3 Kristallographie von Frostgraupeln 


3.31 Wachstum der Einkristalle. Ein Teil der Reifgraupeln vom 28. Juli 
1957 wurde nicht direkt in gekiihlte Gefasse aufgefangen, sondern vom war- 
men, nassen Boden aufgenommen, wo eine Benetzung durch Ansaugen von 
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~ Figur 18 Figur 19 
Kiinstliche Frostgraupel, mit polarisiertem Licht. Frostgraupel, mit polarisiertem Licht. 


fliissigem Wasser in das Eisgeriist hinein stattfand. Das Resultat war eine Art 
kiinstlicher Frostgraupeln (Figur 78), die von natiirlichen (Figur 79) nicht zu 
unterscheiden sind. Die beobachtete mittlere Einkristallgrésse entsprach den- 
jenigen der Reifgraupeln gleichen Datums. Dies bedeutet, dass das vom Eis- 
geriist aufgenommene Wasser beim Gefrieren die Kristallrichtung der Um- 
gebung annahm. Zum gleichen Resultat fiihrten auch Diinnschnitte, bei denen 
Reifgraupeln durch Auftauen und Angefrieren an Objekttrager zubereitet wur- 
den, anstatt wie normal mit der anderswo beschriebenen Phtalsduretechnik [7]. 

Die Einkristallanordnung des Ubergangsgraupels (Figur 20) zeigt weiter, 
dass vereinzelt eine Ausdehnung der Kristallite durch verschiedene Wachs- 
tumsphasen hindurch modglich ist. 


3.32 Anzahl der wirksamen Eishildungskerne. Die Anzahl der wirksamen 
Eisbildungskerne bleibt auf Grund der oben gemachten Angaben ungefahr 
gleich wie bei den urspriinglich gleich grossen Reifgraupeln (350/cm%). Die spe- 
zifische Anzahl pro Liter Wolkenluft wird aber zufolge der grésseren Teilchen- 
dichte wesentlich reduziert. Dies ist leicht erklarlich, da das in die Kapillar- 
systeme eindringende Wasser die beste Gefrierkernauslésung antrifft: Das vor- 
handene Eisgeriist, dessen Einsatzpunkt bei 0°C liegt. Demgegeniiber bendtigen 
jaalle bekannten Kernsubstanzen Unterkihlungen von mindestens — 4°C [1, 2]. 

Ein Vergleich verschiedener natiirlicher Proben ist in gréssenordnungs- 
missiger Ubereinstimmung mit den oben diskutierten Werten. Dabei muss 
aber klar sein, dass eine gréssere Streuung erwartet und gefordert werden 
muss! Die Unterkithlungsfahigkeit der Wolkentropfen, die von der Anwesen- 
heit und der Art der Fremdsubstanzen abhangt, beeinflusst die Wachstums- 
bedingungen atmospharischer Eisteilchen in ausserordentlichem Masse und 


schafft erst durch ihre Variation die Vielfalt des Aufbaues natiirlicher Graupel- 
und Hagelkérner. 
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4. Klareispartikeln 


Die oni eigenstandig oder mit Frostgraupeln auftretenden Klareis- 
- partikel machen eine weitere Spezifizierung ndtig: 


Figur 20 
Ubergangsgraupel, mit polarisiertem Licht. 


4.1 Klareispartikeln Typ I 


Aussehen: Klareispartikeln Typ I sind klare, kugelige Teilchen mit verein- 
_zelt ins Eis eingebauten Luftblasen. 

Entstehung: Sie entstehen ahnlich wie Frostgraupeln, mit denen sie beob- 
achtet werden kénnen. Erforderlich sind neben spezifischen Gefriergeschwin- 
digkeiten ganz zusammenhangende Kapillarsysteme der urspriinglichen Reif- 
graupeln, die es der Luft erlauben, beim Eintreten von Wasser bei der Frost- 
Vergraupelung praktisch komplett zu entweichen. 

Kristallaufbau: Klareispartikeln Typ I sind meist polykristallin. 


4.2 Klareispartikeln Typ II 


Aussehen: Klareispartikeln Typ II («ice pellets») sind ebenfalls klare, 
kugelige Teilchen mit vereinzelt ins Eis eingebauten Luftblasen. 

Entstehung: Sie bilden sich durch Gefrieren von Wassertropfen, was zur 
Folge haben kann, dass das Innere unter Umstanden aus fliissigem Wasser be- 
stehen kann. 

Kristallaufbau: Klareispartikeln Typ II kénnen mono- oder polykristallin 


sein. 
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4.3 Klareispartikeln Typ III 


Aussehen: Klareispartikeln Typ III sind klare, teilweise oder ganz eben- 
flachige und kantige Kristalle mit einzelnen eingebauten Luftblasen. 

Entstehung: Typ III wurde mit Frostgraupeln zusammen beobachtet. 
Trotzdem ist seine Entstehung unklar. 

Kristallaufbau: Diese kristallartigen Klareispartikeln Typ III kénnen 
mono- oder polykristallin sein. 

(Beobachtungen an Klareispartikeln sind bereits in der Mitteilung Zur 
Struktur von Hagelkérnern von R. List und M. DE QUERVAIN [9] publiziert 
worden.) 

Fir die Unterstiitzung und Foérderung dieser Arbeit méchte ich vor allem 
Herrn Dr. M. DE QUERVAIN sowie auch Herrn Dr. Tu. Zincce danken. Die 
Untersuchungen wurden im Rahmen des Forschungsprogrammes der Eidge- 
néssischen Kommission zum Studium der Hagelbildung und der Hagelabwehr 
an der Forschungsstelle Weissfluhjoch durchgefiihrt ; zur Finanzierung standen 
teilweise Kredite des Schweizerischen Nationalfonds zur Verfiigung. 
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Summary 


Observations and measurements of more than 15 samples of soft hail and 
small hail disclose the nature and the growing process of these particles. 


(Eingegangen: 21. Januar 1958.) 
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Normal Modes Methods for Boundary-Excited Wave Guides 
By JEAN VAN BLADEL, Madison, Wisconsin, USA?) 


1. Introduction. The electromagnetic field within a cavity is uniquely deter- 
mined by the current density /(x, y, z) in the cavity volume, and the tangential 


component of the electric field E at the boundary surface S (Figure 1). This 


E tangential 


E 


normal 


Figure 1 


Cavity with hole in the wall. 


component vanishes along those parts of S which consist of a perfectly con- 
ducting wall, and is different from zero across the surface of any apertures which 
couple the cavity with outside space. Explicit formulas giving the cavity field 
in terms of ‘volume’ excitation J and ‘boundary’ excitation E,,,,, are impor- 


tant in a number of problems, such as the determination of the resonant 


frequencies of coupled cavities. Formulas of that type can be obtained by 
using either GREEN’S dyadic ({1], p. 1782)?) or normal mode [2] methods. The 
two formulations are essentially equivalent, since GREEN’s dyadic can be 
expressed in terms of the normal modes of the cavity ({1], p. 1777). In the 
latter form, which is used in most practical applications, GREEN’s dyadic 
appears to be a mere intermediate step in the normal mode formulation. A 
similar situation arises in the study of boundary excited wave guides. Here 
again, GREEN’s dyadic can be used in principle, but its analytical expression 
always involves the normal modes of the duct ({1], p. 1825). It is the purpose 
of the present paper to solve the boundary excitation problem directly by 
normal mode procedures. It is the author’s experience that such procedures 
are eminently suitable for practical computations. The formalism is clearer, 
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the probability of mixing up coefficients is lower, and a very desirable flexi- 
bility is afforded by the independence of the method on the type of boundary 
conditions existing at the terminal planes S, and S, (Figure 2). In the mathe- 


Figure 2 


Boundary-excited wave guide. 


matical treatment which follows, volume excitation J is added for reasons of 
completeness, although its effects are quite satisfactorily analyzed in the existing 
literature. 


2. Consider the following two-dimensional problems 


: LC eee 
a) Vi,9 + Bo = 55 + ayet Mp =0, | 


(2g 
with » = 0 on the contour C of the wave guide, | 


(b) Vey yp + ky =0 with su =0 on C (nm isthe outward normal toC). (2.2) 


The solutions of (2.1) form a denumerable infinite closed and complete set of 
orthogonal eigenfunctions, which will be denoted’) by E,,,(x,y) and the 
corresponding eigenvalues by k?,. Another such set is formed by the solutions 
of (2.2), denoted by H,,,(x, y), with eigenvalues k?. The eigenfunctions E,,, and 
H,, will be used to expand the z-component of the electric and magnetic fields 
respectively. 


(c) Ve V = 0, with V = const on each closed curve of boundary C. (2.3) 


Equation (2.3) has 0,1, —m linearly independent solutions depending on 
whether the wave-guide space is simply, doubly, ... , (# + 1)-fold multiply- 
connected. The solutions of (2.3) will be denoted by V5. 


3) The subscripts m and m stand for a double set of indices j, k, not written explicitly for 
reasons of brevity. 
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Consider now the following two sets of vectors. 
Set 1 1.1 grad,, E,,, , denoted by E,,, (x, y). 
1.2 u, Xx grad,,H,, , denoted by E,, (x, y). 
A 1.3 grad,, Vo, denoted by E,(x, y) 
Set 2 2.1 grad,, H,,,, denoted by H,,(x, y) . 
2.2 u, X grad,, E,m , denoted by Hyya(“, y) . 


2.3 u, x grad Ys, denoted by Hyo(x, ¥) . 
The first set comprises the mutually orthogonal‘) eigenvectors of equation 
V5 A+k®?A=0 (A perpendicular to C, div,, A=0onC). (2.4) 
The second set comprises the mutually orthogonal eigenvectors of 
ees RB=0 (B tangent tai, (dive; B) = 0 on c) (2h 


It will be noticed that Eig H,,, , E,, Hy are irrotational, EE H,,, : Ee H,, SO- 
lenoidal®). Each set of vectors forms a closed set in terms of which a sufficiently 
‘smooth’ two-dimensional vector can be expanded. Set 1 will be used to 


expand Ey transverse component of the electric field, set 2 to expand H,. 
More explicitly: 


E(x,y,z,t) = Di nl (2, t) Evy, (x, 9) +B bal (2, t) Ey, (x, 9) 


siz D ital (z, t) ) Exo( x, ¥) = dD) Omi, t) Ezm(%, 9) Uz, 
0 m 


H(x, y, z,t) = Dy %ml2t Han(*, 9) +S Bre H,,,(x, 9) 


Fd, tol, # H,,(x, ¥) + 2, Yal@ 1) Healt, 9) te. 


(2.7) 


4) The orthogonality property is of the form /| Eig? Ey 4S = Op. 
s 
5) The m, n and o functions are, except for multiplicative constants, the components of the 
TE, TM and TEM modes of classical wave-guide theory. The multiplicative constants have been 
dropped for reasons of clarity. They will be reintroduced automatically in the sequel [see equation 


(3. 6)]. 


196 JEAN VAN BLADEL ZAMP 


The equality sign must, as usual, be understood in the sense of convergence in 

the mean. The problem is to determine the differential equations satisfied by — 
the various (z, t)-dependent coefficients. These equations are listed below 
without proof, a typical proof being given in the appendix. 


(a) For the O Modes 


da da 1 av, 
Ber xg_| one wg | Btn “de, (2.8) 
Cc 
Ox, da te 
me tee = — ag ff] FB dS (2.9) 
(b) For the M Modes 
0am Oem il = = 
pele — Om == 9g | (in X EB)» Hiyy do = vag f B ey eis. 
Cc 
Om 04m 1 ffs 
st +e = | Bm aS. (2.11) 
S 
un hope La 
a ee nz || J: Een aS. (2.12) 
S 
(c) For the N Modes 
OBn Ob, 1 fos 
oe errs tet hears HT -Em4S, (2.13) 
a 
Ob, OB, 1 "7 >\ 
Oz U “Ot a Ne / (u,, x E) - #7, tn dc = ne Be ae dc, (2.14) 
Cc 
ba ene Soe eT 
Be EST ory Ae ae Gee sade Nol ar Bie ges 2-43) 
Cc n C 


Some of the notation should be clarified: 


y= | Eolas, 
s 


N= |f Bini? 4S = ff] Hoy aS = #2, ff BR, a5, 


pe) S s 
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= [[\Bn las = [[1H,,8as = 2 [HB as, 
iS S s 


j Ses En, U U, H.,, =H,,, Uz. 


U, = unit vector pointing out of the guide. (w,, %,, u,) form a right-handed 
frame. 

The second members of equations (2.8) to (2.15) are known functions of z 
and #. One readily recognizes the formal analogy between the volume electric 
current re which couples to the E field of the mode, and the fictive magnetic 


surface current “, x E, which couples to the H field of the moe 


3. At this point the vectorial field problem has been reduced to a set of three 
independent systems of partial differential equations in z and ¢. The solution 
of the latter can now proceed, given suitable initial conditions at some time 
#=0, and boundary conditions*) at terminal planes S, and S,. The usual 
arsenal of mathematical techniques is available for the purpose. Laplace and 
Fourier transforms can be used to eliminate the time dependence in the 
equations. Laplace and Fourier transforms, or finite Fourier transforms in z 
are helpful for guides of respectively infinite, semi-infinite or finite length. It is 
advisable to start from equations involving one function only. Such equations, 
which can be obtained by elimination, are: 


oe op em ae (is a Ir E,.45 — 5 | (iin E)-Hy de], (3.1) 
S Cc 


dz? N? 
2 De oes 
— 159) ee [esr | (iin x py. Hj, de — ~ IJ: E45}, (3.2) 
0 é é 
02b, OM pa 1 ak a 
pa 8B gp Pn On = HR eal Dore? | 
(3.3) 
oma Hy de— Hef (Ba B) HE nde], | 
07% O° am, 2 i! 0 (tip X E)-H, dc 
5 peel arr Bin Sm = N le or | “rn me 
c (3.4) 
— ff F Ey dS + ff J > Bam de] | 


S S 


8) Example: yp; 4m; %, On Vanish at all times on a perfectly conducting terminal plane. 
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Equations (2.8) to (3.4) contain all the familiar [3] results for the normal modes 
in a region free of volume and boundary excitation. For an ‘m’ mode in sinus- 
oidal time-dependence, for instance, equations (2.10) to (2.12) show that ¢m 
must satisfy 
—™ + (wre w— ki) a ae Vin &m = 0. (3.5) 


Ie ake Ripe ees 

ire imi nes (3.6) 
JWeE Dee ae 

bn = i, C= 75 e 1%? | (3) 


yielding the usual proportionality relations between the mode components of a 
progressive wave. 


dol 


role 


Figure 3 


Gap-excited accelerator: (a) top view, (b) side view, (c) cross-sectional view. 


4. To illustrate the preceding analysis by a practical example, consider an 
infinitely long rectangular wave-guide, cut in two by a plane making an angle 0 
with the guide’s longitudinal axis (see Figure 3). The gap is very narrow, and 
an alternating voltage of constant phase and amplitude is applied across it 
Such a physical system can be encountered in particle accelerators, for ingvance 


The fields inside the wave-guide can be calculated by setting J=0 (no volume 
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_ excitation) and 


U, X E=V coswt d(p) (tn X ty) 


in equations (2.8) to (3.4). If one assumes that the frequency of operation is 
less than the cut-off frequency c/2 a of the guide, the various z-dependent 
coefficients %», b,, etc., are required to vanish for z = + oo. If this is taken 


' into account, and noticing that 


_ mI ON Seo Tg b é 
E,, = sin 7 (x+ 7) sin ; (y + 5) (m and n integers) , 


H,,, = cos («+ 5) cos (y+ 7] 


_ for a rectangular wave-guide, the following expressions are obtained for the 


y-component of the magnetic field: 


(a) Within the Gap Region (0 Sz Sa tg 6/2). 


(m-+-n)/2 
WV aes ne = ees _ Max NEY — 
es = =. - 0 sinwt-y Da ae —sin : cos . Me e Sm,n ™ tg6/2 
‘ of E ye. 3n=1,3 
Se) oS pease ; 
I Em ( —1) Mm x NY 
POCOSD Ss. -—— ae cos cos 
oy 1,3 20S 2” e b 
= 6/2 .: 7t z| 
x Ee Sm.n2 692 sinh sn Pay 
(b) Outside the Gap Region (a tg0/2 Sz < ov). 
oo ie) (m +n) /2 
iol ee teint t "7 aah sin“ =* cos ae e~Smnmzla 
A si & m=1,3 n=1,3 WS n,n 
oo ie.) (m+n +1)/ 
atgé Em (— 1) a MAX NEY 
oNera pt Ss AM ee eeL: C08 COS 


x ET smmn 72/8 sinh ee | 


with the following meaning for the various symbols encountered in the formulas: 

(1) » = w/w, with w, = c/a. This means that » is the frequency measured in 
terms of the cut-off value /,. 

(2) Yo= Veo/o 5 

GB) 33, =m? + 1? (a/b)? = 9. 

(4) &m = NEUMANN’S factor, equal to one for m = 0, to two for m + 0. 
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Similar expressions can be derived for the other five fieldcomponents. 

The variation of H,, along the main axis of the guide is depicted in Figure 4 — 
for a few values of the parameters, and shows how much the field penetrates 
into the guide on both sides of the gap region. 


Hy in units of ue sin ot 


D|N 


Figure 4 


H,, along the main axis + = y = 0. The curves are labeled according to the value of tg@. The | 
frequency is equal to half the cut-off value. The aspect ratio a/b of the cross-section is equal to five. 


Appendix 
A relation such as (2.10) finds its origin in the following Maxwell equation. 


OH 
— (4. 1) 


The expansion for H, namely (2.7), can be safely introduced in the second 
member. The first member, however, cannot be obtained by merely ‘curling’ 
equation (2.6) term by term, i. e. summing the curls of the individual terms of 


the expansion for E. It is necessary, on the contrary, to expand curl E sepa- 
rately as: 


curl B= 3 ty Him + 3) On Hin + 3) WnHen + 3 U9 Hy. (4.2) 
m n n 0 


curl B= — p- 
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The coefficient w,, is found to be 

I Flim (x, y) curlE (x, y,2,t) dS 

2 


on PiFimnlenieas ee 
Ss 


This result is obtained by multiplying (4.2) on both sides by H,,,, integrating 
over the cross-section S, and making use of the orthogonality properties of 
the various H,. The next step is to transform the numerator by using the 
following formula of vector analysis: 


curl E = cutle, E, + curl,, E, oe es (ig x ey (4.4) 


where curl,,, indicate that the derivatives with respect to x and y only have 
been kept in the normal formula for the curl. Multiplying both sides of (4.4) 


by Hin yields 


x H,,,) =f H,,, curl ee Y, fee x H.,,) 


H,,, curl E = u,(5* sue 


Be | (4.5) 


su E, curl, Hip, +S CV (Him x E,) : 


If one uses the fact that H,,, = i x Ee Curle, H,,, = — k? E,,,, integrates 
over the cross-section S, and applies the divergence theorem in two dimensions, 
u,, turns out to be 


1 = ao / 

m= ye ae [2 E,-Eim 4S — #2, [* ed se vp E, "Sm de]. (4.6) 
cS 

The first two integrals can be expressed in terms of a,, and c,, (see 2.5), and, 

as a consequence: 


nee Re ae ve, | a (4.7) 


Um Oz Cm 


from which (2.10) follows immediately by equating the coefficients of H,,, in 
both members of (4.1). 
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Résumé 


Le champ électrique dans un guide d’ondes, couplé au milieu extérieur par 
une ouverture dans sa paroi métallique, peut étre calculé en fonction de la 
composante tangentielle du champ dans l’ouverture. Le but de cette note est de 
présenter ce calcul par la methode des modes propres, equivalente a celle du 
tenseur de GREEN, mais douée de certains avantages de souplesse et de clarté. 


(Received: September 30, 1957.) 
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Einschliessung von Eigenwerten und Mohrsches Spannungsdiagramm 


Von ERwIN KREyYSz1IG') 


1. Einleitung 


Es sei A = (a;,) eine n-reihige Hermitesche Matrix. Dann kann man unter 
Verwendung eines beliebigen m-komponentigen Vektors ¥, + 0 und des zugeho- 
rigen Vektors +, = Ax, Aussagen tiber die Eigenwerte 4, = A, = --: = d, dieser 
Matrix [das heisst tiber die L6sungen der Gleichung det (A 6; — aj.) = 0] gewin- 
nen. Es gilt der folgende 

Einschliessungssatz (WIELANDT [1, 2]?)). Jedes Punktepaar auf der reellen 
projektiven y-Geraden, bestehend aus einem beliebig gewahlten Punkt y = u und der 
sugehorigen Nullstelle y = w(u) der lineaven Gleichung 


bevandet zwei (abgeschlossene) Intervalle, deren jedes mindestens einen Eigenwert i, 
dey Matrix A enthailt. 

Diese Einschliessungsintervalle lassen sich geometrisch konstruieren [2]: Man 
bestimmt den Hilfspunkt H in der (y z)-Ebene mit den Koordinaten 


7 ey pre My ee 2 2 

te My’ my ae i (2) 

Der Halbkreis durch H mit dem Mittelpunkt auf der y-Achse und dem einen Rand- 
punkt y= uw, z= 0 hat den Punkt y= w(u), z= 0 als zweiten Randpunkt (Figur a8) 
In Abschnitt 2 der vorliegenden Arbeit wird ein neuer Beweis des obigen Ein- 


schliessungssatzes angegeben und in Abschnitt 3 eine Beziehung untersucht, die 


zwischen dieser Eigenwerteinschliessung und dem Mohrschen Spannungsdia- 
gramm besteht. 


1) Ohio State University, Columbus, Ohio, USA. 
a eee : 2 : ‘ ; 
“) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 206. 


Vol. IXa, 1958 Kurze Mitteilungen — Brief Reports — Communications bréves 203 


Bemerkung 


Bei dem obigen Satz wird die Kenntnis zweier Vektoren %, X,, also eines 
Iterationsschrittes fiir die Eigenwerteinschliessung ausgenutzt (« Einschrittver- 
fahren»). Scharfere Aussagen tiber die Eigenwerte der Hermiteschen Matrix A 
erhalt man, indem man mehr als 2 Vektoren, etwa Poa eNektoren 25.4) = As, 

-,»%,=AxX, 1, 1<p <n, also p Iterationsschritte heranzieht (« p-Schritt- 
verfahven»). Auch hierfiir hat WrELaNpT [3] die bestmégliche Einschliessungs- 
theorie entwickelt. Hinweise zur praktischen Durchfiihrung der Eigenwertein- 
schliessung beim Ein- und Zweischrittverfahren findet man in [4, 5]. 


“| 


> 
w(u) Uu y 


Figur 1 
Zur geometrischen Konstruktion der Einschliessungsintervalle. 


2. Neuer Beweis des Einschliessungssatzes von Wielandt 


Es sei N die Menge der Punktepaare {a, b}, a+ 6, auf der projektiven y-Ge- 
raden, die je zwei (abgeschlossene) Intervalle beranden, deren jedes mindestens 
einen Eigenwert 4, der Hermiteschen Matrix A enthalt. Ein Punkt P der oberen 
(y z)-Halbebene heisse ein ausgezeichneter Punkt, wenn jeder Kreis durch P mit 
dem Mittelpunkt auf der y-Achse diese Achse in zwei Punkten schneidet, die ein 
zu N gehGriges Paar bilden, oder wenn P selbst auf der y-Achse liegt und einem 
der m Koordinatenwerte y = d,, vy = 1, 2,..., m, entspricht. 

Die Menge aller ausgezeichneten Punkte lasst sich leicht durch die (unbe- 
kannten) Eigenwerte der Matrix A kennzeichnen. Offenbar ist ein Punkt genau 
dann ausgezeichnet, wenn er dem abgeschlossenen Bereich B¥ der oberen (y z)- 
Halbebene angehort, der von dem einfachen Halbkreisbogen-n-Eck mit den Spit- 


zen in den Punkten /,,...,/, berandet wird (Figur 2 zeigt B¥). Ein Punkt ist 
demnach genau dann ausgezeichnet, wenn seine Koordinaten y, z den x Un- 
gleichungen 

22+ y2— (Ay tA) yA An S9, | (3) 


Pee A) yt 4420 (v¥=1.2,...0m—1) | 


genugen. : ; 
Auf Grund des Entwicklungssatzes lasst sich der vorstehend genannte Vektor 


%, als Linearkombination der (orthonormierten) Eigenvektoren e,; der Matrix 4, 


also in der Form 
n 
ai = Deve, (4) 


s=1 
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darstellen. Demnach gilt 


und 


Hieraus folgt wegen 4, => 1, = --: = 4, unmittelbar 


n—1 


My — (Ay + Ay) y+ Ay An Mo = ne G (Ay —:) (An — 4) SO, 


= 


n 
4 2 A 2 . 
Mg — (Ay Ay) My + Ay Apps M9 = y G: (Ay — Ax) (Apa — As) 20 
i=1 
7 +v,v4+1 


es Re oe = 1) 


Dies bedeutet, dass die Ungleichungen (3) erfillt sind, wenn man y und z gemass 
(2) wahlt; ein Punkt mit diesen Koordinaten ist ein ausgezeichneter Punkt. Aut 
Grund der Definition der ausgezeichneten Punkte ist damit der Einschliessungs- 
satz bewiesen. 


=A 


As A, A, a, A, y 


Figur 2 


Zum Begriff der ausgezeichneten Punkte. 


3. Beziehung zwischen dem Mohrschen Spannungsdiagramm und der 
vorstehenden Eigenwerteinschliessung 


Das Mohrsche Spannungsdiagramm (Figur 3) gibt einen Uberblick iiber die 
Verteilung der Zug- (bzw. Druck-) Spannungen o und der Schubspannungen r, 
die in einem Punkte Q eines elastischen Kérpers herrschen, wenn dieser Kérper 
ausseren Kraften unterworfen ist. Es sei ¢, der zu einer Flache durch Q mit dem 
Normalenvektor v gehérige Spannungsvektor. Dieser Normalenrichtung ent- 
spricht genau ein Punkt des Mohrschen Diagramms mit der Normalkomponente 
%, von t, als Abszisse und der Tangentialkomponente 1, von #, als Ordinate. 
Beriicksichtigt man jede Normalenrichtung im Punkte Q, so erhalt man genau 
alle Punkte des (abgeschlossenen) Bereichs B in der (o t)-Ebene, der durch das 


il Bile aa 
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Halbkreisbogendreieck mit den Spitzen in den Punkten oj, o, und o, der Abszis- 
_ Senachse berandet wird ; die Spitzen entsprechen den Hauptspannungen (Figur 3) 


5 
3 
= 
he 


: T 
| } 


Figur 3 
Mohrsches Spannungsdiagramm. 


: Die Komponenten des Spannungstensors bilden eine dreireihige symmetrische 

Matrix S, und die drei Hauptspannungen o,, 6, und o, sind die Eigenwerte des 
_ Problems 
SOL 


Der Spannungsvektor ¢,, lasst sich bekanntlich in der Form 


3 
j og Te (7) 


darstellen, wobei der Einheitsvektor e; die zu o; gehorende Hauptspannungsrich- 

tung besitzt und C; der Kosinus des Winkels zwischen dem Vektor e; und dem 

Normalenvektor v ist. Der Vektor ¢, entspricht dem Vektor +, beim Einschlies- 
’ sungsverfahren, und die Darstellung (7) entspricht der Darstellung (5). Fiihrt 
* man vermdoge ¢, = S x, den Hilfsvektor %, ein, so ist wegen (7) 


Daher gilt nun 


3 3 3 
Wig = tq to = x Cr=1, (b) m= *%t,= DG 4: (c) my = tyty = DC g;. (8) 


?=1 1=1 1=1 


Bekanntlich ist (8b) eine Darstellung der Normalkomponente o, des Vektors 1, 
Demnach hat dessen Tangentialkomponente t, die Darstellung 


Ty, = Vt, ty — a. = ms = mi? : (9) 


Wegen (8a) entspricht dies dem in (2) fiir die Koordinate z angegebenen Aus- 
druck. Jeder Punkt des Mohrschen Bereichs B (Figur 3) ist also ein ausgezeich- 
neter Punkt im Sinne der obigen Definition. In der Tat ist der Bereich der aus- 
gezeichneten Punkte bei einer dreireihigen Matrix ein von einem Halbkreisbogen- 
dreieck begrenzter abgeschlossener Bereich, dessen Spitzen in denjenigen Punkten 
der Abszissenachse liegen, die den 3 Eigenwerten der Matrix entsprechen. 
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Summary 


A new proof of a theorem by WIELANDT for inclusion of eigenvalues of Hermi- 
tian matrices is given. Relations between that inclusion method and the Mohr 
stress diagram are investigated. 


(Eingegangen: 8. November 1957.) 


Elementary Solution of Some Plate Problems 


By VActav VopiéxKa, Plzen, Czechoslovakia 


There are many possibilities of generalizing the well-known direct methods of 
solving problems on bending of elliptic plates. The present paper treats from a 
general viewpoint the case where the right-hand side of the basic differential 
equation is represented by a polynomial of any degree. Elementary results cited 
usually in the literature can be obtained as special cases of our more general 
conception. Probably, several of the following formulas are new. 


1. Mathematical statement of the problem. Let 
{(*, ¥) = pL lik x yh (1) 
Osi+ksn 


be a polynomial of any degree n = 0 and let 0/0n denote derivative in the normal 
direction. Then we seek the solution w = w(x, y) of the problem: 


A’w = f(x, y) in the domain b? x? + a yeaa bo OS (2) 
Ow 2 29 9 9 2,2 
w= ie 0 along the contour b° #° + a° y*— a*b?=0; (3) 
herein is 
4 4 4 
ee 0 0 sie 0 


eae 2 2 oe 
ox* Ox” Oy" ov" 


2. Solution of (1) — (3), Assuming 
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and taking into account the elementary fact 


we easily get by substitution of (4) into (2) 


Osi+ksn #) 


2, on!3 (4) bt af t4 yh—4 4 Be [° es 4 ip he: 4 


PCT eral SCL ewes (topo 
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A2(x* y2) = 8 [3 (4) hy 4 ie (5) Pe a me (E\ 4 yet] 


207 


This relation must be satisfied in the whole region b? 72+ a®y?— a?b?<0 


and the comparison of coefficients leads to the system 


ar | eee [6 (ea ae ee te al re. 


ef? hare (5) (ME?) ror 4| 


2 


a 


(eae ae eee 


k+4 


4 274 3 
+3{ 4 )a utes *| 4 ja b C;_oeh +4 


S206" [3 ‘cal a? + ( 4 bea o*] Btu 


—— 
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TRAN Lye eels per : 
—6( 4 ja b Ces pag ee bY 3 4 CiRt4 (5) 


Weviewinswieou 


(j,k =0,1,2,...,2; i +k Sn) 


of equations for the unknowns ¢,,. Of course, in (5) we have to put cap = 0 when- 
ever one of the indexes A, 9 becomes negative or greater than m or when 4+ @ > 7. 

Calculating c,,, from (5) and introducing into (4) leads to the wanted solution 
w(x, vy) of our problem (2), (3). General results are very complicated and so we 
give the complete formulae at least for the case m = 2 which is also of practical « 
significance and can hardly be found in literature. 


3. Special case m = 2. The basic equations (5) become 
1 . 
%1 "90 %a2%19 + %ag%o1 + a4 “20 F Mas 11 + 26 “02 = B, (A=1,2,..., 6); 
p= sa + 2a b4+30*, wo =a,,=0, y= — 268 @ 3b), «0, 


Wig = — 2.0767 (3 a? + 87), BL =foo> 


R 
~) 
= 
I 
i=) 
R 
we 


949 ae 

o9= 3 (a+ 2.4" b° + 5°), Oy = Oyy = M5 =%_e=O0, B= fy, 
4 ee 9. pl 

Os, = %yg= 9, Ayg—=3(Sa°+2a°d +"), Otay = %yzx=%e—e=0, B3=fo,, 


Beg = Og = Cia = Ol, yp 3 (a 4a b° 4 155"), a= 0, 


4 2572 4 
O1 = Oo = AK g=a,,=0, a.=3(5a OG 0°“ 51b )\, ea 0), Bs=hiyp 


» 9 ° 
Oe; =%o—=%g=0, a,=—6a (a+b), a. =0 


> 


65 
gg = S15 a +40"? 4b, Boe 


The solution is quite easy and gives 


1 
8 (3at+2a7b?+3d4) 


00 


2a" b°[(9a° + 41 a* b> +11 a? b*4 3b) f+ (3 a°+4 11440? + 4102 b* +9 0°) f 
3[(a*+ 4a" b+ 154) (1504+ 4a" b?+ bY) — 4 a? b? (a2 + bY)? 


ce 
x io L 


ho 


C, = ea if for 
10 BAe Oh he 4 2 ne ae 4 272 , 
24 (a* +20” b? + 554) 24 (5a*+ 2a®b?+ 54) 


=~" 


. 


—r-lU Ue ee 
\ 


a ee eS 


i 


Vol. IXa, 1958 Kurze Mitteilungen — Brief Reports —- Communications bréves 209 


(15 a* + 4a” b+ 4) f,,— 2b? (a? +B?) fry 


Cy = one oo): Sleaol 2 o.as 2 
Pea 24 (a da’ bP 4.15 64) (15 a*+ 4a? b?+ bY) — 4a? Bb? (a? + B27] 
repels ha 
11 4 2 72 = p4y ” 
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ig 24[(at+4a2b2+ 15 b*) (15 at +4 a2? + bt) 4472? (a? + b2)7] q 


Substituting these values of c,, into (4) leads to the following solution of our 
beginning problem (1) — (3) for the case n = 2: 


1 2 29 >) io 9,92 2 | il 
w(%,y)=—— (b « +a y"-—a db) 5 — > 
(4,9) =>, ( i ) \sa*4 2a 884 ao 
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For a > oo the result (6) reduces to the formula for cylindrical bending 
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whereas for b = a our solution (6) transforms into the formula for a circular plate 
of radius a: 


w(x, 7) = “= (4? + y? — a”)? {4[9 fy +2" (foo + foa)] +12 tio ¥ 
+12 fo) ¥ + (5 fo ~ fox) 46,479 +(— foot 5 fo2) V7} - | 
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If fin = for = fur = 0) fos = Foo then (6. 2)\ gives the formmala 


1 


Pica (a? — 9") (9 fog + 24" foot foo”), = sy y? (6. 3) 


w(r) = 


for symmetrical bending of a circular plate with clamped edge. 


4. Final remarks. (a) If all the f;, = 0, except fo), then our general relations 
(6)—(6.3) lead to the well-known elementary results as given occasionally in 
literature. , 

(b) A detailed study of the basic system (5) allows various interesting conclu- 
sions which are of value in the theory of bending of elliptic plates. 

(c) The above method can be successfully used in solving other boundary 
value problems not only of the biharmonic but also of the polyharmonic and 
harmonic equations!). 


Zusammenfassung 


Die Arbeit behandelt das Biegungsproblem einer am Rande fest eingespannten 
elliptischen Platte unter einer durch allgemeines Polynom u-ten Grades bestimm- 
ten Belastung. In alle Einzelheiten geht die Rechnung fiir den in der Literatur 
kaum zu findenden Spezialfall n = 2. 


(Received: November 25, 1957.) 
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Herbst-Tagung der SPG 1958 


Die nachste Tagung der SPG (Schweizerische Physikalische Gesellschaft) 
wird am 13. und 14, September in Glarus stattfinden (Sekretariat: Basel, Klingel- 
bergstrasse 22). Die SPG wird als Sektion und im Rahmen der 138. Jahresver- 
sammlung der Schweizerischen Naturforschenden Gesellschaft tagen (13.-15. Sep- 
tember, Glarus). K. P. MEYER 


1 r Caer, ~ ca) 
) V. VopiéKa, Elementare Faille des Dirichletschen Problems jiir elliptische Gebiete der Ebene, 
Z. angew. Math. Phys. 8, 309-313 (1957) 
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Progress in Nuclear Energy. Serie 1: Physics and Mathematics, Bd. 1 
(Pergamon Press Ltd., London 1956). 398 S., 195 Fig.; 84s. 

International bekannte Fachleute haben zu diesem ersten Band eines Sammel- 
werkes beigetragen, dem noch weitere sieben iiber Reaktoren, Reaktor-Chemie, 
-Technologie, -Metallurgie sowie biologische, medizinische und é6konomische Aspekte 
der Kernenergie folgen sollen. Die darin enthaltenen Zusammenfassungen gehen 
zum grossen Teil auf die Genfer Atomenergiekonferenz vom August 1955 zuriick, 
an der zum ersten Male ein freier Meinungsaustausch zwischen den Wissenschaftern 
aus 76 Landern iiber alle Probleme und Anwendungen der Kernspaltung stattfand. 
Nicht nur werden die iiber 1000 in Genf vorgetragenen Berichte in der vorliegenden 
Buchreihe gedrangt zusammengefasst; der Inhalt wurde ausserdem von den 24 Mit- 
arbeitern kritisch gesichtet und zu einer in sich geschlossenen Darstellung ver- 
arbeitet. 

Zundchst werden die experimentell bestimmten Daten der Kernspaltungs- 
prozesse, wie Wirkungsquerschnitte, Spaltungsneutronenausbeute und Resonanz- 
parameter, fiir die in Betracht fallenden Uran- und Plutonium-Isotope in Tabellen 
und Kurven dargestellt. Die Theorie dieser Vorgange wird anschliessend von H. A. 
BeETHE behandelt. Drei weitere Artikel befassen sich mit der Streuung und Ab- 
sorption von Neutronen in anderen Materialien als den Kernbrennstoffen, wobei 
auch die Messtechnik und die theoretische Behandlung beriicksichtigt werden. Der 
Erscheinung der verz6gerten Emission von Neutronen ist ein weiteres Kapitel ge- 
widmet. Vier Kapitel tiber die experimentellen Daten und die theoretische Be- 
rechnung von homogenen und heterogenen Reaktor-Anordnungen beschliessen das 
Buch, dessen Studium fiir jeden, der sich mit Planung und Bau von Kernenergie- 
anlagen befasst, eine unumgangliche Notwendigkeit sein wird. D. MAEDER 


Physik diinner Schichten. Teil 2: Stvwktur, elektrische Lettfahighett, magne- 
tische Eigenschaften. Von H. MAYER (Wissenschaftliche Verlagsgesellschaft mbH, 
Stuttgart 1955). 392S., 250 Abb., 28 Tab.; DM 78.-. 

Wahrend im Teil 1 die praktische Herstellung, die Dickenmessung und das 
optische Verhalten diinner Schichten behandelt werden, befasst sich der vorlie- 
gende Teil 2 mit deren strukturellem Aufbau, dem elektrischen Widerstand und 
den magnetischen Eigenschaften. Bei den hier besprochenen diinnen Schichten 
handelt es sich im allgemeinen nicht um einfache Grenzschichten, sondern um 
aus zwei Grenzschichten bestehende Gebilde; beide Begriffe sind durch die in der 
Gréssenordnung der Gitterkonstanten liegende dritte Dimension gekennzeichnet; 
in manchen Fallen lasst sich deshalb unter gewissen Bedingungen das iiber die 
diinne Schicht Gesagte auch auf Grenzschichten (Oberflachen) tibertragen. 

Im ersten, der Struktur gewidmeten Kapitel werden die einzelnen Elementar- 
prozesse beim Schichtaufbau analysiert und sodann ihr Zusammenspiel unter- 
sucht; ferner wird die Strukturanalyse (Réntgenographie, Elektronendiffraktion) 
eingehend behandelt. Die grossen Unterschiede im elektrischen und magnetischen 
Verhalten diinner Schichten gegeniiber dem normalen Festkérper sind Gegen- 
stand des 2. und 3. Kapitels. Nebst den vielen experimentellen Daten wird die 
Festkérper-Elektronik (Leitfahigkeit, Thermokraft, Hall-Effekt, Supraleitfahig- 
keit, Ferromagnetismus) theoretisch dargestellt fiir den speziellen Fall der diinnen 
Schicht, deren Dicke ja mit der freien Weglange der Elektronen vergleichbar ist. 
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In diesem sowohl als Lehrbuch wie auch zum Nachschlagen geeigneten Werke 
werden Theorie und Experiment meist parallel behandelt und miteinander ver- 
glichen und die Giite der verschiedenen Modellvorstellungen diskutiert. Eine 
dominierende Stellung, besonders hinsichtlich der experimentellen Daten, nehmen 
die Metall-Aufdampfschichten ein; weitere Schichttypen sollen in einem zukiinf- 
tigen Band 3 behandelt werden. U. Katz 


American Institute of Physics Handbook. Herausgeber: The American 
Institute of Physics. Coordinating Editor: Dwicut E. Gray (McGraw-Hill Book 
Company, New York, Toronto, and London 1957). 1524 S., 328 Fig., 588 Tab.; 
GAO UZS Old) bos 

Das Buch ist in acht Kapitel gegliedert, die der Reihe nach mathematische 
Hilfsmittel, Mechanik, Akustik, Warme, Elektrizitat und Magnetismus, Optik, 
Atom- und Molekiilphysik und schliesslich Kernphysik behandeln. Die einzelnen 
Kapitel sind weiter unterteilt, so dass insgesamt iiber hundert Teilgebiete der 
Physik umfasst werden. 

Gegeniiber andern Handbiichern ahnlichen Umfangs wirkt das vorliegende 
Buch erstaunlich neuzeitlich, gibt es doch Auskunft iiber neue, ungewohnliche 
oder in der Literatur schwer zugangliche Fachgebiete, wie Geophysik einschliess- 
lich Geodasie, Seismologie, Ozeanographie und Meteorologie; Rheologie; Aero- 
dynamik und Hochdruckphysik. Wir finden Abschnitte tiber Schockwellen und 
auch iiber Wirkungsquerschnitte, Zerfallsprodukte und Mesonen. Der Tieftem- 
peratur-Physiker freut sich an einer Ubersicht iiber die Eigenschaften para- 
magnetischer Salze. 

Jedes der von zustandigen Fachleuten behandelten Teilgebiete wird einge- 
leitet durch erlauternde Erklarungen und Definitionen, die das Lesen der iiber- 
sichtlich angeordneten Tabellen stark erleichtern. 

Man kann dem Buch vielleicht vorwerfen, dass es gelegentlich mehr die Form 
eines Nachschlagewerks fiir Ingenieure als die eines Hilfsmittels fiir Physiker 
aufweist. Trotzdem diirfte es aber der Experimentalphysiker als wertvolle und 
handliche Informationsquelle begriissen. E. MoosER 


Elements of Gasdynamics. Von W. W. LiEPpMANN und A. RosHKo (John 
Wiley & Sons, New York 1957). 439 S., 152 Fig.; $11.00. 

Dieses Buch gibt eine ausgezeichnete Zusammenfassung des heutigen Standes 
der Aerodynamik der kompressiblen Medien, einer Wissenschaft, die sich noch 
in rascher Entwicklung befindet. Zur Einfiihrung dient eine knappe, aber sehr 
klare Darstellung der grundlegenden Begriffe der Thermodynamik. Es folgt eine 
klassische Behandlung der eindimensionalen und ebenen Gasdynamik, mit spe- 
ziellem Hinweis auf die Strémung in Uberschall-Windkandlen. Das Kapitel iiber 
Messmethoden beschrankt sich auf eine knappe Beschreibung der wichtigsten 
Hilfsmittel. Die Abschnitte 7 bis 12 sind den Problemen der 2- und 3dimensionalen 
reibungslosen Strémung gewidmet, wobei die Diskussion sich hauptsachlich auf 
die Theorie der kleinen Stérung im Zusammenhang mit diinnen Profilen und 
schlanken Drehkérpern beschrankt. Der Einfluss der Reibung und der Warmeleit- 
fahigkeit sowie die Ahnlichkeitsgesetze der Unter- und Uberschallstrémung sind 
Gegenstand des nachsten Kapitels. Das Buch tragt den modernen Entwicklungs- 
tendenzen Rechnung, indem auf die Bedeutung der kinetischen Gastheorie und der 
Grundbegriffe der physikalischen Chemie fiir die Gasdynamik hingewiesen wird. 
Eine Sammlung von Ubungen und Zahlentabellen bildet den Schluss des aus- 
gezeichnet gedruckten Buches. P. DE HALLER 
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Applied Probability. Vol..VII. Proceedings of the Seventh Symposium in 
Applied Mathematics of the American Mathematical Society. Gehalten am Poly- 
technic Institute of Brooklyn, 14-15. April 1955. Herausgeber L. A. MacCoLyi 
ae ae Book Company, New York, Toronto, London IOS a WEES: 

s. 6d. 

Beitrage von P. Livy, J. L. Doon, W. Fetter, E. Hopr, G. Mtncu, G. K. 
BATCHELOR, M. Kac, S. M. Utam und B. O. Koopman betreffend Anwendungen 
der modernen Wahrscheinlichkeitsrechnung auf physikalische Probleme, wie 
Brownsche Bewegung, Turbulenz, Quantentheorie, Randwertprobleme usw. Der 
Leser muss die moderne Wahrscheinlichkeitsrechnung gut kennen, wenn er diese 
Vortrage fiihrender Forscher auf diesem Gebiete verstehen will. W. SAXER 


Advances in Applied Mechanics. Band IV. Herausgegeben von H. L. 
DRYDEN und TH. von KArmAN; Herstellungsredaktor G. Kurrti (Academic 
Press Inc., New York 1956). 413 S., 135 Fig.; $10.00. 

Der vierte Band dieser Sammlung enthalt Arbeiten aus dem Gebiet der 
Grenzschichttheorie, der Elastizitat und Plastizitat und der technischen Dyna- 
mik. Francis H. CLAusER gibt eine sehr gute Zusammenfassung des heutigen 
Standes der Theorie der turbulenten Grenzschicht. Unter Annahme einer kon- 
stanten scheinbaren Viskositat im ausseren Teil der Grenzschicht ist es méglich, 
ein gutes Bild der turbulenten Geschwindigkeitsprofile zu geben, unabhangig 
von Druckgradient, Reynoldsscher Zahl oder Wandrauhigkeit. FRANKLIN K. 
Moore gibt einen Uberblick iiber die technisch so wichtige dreidimensionale 
Grenzschicht, mit spezieller Beriicksichtigung der sekundaren Stromungen beim 
Zusammentreffen zweier Grenzschichten. 

Im Kapitel Nonlinear Elasticity von T. C. DoyiEe und J. L. ERIKSEN werden, 
einleitend, einige Begriffe und mathematische Grundlagen dargestellt. Die fol- 
genden Abschnitte behandeln die elastischen Grundgleichungen, Grenzwertpro- 
bleme, iso- und orthotrope Materialien, allgemeine Losungen der Gleichgewichts- 
bedingungen sowie Approximation der Deformationsenergie durch Polynome. 

Der Abschnitt Physical and Statistical Aspects of Fatigue von A. M. FREUDEN- 
THAL und E. GumBEL behandelt auf den ersten Seiten die physikalischen Ur- 
sachen der Ermiidung, welche anschliessend als Grundlage zum Aufbau einer 
statistischen Theorie dienen. Zum Abschluss werden noch einige Fragen der 
Schadenakkumulation diskutiert. 

Im Kapitel Dislocation Theory of Plasticity of Metals beschreibt G. SCHOECK, 
nach einer kurzen Einleitung iiber die plastische Verformung von Kristallen, zuerst 
die Versetzungen in ihrer Geometrie (Stufen- und Schraubenversetzungen, Burgers- 
Vektor, Versetzungsquellen, Stapelfehler usw.) und in ihrem Spannungszustand. 

Ein weiterer Abschnitt behandelt, allerdings sehr kurz, den Nachweis von 
Versetzungen durch Beobachtungen an Kristalloberflachen (Wachstumsspiralen, 
Atzgriibchen an Korngrenzen oder Subkorngrenzen). 

Schliesslich wird ausfiihrlicher auf die mechanischen Eigenschaften von Ein- 
kristallen im elastischen und plastischen Bereich eingegangen, wobei die bekann- 
ten Begriffe der Relaxation, Dampfung, Plastizitat und Verfestigung vom Ge- 
sichtspunkt der Versetzungstheorie aus beleuchtet werden. Dest 

Zum Schluss werden kurz auch einige fiir technische Werkstoffe wichtige 
Phanomena, wie Ausscheidungshartung und das Fliessen von Kristallhaufwerken, 
besprochen. ; 

Der Aufsatz von H. S. Tsten behandelt eine Erweiterung der Poincaréschen 
Methode fiir kleine Stérungen und zeigt, wie diese, sinngemass angewendet, fiir 
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viele Probleme der Gasdynamik gute Dienste leisten kann. Schliesslich sind die 
den Lesern der ZAMP gut bekannten Arbeiten von H. ZIEGLER iiber den Begriff — 
der elastischen Stabilitat in extenso wiedergegeben. Dieses schwierige und tech- 
nisch wichtige Problem ist hier meisterhaft und erschépfend behandelt. 

P. DE HALLER 


Aeroelasticity. Von R. L. BISPLINGHOFF, H. AsHLEy, und R. L. HALFMAN 
(Addison Wesley Publishing Company, Reading 1955). 860 S., 303 Fig;; $ 14.50. 

Nach einem kurzen historischen Uberblick behandelt das Buch in seiner 
ersten Halfte die theoretischen Grundlagen der Aeroelastizitat. Die verschie- 
denen Methoden (unter anderem Lagrangesche Gleichungen, Rayleigh-Ritz- 
Methoden usw.) zur Berechnung der Verformung unter statischer und dynamischer _ 
Belastung, der Eigenfrequenzen und der zugehérigen Schwingungsformen sind 
beschrieben und durch ausfiihrliche Zahlenbeispiele illustriert. Ein grosser Teil 
des Buches ist der Aerodynamik gewidmet: zwei- und dreidimensionale, statio- 
nare und instationdare Stroémungen, Einfluss der Kompressibilitat, werden be- 
handelt. 

Im Abschnitt iiber die eigentlichen Fliigelschwingungen sind einige Dia- 
gramme wiedergegeben, die den Einfluss der verschiedenen Parameter zeigen. 
Die Wirkung der Dampfung, der Kompressibilitat, der Fliigelstreckung usw. 
werden beriicksichtigt. 

Es folgt ein ausfiihrliches Kapitel tiber die Erregung der Schwingungen durch 
Boen und atmospharische Turbulenz, und zuletzt werden die Ahnlichkeitsgesetze 
und die Technik der Modellversuche behandelt. Eine kurze Zusammenfassung der 
notwendigen mathematischen Begriffe ist im Anhang enthalten. 

Dank dem klaren Aufbau des mit vielen Figuren und Zahlenbeispielen ausge- 
statteten Buches und dem ausfiihrliichen Literaturverzeichnis kann dieses Werk 
dem Flugingenieur bestens empfohlen werden. P. DE HALLER 


Linearized Theory of Steady High-Speed Flow. Von G. N. Warp (Cam- 
bridge University Press, London 1955). 243 S., 28 Fig.; 30s. 

Die Schwierigkeit, die nichtlinearen Bewegungsgleichungen der kompres- 
siblen Str6mungen zu lésen, hat zu einer weiten Anwendung der linearen Appro- 
ximationen gefiihrt. Dieser sogenannten linearisierten Theorie und den physi- 
kalischen und mathematischen Annahmen, auf denen sie beruht, ist das vorlie- 
gende Buch gewidmet. 

Der Verfasser wendet sich dabei an eine Leserschaft, welche erstens mit den 
Grundlagen der Gasdynamik und zweitens mit der vektoriellen Schreibweise 
vertraut ist. Letztere Darstellung erlaubt eine physikalisch durchsichtigere 
Schreibweise, erschwert wohl aber dem die orthodoxe Schreibweise gewohnten 
Leser die Lektiire. 

Der erste Teil des Buches enthalt die Herleitungen der linearen Differential- 
gleichungen, fiir stationare Str6mungen, die Lésungen dieser Gleichungen, die Dis- 
kussion der Randbedingungen und die Berechnung der aerodynamischen Krafte. 

Der zweite Teil ist speziellen Randwertproblemen géwidmet, wie der Um- 
stromung diinner Kérper im Unterschallbereich, der Uberschallstrémung um 
nahezu ebene Fliigel, den konischen Uberschall-Stromungsfeldern sowie der so- 
genannten «slender-body »-Theorie. 

Besonderer Dank gebiihrt dem Verfasser fiir das ausfiihrliche Literaturver- 
zeichnis, welches dem Leser gestattet, die entsprechende Fachliteratur zur Er- 
ganzung des oft recht konzentrierten Textes rasch aufzufinden, 
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Das Buch wird als erste zusammenfassende Darstellung der linearisierten 
_ Theorie der Uberschallstrémungen dem fortgeschrittenen Studenten oder werk- 
tatigen Ingenieur wertvolle Dienste leisten. H. R. VoELtLMy 


The Calculation of Atomic Structures. Von Dovucras R. HaRtTREE 
(John Wiley & Sons, New York 1957). 181 S., 14 Fig.; $ 5.00. 

Das vorliegende Buch gibt von berufenster Seite eine griindliche Darstellung 
der Methode des « Self-Consistent Field» zur Bestimmung der Atomstruktur. 
Das Hauptgewicht liegt dabei auf der tatsachlichen Durchfiihrung der Methode 
bis zur numerischen Behandlung. Die Darstellung ist bewusst elementar gehalten. 
Es werden zum Beispiel keine gruppentheoretischen Methoden verwendet. Aus- 
gangspunkt bildet das Variationsprinzip zur Schrédinger-Gleichung mit einem 
Determinantenansatz fiir die Testfunktion. Die Focksche Integro-Differential- 
gleichung wird durch sukzessive Approximation gelést, wobei die Wahl der ersten 
Naherung mit besonderer Sorgfalt diskutiert wird. Uberall kommt die grosse 
Erfahrung des Autors auf diesem Gebiet zur Geltung. Ein besonderer Abschnitt 
ist der numerischen Behandlung von Differentialgleichungen gewidmet. Der 
Schluss des Buches diskutiert relativistische Verfeinerungen und Approximatio- 
nen, die tiber die Methoden der « Self-Consistent-Fields » hinausgehen. Zwei An- 
_hange geben Literaturhinweise auf Resultate und reduzierte radiale Wellen- 

funktionen. 
; Die Tatsache, dass das Buch aus Vorlesungen vor einem breiteren Auditorium 
entstanden ist, ist angenehm spiirbar und erleichtert dem Aussenstehenden die 
Lektiire. R. Jost 


Principles and Techniques of Applied Mathematics. Von B. FRIEDMAN 
(John Wiley and Sons, New York; Chapman & Hall, London 1956). 315 S., 6 Abb.; 
$8.00. 

In einem Werk mit diesem Titel erwartet der Leser eine Behandlung der 
theoretischen Grundlagen, auf denen die numerischen Methoden der angewandten 
Mathematik letzten Endes beruhen. Es beginnt denn auch vielversprechend mit 
der Theorie der linearen Operatoren, welche recht anschaulich und mit zahl- 
reichen Ubungsbeispielen behandelt wird. Es ist durchaus zu begriissen, dass der 
Verfasser unter Hinweis auf die entsprechenden Literaturstellen oft auf die 
exakten Beweise verzichtet; dass er aber in Kapitel 1 und 2 auf die fiir die 
- angewandte Mathematik so wichtigen Extremalprinzipien tiberhaupt nicht ein- 
geht, wahrend andererseits die nichtlinearen Elementarteiler in unndtiger Aus- 
fiihrlichkeit behandelt werden, wird doch als Mangel empfunden. 

Im Kapitel 3 wird die Greensche Funktion unter Zuhilfenahme der Dirac- 
schen Impulsfunktion 6(*) eingefiihrt, was dem iiblichen Vorgehen entschieden 
vorzuziehen ist. Freilich ist das ganze Kapitel stark auf Differentialoperatoren 
L = — (d/dx)? + g(x) zugeschnitten, ebenso das folgende Kapitel iiber Eigen- 
wertprobleme bei gewohnlichen Differentialgleichungen. In diesem Kapitel 
kommt auch die Spektraldarstellung von Operatoren zur Sprache ; als spezielle 
Anwendung davon erhalt man Darstellungen von 6(¥) durch divergente Summen 


und Integrale wie 
+ CoO 


220(x) = / exp(ik x) dk, 
—0oo 


die hier allerdings nicht exakt begriindet werden. 
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Die Betrachtung iiber die Greensche Funktion fiir die genannten Differential- 


operatoren liefern unmittelbar auch die Greensche Funktion fiir Potentialpro- — 


bleme in Rechteckbereichen sowie fiir die Wellengleichung und gestatten, die 


myp eee!) 


bekannten expliziten Lésungen sofort anzugeben. _ H. RuTISHAUSER - 


Technische Hydro- und Aeromechanik. Von WALTER KAUFMANN © 


_ (Springer-Verlag, Berlin 1954). 352 S., 254 Abb.; DM 36.-. 


Das vorliegende Buch enthalt in grossen Ziigen den Stoff der Vorlesung, . 
welche der Verfasser an der Technischen Hochschule Miinchen halt. Entsprechend ~ 


der Absicht des Verfassers, ein Lehrbuch fiir Studierende zu schreiben, zeichnet 
sich auch diese Neubearbeitung des 1931 und 1934 erschienenen Buches Ange- 
wandte Hydromechanik durch einen klaren, methodischen Aufbau und durch 
meist ausfiihrlichen Text aus. Besonders alle mit der Grenzschichttheorie zusam- 
menhangenden Fragen werden griindlich erlautert. Die einzelnen Abschnitte des 
Buches tragen folgende Titel: Eigenschaften der Fliissigkeiten und Gase, Gleich- 
gewicht (Hydro- bzw. Aerostatik), Bewegung der Fliissigkeit (Hydro- bzw. Aero- 
dynamik) mit den Unterabschnitten: Eindimensionale Str6mung (Reibungsfreie 


Stromung, StrOmung mit Energieverlusten), Allgemeine Theorie der zwei- und © 


dreidimensionalen Strémung, Grundlagen der Dynamik kompressibler Fliissig- 
keiten (Gasdynamik). 

Ein Sachverzeichnis erleichtert das rasche Auffinden des Gesuchten. Die vie- 
len im Text eingestreuten klaren Diagramme und Zeichnungen zusammen mit 
dem sauberen Druck des Verlages ergeben ein Lehrbuch, das nicht nur dem 
Studierenden, sondern auch dem praktisch arbeitenden Ingenieur als Nach- 


schlagewerk warm empfohlen werden kann. H. R. VOELLMY © 


Water Waves. Von J. J. STOKER ([Pure and Applied Mathematics, Bd. 4] 
Interscience Publishers, New York and London 1957). 567 S., 189 Abb.; $12.00. 

Obwohl zwischen den Gravitationswellen an der Oberflache einer Fliissigkeit 
und den Wellen der Gasdynamik vielfach eine enge Analogie besteht, sind die 
Oberflachenwellen weniger spektakular. Viele hierher gehérende Probleme sind 
erst in neuerer Zeit von der hydraulischen Ebene auf das Niveau der exakten 
Hydrodynamik gehoben und einer strengen Lésung zugefiihrt worden, wobei ins- 
besondere die New-Yorker Schule mit STOKER, FRIEDRICHS und anderen wesent- 
liche Beitrage geleistet hat. Um so erfreulicher ist es, dass es STOKER, ein Autor, 
der ebensosehr fiir seine mathematische Strenge wie fiir seine Aufgeschlossenheit 
physikalischen und technischen Problemen gegeniiber bekannt ist, unterhnommen 
hat, das Gebiet umfassend darzustellen. 

Das Werk, das uneingeschrankte Anerkennung verdient, behandelt in vier 
Hauptteilen: 

iS Die Grundlagen der Dynamik reibungsloser, inkompressibler Fliissigkeiten 
samt einer exakten Begriindung der Naherungstheorien fiir kleine Amplituden 
bzw. geringe Wassertiefe ; 

II. Lésungen der Theorie kleiner Amplituden (Wellen auf geneigtem Unter- 
grund und in der Nahe von Hindernissen, Wellen auf einem Strom, Bewegung von 
Schiffen im Seegang) ; i ae 

III. Lésungen der Seichtwassertheorie (Fortpflanzung von Stérungen, Bre- 
chen von Wellen, schwimmende Wellenbrecher, Unstetigkeitsfronten in der 
Atmosphare, Bewegung in Kanalen mit rauhen Wanden, Flutwellen in Stré6men) ; 

IV. Lésungen der exakten Theorie (Dammbruch). H. ZIEGLER 
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der Laplace- Transformation 


Ein tc bas fiir Studierende der Mathematik, _ Physik und 
Tngenieurwissenschatt. FE 


Von Prof. Dr. Gustav Doetscu, Professor an der Universitat Freiburg i.B. 


. 


-Mathematische Reihe, Band 24 — Sammlung «Lehrbiicher und Monographien ‘ 


aus dem Gebiete der exakten Wissenschaften» 
304 Seiten mit 40 Figuren. (1958). Ganzleinen Fr. 39.40 (DM 39.40). 


Inhalt (gekiirzt): Das Laplace-Integral von physikalischen und mathema- 
tischen Gesichtspunkten aus — Die Konvergenzhalbebene und die Holomor- 
phie der Laplace-Transformierten — Die Abbildung der fundamentalen Opera- 
tionen — Gew6éhnliche Differentialgleichungen und Systeme von solchen — 
Differenzengleichungen — Das Verhalten der Laplace-Transformierten im 
Unendlichen — Das komplexe Umkehrintegral, Deformation des Weges, Aus- 
wertung des Integrals — Bestimmung der Originalfunktion durch Reihenent- 
wicklung— Die Parsevalsche Gleichung — Asymptotische Entwicklung der Bild- 
und der Originalfunktion — Gewdhnliche Differentialgleichungen mit Poly- 
nomkoeffizienten — Partielle Differentialgleichungen — Integralgleichungen. 


Das dreibandige «Handbuch der Laplace-Transformation» ‘des Verfassers 
stellt das gesamte heute vorliegende Material iiber die Laplace-Transforma- 
tion zusammen und ist in erster Linie als Grundlage fiir die wissenschaftliche 
Arbeit auf diesem Gebiet anzusehen. Im Gegensatz dazu bringt das vorlie- 
gende kurze Lehrbuch nur den Stoff, der fiir den «Normalverbraucher» un- 
erlasslich ist, wenn er die Literatur, die sich der Laplace-Transformation 
bedient, wirklich verstehen oder die Laplace-Transformation bei seinen 
eigenen Arbeiten benutzen will. Wahrend das kirzlich erschienene Buch 
des Verfassers «Anleitung zum praktischen Gebrauch der Laplace-Trans- 
formation» nur die Technik des Kalkiils ohne Beweise der benutzten Satze 
lehrt, behandelt das vorliegende Werk die Laplace-Transformation lehrbuch- 
mAssig in voller Strenge und mit ausfiihrlichen Beweisen 


Zu beziehen durch Ihre Buchhandlung — Obtainable from your bookseller 
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